Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,423,225 G→A 91.8% pseudogene (25/444 nt) rzpR → pseudogene, Rac prophage; Bacteriophage Rz lysis protein family;Phage or Prophage Related; putative Rac prophage endopeptidase; completely contained in another CDS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,423,2250GA91.8% 758.7 / 62.9 566pseudogene (25/444 nt)rzpRpseudogene, Rac prophage; Bacteriophage Rz lysis protein family;Phage or Prophage Related; putative Rac prophage endopeptidase; completely contained in another CDS
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base A (318/197);  total (345/221)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 640 total aligned reads displayed.

GTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCT  >  NC_000913/1423079‑1423370
                                                                                                                                                  |                                                                                                                                                 
    cTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCt                                                                                                                                             >  3:115096/1‑149 (MQ=255)
        aaCGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCt                                                                                                                                         >  8:247928/1‑149 (MQ=255)
          cGCAGGCAATGGCAATCAGAAAAGGTCGTCTGATTAGCGTAAATCAGGGCGGTAAATTACGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGCTATCTCCGCTCTGGTTAACTGCATCGTCGTCAGCCAGt                                                                                                                                       >  2:40086/1‑149 (MQ=255)
           gCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGt                                                                                                                                       >  4:329352/1‑148 (MQ=255)
           gCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTc                                                                                                                                      <  3:142592/149‑1 (MQ=255)
                   tggcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCAt                                                                                                                                    <  5:190938/143‑1 (MQ=255)
                                        tGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACg                                                                                                         >  3:336071/1‑149 (MQ=255)
                                                            ggTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGgt                                                                                                                               >  3:168224/1‑107 (MQ=255)
                                                                               tACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCa                                                                     >  4:18097/1‑146 (MQ=255)
                                                                                                                  ccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAACAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  5:35270/8‑149 (MQ=11)
                                                                                                                  ccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  8:78231/8‑149 (MQ=11)
                                                                                                                  ccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  5:366946/8‑149 (MQ=11)
                                                                                                                  ccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  1:203707/8‑149 (MQ=11)
                                                                                                                   cgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGACCTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  2:248582/7‑148 (MQ=11)
                                                                                                                   cgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACt                                     <  7:366211/136‑1 (MQ=11)
                                                                                                                   cgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgc                              <  8:338259/143‑1 (MQ=11)
                                                                                                                   cgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  8:190671/7‑148 (MQ=11)
                                                                                                                    gcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGCGAACTGAAGCTGGCGACCGCGGCAATTACTGacatg                               >  5:301374/6‑147 (MQ=11)
                                                                                                                    gcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgca                             <  2:22153/144‑1 (MQ=11)
                                                                                                                    gcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  5:175785/6‑147 (MQ=11)
                                                                                                                    gcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  7:199897/6‑147 (MQ=11)
                                                                                                                      gattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  7:169055/4‑145 (MQ=11)
                                                                                                                      gattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  1:383493/4‑145 (MQ=11)
                                                                                                                      gattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                               >  6:202488/4‑145 (MQ=11)
                                                                                                                       attATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagat                          <  3:357350/147‑1 (MQ=11)
                                                                                                                       attATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagat                          <  7:64414/147‑1 (MQ=11)
                                                                                                                          aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                        <  2:9790/148‑1 (MQ=11)
                                                                                                                          aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                        <  1:62272/148‑1 (MQ=11)
                                                                                                                          aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                        <  3:188675/148‑1 (MQ=11)
                                                                                                                          aTCTCCGCTCTGGTTATCTGCATCATCGGCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAACCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                        <  6:357120/148‑1 (MQ=11)
                                                                                                                          aTCTCCGCGCTGGTTATCTGCATCATCGTCTGCCGGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                        <  6:97297/148‑1 (MQ=11)
                                                                                                                           tctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatg                         >  6:178912/1‑146 (MQ=11)
                                                                                                                           tctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      >  8:217569/1‑149 (MQ=11)
                                                                                                                           tctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      >  2:60048/1‑149 (MQ=11)
                                                                                                                           tctcCGCTCTGGTTATCTGCATCATCGTCTGACTGTCATGGGCTGTTAATCATTACCGTGATAACGACATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGAAGATGCg                       >  5:333802/1‑148 (MQ=11)
                                                                                                                            ctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      >  5:92621/1‑148 (MQ=11)
                                                                                                                            ctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      >  1:175083/1‑148 (MQ=11)
                                                                                                                            ctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACAGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      >  5:149702/1‑148 (MQ=11)
                                                                                                                             tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGaa                                                           <  6:288106/110‑1 (MQ=11)
                                                                                                                             tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGGaa                                          >  6:37500/1‑127 (MQ=11)
                                                                                                                             tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      >  1:283063/1‑147 (MQ=11)
                                                                                                                             tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      >  1:304596/1‑147 (MQ=11)
                                                                                                                             tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      >  6:80904/1‑147 (MQ=11)
                                                                                                                              ccGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCa                                           <  3:273246/125‑1 (MQ=11)
                                                                                                                              ccGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAg                   >  3:98499/1‑149 (MQ=11)
                                                                                                                               cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATACCGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCGGATGCGTCAGc                  >  4:119446/1‑149 (MQ=11)
                                                                                                                               cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                      <  6:236505/145‑1 (MQ=11)
                                                                                                                               cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTc                     <  6:162402/146‑1 (MQ=11)
                                                                                                                               cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGc                  <  1:256696/149‑1 (MQ=11)
                                                                                                                               cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGc                  >  4:147169/1‑149 (MQ=11)
                                                                                                                               cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGc                  <  3:82908/149‑1 (MQ=11)
                                                                                                                                 ctctGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGt                >  4:26552/1‑149 (MQ=255)
                                                                                                                                 ctctGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGt                >  7:25355/1‑149 (MQ=11)
                                                                                                                                  tctGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGc                  >  7:319328/1‑146 (MQ=11)
                                                                                                                                  tctGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGt                >  8:363982/1‑148 (MQ=11)
                                                                                                                                    tGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTg               <  1:33154/147‑1 (MQ=11)
                                                                                                                                      gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAACATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt           <  2:326804/149‑1 (MQ=11)
                                                                                                                                      gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                        <  2:310040/136‑1 (MQ=11)
                                                                                                                                      gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATg            <  3:158313/148‑1 (MQ=11)
                                                                                                                                      gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATg            <  6:128190/148‑1 (MQ=11)
                                                                                                                                      gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACAGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt           >  8:375658/1‑149 (MQ=11)
                                                                                                                                      gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAAGCATTACCGTGATAACGCCATGACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGACCGCGGCAATTACTGCCATGCAGATGCGTCAGCGTGATGt           >  2:184144/1‑149 (MQ=11)
                                                                                                                                       ttATCTGCATCGTCGTCTGCCTGTCATGGGGTGTGAATCATTACCGTGATAACGCCATCGCCTACACAGAACAGCGCGACATAAATGCCAGAGAACTGCAGCTGGCGAACGCGGCAATTAATGAGATGCAGATGCGTCAGCGTGATGtt          >  4:207655/1‑149 (MQ=17)
                                                                                                                                       ttATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGtt          >  4:163256/1‑149 (MQ=11)
                                                                                                                                       ttATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGtt          <  5:376535/149‑1 (MQ=11)
                                                                                                                                         aTCTGCATCGTCGTCTGCCTGTAATGGGGTGTTAATCATTACCGTTATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTtgc        >  8:101421/1‑149 (MQ=255)
                                                                                                                                         aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTGATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCAAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGAtgc        >  1:106056/1‑149 (MQ=11)
                                                                                                                                         aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgc        >  8:13792/1‑149 (MQ=11)
                                                                                                                                         aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgc        >  1:8328/1‑149 (MQ=11)
                                                                                                                                         aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgc        >  3:227714/1‑149 (MQ=11)
                                                                                                                                         aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtg         <  7:325450/148‑1 (MQ=11)
                                                                                                                                          tCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTa                                       >  8:29660/1‑117 (MQ=11)
                                                                                                                                          tCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagat                          <  4:300550/130‑1 (MQ=11)
                                                                                                                                           cTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgctg      >  6:250448/1‑149 (MQ=11)
                                                                                                                                           cTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct       <  8:15499/148‑1 (MQ=11)
                                                                                                                                            tGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATg                                                                        <  2:358428/82‑1 (MQ=11)
                                                                                                                                            tGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATg                                                                        >  1:358428/1‑82 (MQ=11)
                                                                                                                                            tGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgctgc     >  5:281766/1‑149 (MQ=11)
                                                                                                                                            tGCATCATCGTCTGCCTCTCATGGGCTGTTAATCATTACCGAGATAACGCCATTACCTACAACGCCCCGCGCGCCCAAAATGCCAGAGAACTGAAGCAGGCGCACGCGGCAATTACTGACATTCTGATGCGACAGCGTGATGTtgctgc     >  2:104587/1‑149 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCCGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcg    >  1:31594/1‑149 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGc                                                                       <  3:55253/82‑1 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                     >  7:13958/1‑84 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTgg                                                      <  2:276561/99‑1 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGc                                                     >  6:323356/1‑100 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGc                                            <  4:367592/109‑1 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcg    >  2:1412/1‑149 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcg    >  6:62185/1‑149 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcg    >  7:283988/1‑149 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcg    >  5:1264/1‑149 (MQ=11)
                                                                                                                                             gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcg    >  4:83321/1‑149 (MQ=11)
                                                                                                                                              catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTTATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgctgc     <  6:10148/147‑1 (MQ=11)
                                                                                                                                              catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgc                              <  3:295792/122‑1 (MQ=11)
                                                                                                                                              catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgctgc     <  4:229099/147‑1 (MQ=11)
                                                                                                                                              catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct       <  7:44678/145‑1 (MQ=11)
                                                                                                                                              catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct       <  8:389946/145‑1 (MQ=11)
                                                                                                                                              catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcgc   >  3:316226/1‑149 (MQ=11)
                                                                                                                                              catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcgc   <  5:264835/149‑1 (MQ=11)
                                                                                                                                               atcatcGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCt  >  8:69216/1‑149 (MQ=11)
                                                                                                                                               atcatcGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCt  >  7:247441/1‑149 (MQ=11)
                                                                                                                                               atcatcGTCTGCCTGTCATGGGCTGTTAATAATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGACCAGAAGCTGGCGCACCCGCCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCt  >  8:320787/1‑149 (MQ=11)
                                                                                                                                                  |                                                                                                                                                 
GTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCT  >  NC_000913/1423079‑1423370

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: