Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,173 C→A 7.0% intergenic (‑107/+488) lacZ ← / ← lacI pseudogene, truncated/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,1730CA 7.0% 488.7 / 8.0 178intergenic (‑107/+488)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (45/115);  new base A (9/9);  total (54/124)
Fisher's exact test for biased strand distribution p-value = 6.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

GCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGA  >  CP009273/360027‑360321
                                                                                                                                                  |                                                                                                                                                    
 cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                                   <  2:1925619/149‑1 (MQ=255)
 cGGCTCACCATCCCGCGCCACCAGCCAGTGCTGGCGACCCGTTGATACGCTGTAAAACATCATAAGTCCGACCCCAATTTCCTGGCTGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                                   <  1:3797355/149‑1 (MQ=17)
 accgacaaacaacagatacaacgaaaggcccagtctttcgcatgagcctgtcgttttatttgatgcctggggaatcAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                                   >  2:3865899‑M1/77‑149 (MQ=255)
  ccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgccgggggaattATTTCCTGGCAGGTTATGGCGGGCGTCCTGCGCGCCACCCTCCGGGCCGTGGGGTGGCACCGTTAAAATCCgcg                                                                                                                                                  <  2:1496339‑M1/74‑2 (MQ=255)
    gacaaacaacagataaaacgaaaggcccagtctttcgactgagccttgcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACCTGCAAATCCCCTcc                                                                                                                                                <  1:2215705‑M1/76‑1 (MQ=255)
     ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                               >  2:707465/1‑149 (MQ=255)
     taacaaTCCAGCCCCACCATCCATTGCAGGCGCTCCGATTATATGATTTAAAACATGATAATTCCGACCCCCATTTCCTGGATGTTTTTGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                               <  1:1991423/144‑1 (MQ=18)
      caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCGCGTTCAAagccgcggccg                                                                                                                                              <  2:2574738‑M1/78‑6 (MQ=255)
      caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAAACCGCTCCcg                                                                                                                                              <  2:2539442‑M1/78‑1 (MQ=255)
          caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTCAATTCCGCTCCcggcgg                                                                                                                                          <  2:4169247‑M1/82‑1 (MQ=255)
           aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCACCCCGCGGa                                                                                                                                         <  1:2470500‑M1/83‑1 (MQ=255)
            ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                        <  2:918077/149‑1 (MQ=255)
            acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcccggcggct                                                                                                                                        <  1:651853‑M1/84‑3 (MQ=255)
            acagataaaacgaaaggcccagtctttcgactgagcctttcgtggtatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGgag                                                                                                                                        <  2:3228445‑M1/84‑2 (MQ=255)
              aGCGCCACCATCCAGTGCAGTAGCCCCGATTATAAGCTGTCAAAAATCATAATTCCCATCCCCATTTCCTGGCTTTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                      <  1:3981567/149‑1 (MQ=18)
               gcgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCCCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                     >  2:1275968/1‑149 (MQ=255)
               gataacccgaacggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                     >  2:3125557‑M1/63‑149 (MQ=255)
                 gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGATAATTCCGACCCCCAAGTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                   <  1:2406318/149‑1 (MQ=255)
                 gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGGTCGCAACCTTCAAATCCCCGCCCCGCGGCttttta                                                                                                                                   >  1:4093574/1‑146 (MQ=25)
                  aaaacgaaaggcccagtctttcgactgagcctttcgttgtatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTACTGCCCGCCACCCTCCGGGCCGTTGGTTAGCCACGTTGAAATACGCTAccggcggaagtttac                                                                                                                                  <  1:2044293‑M1/90‑8 (MQ=255)
                    accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                >  2:3365114/1‑149 (MQ=255)
                    aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTACCGGCGGAGTTGTCCTa                                                                                                                                <  2:321232‑M1/92‑1 (MQ=255)
                       aTCCAGTGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                                             <  1:3228888/149‑1 (MQ=255)
                        tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                            >  2:3562388/1‑149 (MQ=255)
                         ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                           <  2:1994886/149‑1 (MQ=255)
                           ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                                         <  2:2790113‑M1/99‑1 (MQ=255)
                               cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                     <  2:2863054/149‑1 (MQ=255)
                                 ggAGCTCGGATGATATGCTGTCAAACATGATAATTCCGATCCCCAATTCGTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                   <  1:1678907/149‑1 (MQ=18)
                                  tatttcgattgagcctttctttttatgtgatgccggtggaattAATTCCGGGCAGTTTTTGGGGGGCGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                  >  1:289283‑M1/44‑149 (MQ=255)
                                  gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                  <  1:1772533/149‑1 (MQ=255)
                                   ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                 <  2:2115110‑M1/107‑1 (MQ=255)
                                    tatcttatgatACGCTGTCAAACTTTCTAATACCGTTCCCAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                                <  2:1605051/143‑1 (MQ=14)
                                     ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                               <  2:3512162‑M1/109‑1 (MQ=255)
                                      tcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCTTAAAATAAGCACCCCGCGGATTTGTACGACACACGATCGCGTTCAcc                                                                                                              <  1:4179074‑M1/110‑1 (MQ=255)
                                      tCCGCTGATAAGCTGTCACACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                              <  2:998227/149‑1 (MQ=255)
                                      tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                              <  2:1756187/149‑1 (MQ=255)
                                      tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCCGCGGCTTTGTCCGACTCAGGAGAGCGTTCAcc                                                                                                              >  1:3908196/1‑149 (MQ=255)
                                          atgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAAACCCACCCCGCGGCATTGTACTAAACAGGGGTGCGTTAACCGaaa                                                                                                          >  1:1605051/1‑147 (MQ=18)
                                           tgagcctttcgttttatttgatgcctggggaattAATTCATGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                         <  2:1773922‑M1/115‑1 (MQ=255)
                                             atattttctttttatttgatgcctgggtaattAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                       >  2:92982‑M1/33‑149 (MQ=255)
                                             atAAGCTGTCAAACATGTTATGTCCGAGCCCCAATTCCTGGCAGTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                       <  2:815417/149‑1 (MQ=14)
                                             atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                       <  2:1656225/149‑1 (MQ=255)
                                              tatgCTGTCAAACATGAGATTTCCGATCCCCTTTTCCTGGCTGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCATCGacaaaca                                                                                                      <  1:163368/146‑1 (MQ=17)
                                                aGCTTTCAAAATTGACAATTCCGATCCCCTAGTCCTGGATGTTTAGGGCGGGGGTCGTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGGGAGCGTTCACCGACAaacaac                                                                                                    <  1:1564547/149‑1 (MQ=14)
                                                 tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                   <  2:3806454‑M1/121‑1 (MQ=255)
                                                 tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTAAAATCCGCACCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACcgaccaacaaaa                                                                                                   <  2:1138657‑M1/121‑3 (MQ=255)
                                                   tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCCACCTTCAAATACCCACCCGGCGGGTTTGTACTACGCAGGAGAGCGGTCACCCACCAacaacaca                                                                                                 >  2:2527781/1‑147 (MQ=17)
                                                     catagttttgatgcccggtgatttAATTCCGGGCTGTGTATGGCGGGCGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGatc                                                                                               >  1:992957‑M1/25‑148 (MQ=255)
                                                      cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                              >  2:342918/1‑149 (MQ=255)
                                                      cAAACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                              <  2:3375951/149‑1 (MQ=255)
                                                        ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                            >  2:3635160‑M1/22‑149 (MQ=255)
                                                        aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                            <  1:276446/149‑1 (MQ=255)
                                                         aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCCAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAACCACCAGATAAAAc                                                                                           >  2:3014651/1‑149 (MQ=17)
                                                            tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                                        <  2:3957263/149‑1 (MQ=255)
                                                             gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGCAAAGTTAAAAAACGCTCCCCGCGGGATTTTCCAACAAAGGAGAACGTTAACCGACAAACAAAACATaaaaaaaaa                                                                                       >  1:2073859/1‑144 (MQ=9)
                                                              agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                      <  2:2622703/149‑1 (MQ=255)
                                                                tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACCTTCCAATCCCCGCCCGGCGGATTTGTCCTACTCAGGAGAGCGGTCACCCACAAACAACAGAAAAAaccaaacgc                                                                                    <  2:1248486‑M1/136‑4 (MQ=255)
                                                                 aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCTTAAAAAACGCACCCCGCGGATTTGTACAAAACCGGAGCGGGTTCACCCCCAAAAAAAAGATAAAAAGAAAAGcc                                                                                   >  1:3029238/1‑149 (MQ=9)
                                                                  ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGccc                                                                                  <  2:2153542/149‑1 (MQ=255)
                                                                  gtCCGATCACAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                  <  2:1936796/148‑1 (MQ=255)
                                                                  cctgggtatttAATTCCGGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                  >  2:811453‑M1/12‑149 (MQ=255)
                                                                    ccgtgcaccAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                                <  2:1652869/142‑1 (MQ=12)
                                                                      gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                              >  2:1341584‑M1/8‑149 (MQ=255)
                                                                       ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACTTCAAATCCCCTCCCGGCGGCTTTGTCCTACACAGGCGAGCGTTCACCGCCAAACAACAGAAAAAACGAAAGGCCCAGTCt                                                                             <  2:1208690‑M1/143‑1 (MQ=255)
                                                                       aTCCCGAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                             <  2:3286893/149‑1 (MQ=255)
                                                                        tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                            <  2:1168469/149‑1 (MQ=37)
                                                                          cccATTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                          <  2:1882193/149‑1 (MQ=17)
                                                                             aaTTCCTGGATGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                       <  1:942665/149‑1 (MQ=17)
                                                                              aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTAACCGCCAAACAACAGATAAAAAGAAATGCCCAGTCTTTCGACt                                                                      >  2:537936/1‑149 (MQ=11)
                                                                               ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                     <  2:2727904/149‑1 (MQ=31)
                                                                               ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                     <  2:2727966/149‑1 (MQ=31)
                                                                               ttCCGGGCTGTTTTGGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                     <  1:4050746/149‑1 (MQ=11)
                                                                                     gctgTTTTTGGGGGGCGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCtt                                                               <  1:2536455/146‑1 (MQ=1)
                                                                                          ttggggggggcgGCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGtt                                                          <  1:731675/142‑1 (MQ=1)
                                                                                                    cgcccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATGTGTCCCACTCACGAGAGCGTTCACCGACAAACAACAGATACACCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAATTGATg                                                <  2:2905856/146‑1 (MQ=1)
                                                                                                       ccggCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCt                                             <  1:469885/146‑1 (MQ=1)
                                                                                                           cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                         <  1:3806454/149‑1 (MQ=31)
                                                                                                               ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAACGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAACCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                     <  1:1299132/149‑1 (MQ=255)
                                                                                                               ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGACAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                     <  1:1275968/149‑1 (MQ=21)
                                                                                                               ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                     <  1:3224896/149‑1 (MQ=34)
                                                                                                                cACCCTCCGGGCCGTTGCTTCGCAACGGTAAAAAACGCTCCCCGGGGATTTGTCCTACACACGAGAGCGTTAACCGACAAACAAAAAATAAAACGAAATGCCCAGTAATTCGACTGTGCCTTTCGTTTTATTTTTTTCCTGgggaaatt                                    >  2:1795884/1‑146 (MQ=11)
                                                                                                                   ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggcattggga                                 >  2:3225048/1‑149 (MQ=25)
                                                                                                                   ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga                                 >  2:3224185/1‑149 (MQ=38)
                                                                                                                   ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga                                 >  2:1043612/1‑149 (MQ=38)
                                                                                                                    cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaa                                <  1:1051962/149‑1 (MQ=38)
                                                                                                                     tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat                               <  2:2428601/149‑1 (MQ=39)
                                                                                                                          gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt                          <  1:1664825/149‑1 (MQ=255)
                                                                                                                            cGTTGCTTCGCAACGTACAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTACACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt                        <  1:13575/149‑1 (MQ=255)
                                                                                                                                  ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGa                  >  1:3212741/1‑149 (MQ=255)
                                                                                                                                  ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa                  <  1:1334339/149‑1 (MQ=255)
                                                                                                                                     gCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTGTTATTTGATGCCTGGGGACTTGGGAAGTAAGTCTTGAAGACGa               >  2:2457557/1‑149 (MQ=255)
                                                                                                                                         cGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTATCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAAGTGGGAATTAATTCTTGAAGACGAAAgg           >  2:181715/1‑149 (MQ=255)
                                                                                                                                         cGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGTTTTATTTGATGCCTGGGGAATTGGGACTTATTTCTTGAAGACGAAAgg           >  2:1279796/1‑149 (MQ=17)
                                                                                                                                           ttCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc         <  2:845102/149‑1 (MQ=255)
                                                                                                                                           ttCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc         <  2:845655/149‑1 (MQ=255)
                                                                                                                                              aaaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTc      >  1:970894/1‑149 (MQ=255)
                                                                                                                                               aaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCg     <  1:3179228/149‑1 (MQ=255)
                                                                                                                                                 tCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg   >  2:3570698/1‑149 (MQ=255)
                                                                                                                                                  ccGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGCGCCTTTCGTTTTATTTGTTGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGa  <  1:3367363/149‑1 (MQ=255)
                                                                                                                                                  ccGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGa  <  1:1043612/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGA  >  CP009273/360027‑360321

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: