Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,157 T→G 5.8% intergenic (‑91/+504) 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,1570TG 5.8% 631.5 / 3.2 280intergenic (‑91/+504)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (63/188);  new base G (11/15);  total (76/204)
Fisher's exact test for biased strand distribution p-value = 6.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGA  >  CP009273/360009‑360305
                                                                                                                                                    |                                                                                                                                                    
  aCTTCACCGCTTGCCAGCGGCTTACCATCCCGCGCCACCCTCCAGTGCTGGAGCTCCGATGATATGCTGTAAAACATGAGAATTCCGATCCCCAATGCCTGGCTTTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                                   <  1:1465415/150‑1 (MQ=37)
     tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                                <  2:4051278‑M1/60‑1 (MQ=255)
       aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                              <  2:2025791/150‑1 (MQ=255)
       aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                              >  2:2191276/1‑150 (MQ=255)
         cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACggt                                                                                                                                            >  2:4357334/1‑148 (MQ=255)
          agagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                           <  2:963525‑M1/65‑1 (MQ=255)
            agcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCaa                                                                                                                                         <  2:3755133‑M1/67‑1 (MQ=255)
                ttcaccgacaaacaacagataaaacgaacggcccagtctttcgaatgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                     >  2:1988290‑M1/80‑150 (MQ=255)
                cAGCGGCTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                     <  1:1419595/150‑1 (MQ=255)
                cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                     <  2:698014/150‑1 (MQ=255)
                cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                     <  2:1329009/150‑1 (MQ=255)
                 aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGATGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                    <  1:4502262/150‑1 (MQ=255)
                   accgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCaaatcagca                                                                                                                                  <  2:1161066‑M1/74‑5 (MQ=255)
                    ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                 <  2:3158462/150‑1 (MQ=255)
                     gCTTACCATCCCGCGCCACCATCCAGTGCTGGAGCTCCGATGATAAGCTGTCAACCATGAGATTTCCGATCCCAAATTCCTGGAAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                                <  1:1402791/150‑1 (MQ=25)
                     gCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                                <  2:2136250/150‑1 (MQ=255)
                      gcccaacaaccgctaacacgcccggcccagtcggtcgcctgagccgttcgttatattggatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                               >  2:4386164‑M1/74‑150 (MQ=255)
                      gacaaacaacagataaaccgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                               >  2:3264377‑M1/74‑150 (MQ=255)
                      gacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTccc                                                                                                                               <  2:2435006‑M1/77‑1 (MQ=255)
                        caaacaccagagaaaccgcccggcccagtctttcgactgagcctttcgttttatttgctgccgggggaattAATGCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                             >  2:2685011‑M1/72‑150 (MQ=255)
                          ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTAAAATACGCTCCccgcg                                                                                                                           >  1:2366591/1‑150 (MQ=37)
                              ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                       <  2:1850670/150‑1 (MQ=255)
                              acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCggagt                                                                                                                       <  2:2762956‑M1/85‑3 (MQ=255)
                                cgcgcCACCATCCAGTGCTGGCGCTCCGTTGATATGCTGACAAACATGACATTCCCGATCCCCATTTCCTGGCTGTTTTGGGGGGGCGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                     <  1:1996195/149‑1 (MQ=18)
                                  cgcCACCCTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                   <  2:4089443/150‑1 (MQ=255)
                                  ataaaccgaaaggcccagtctttcgcctgagcctttcgttttatttgatgccgggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                   >  2:1580370‑M1/62‑150 (MQ=255)
                                   taaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                  <  2:3333334‑M1/90‑1 (MQ=255)
                                    ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                 >  1:713687/1‑150 (MQ=255)
                                         gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                            <  2:3645249‑M1/96‑1 (MQ=255)
                                          gccAGTGCAGGAGCTCCGATGATAAGCGGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                           <  1:3739751/149‑1 (MQ=255)
                                           ccggggCAGGTGCTCCGATGATACGCTGTCAACCCTGAGAATTCCGAGCCCCACTTCCGGGCTGTTTTGGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                          <  1:1581417/145‑1 (MQ=17)
                                           aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGTCGGTCGTCCTGCCAGCCACCCTCCTGGCCGTTGCTTAGCAACGTTCAAATCAGCTACCGGCGGTATTGTCCTACTCAgg                                                                                                          <  2:89130‑M1/98‑1 (MQ=255)
                                             aGTGCAGGAGCTCCGATTATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                        <  1:650054/150‑1 (MQ=255)
                                              gggCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                       <  2:615790/148‑1 (MQ=255)
                                                gCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                     >  2:4469423/1‑150 (MQ=255)
                                                   gtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAAACCGCTCCCCGCCGATTTGTCCTACTCAGGAGAGCGtt                                                                                                  <  2:1059281‑M1/106‑1 (MQ=255)
                                                   gtaGCTCCGATGATAAGCTGTCAAACTTGAGAATTCCTATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                  <  2:3470288/148‑1 (MQ=255)
                                                    tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCCGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGagcggtc                                                                                                 <  2:3606388‑M1/107‑4 (MQ=255)
                                                       cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                              <  2:2991661/150‑1 (MQ=255)
                                                          gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                           <  2:3299151‑M1/113‑1 (MQ=255)
                                                           gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                          <  1:348067/150‑1 (MQ=255)
                                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                      >  2:2548166‑M1/33‑150 (MQ=255)
                                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                       >  1:4377006‑M1/33‑149 (MQ=255)
                                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCCGGCGAGCGGTCACCGaccaaca                                                                                      <  2:4527760‑M1/118‑1 (MQ=255)
                                                                gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCTGGCGTCCTGCCCGCCACCCTCCTGGCCGTTGCGTGGCAACGTTAAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGGGAGCGTTCTCCGacaaacaa                                                                                     <  1:2904306‑M1/119‑1 (MQ=255)
                                                                   gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                  <  1:282420/150‑1 (MQ=255)
                                                                     tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATCCGCTCCCGGCGGATTTTTCCTACTCCGGAGAGCGGTCACCGACAAACAAAAGAt                                                                                >  1:4155269/1‑150 (MQ=17)
                                                                      gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                               <  2:647818/150‑1 (MQ=255)
                                                                       tCAAACTTTAGAATTCCGATCCCCATTTCCTGGCAGTTTTTGGGGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                              <  1:599571/150‑1 (MQ=14)
                                                                       tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCAACTCAGGCGAGCGGTCACCGACAAACAACAGATaa                                                                              >  2:4257707/1‑150 (MQ=255)
                                                                       gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTAAAAAACCCACCCCGCGGATTTGTACTACTCAAGAGAGCGTTCACCGACAAACAAAAGATaa                                                                              <  1:1026189‑M1/126‑1 (MQ=255)
                                                                        cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                             >  2:2642879/1‑150 (MQ=255)
                                                                         aaaaATGAGAATTCCGATCCCCAATTCCTGGCTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                            <  2:406387/150‑1 (MQ=255)
                                                                         aaaCTTGATAATTCCTATCCCATATTCCTGGATGTTTTGGGGGGGGGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                            <  1:175606/150‑1 (MQ=11)
                                                                         aaaCAGGCGAAGTCCGCGCCCCCATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                            <  2:2671502/150‑1 (MQ=11)
                                                                          aaCATGCGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACGCCGGAGAGCGGTCACCGACCAACAACAGATAAAAc                                                                           >  2:1768846/1‑150 (MQ=14)
                                                                           aCATGATATTTCCGATCCCCAAGACCTGTCTGTTTTTGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                          <  1:2401320/150‑1 (MQ=14)
                                                                            atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                     >  2:3740727‑M1/20‑74 (MQ=255)
                                                                            atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                     <  1:3742836‑M1/55‑1 (MQ=255)
                                                                            atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                           <  2:1080404‑M1/129‑1 (MQ=255)
                                                                              ttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGGGCGTTCACCGACAAACAACAGATAAAAcgata                                                                       <  2:1288483‑M1/133‑3 (MQ=255)
                                                                                    cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAAAGTTCAAAAACGCTCCCCGCGGATTTGTCCTACACAGGTGAGCGTTCCCCGCCCAACAACAGATAAAACGAAAGgcccc                                                                 <  2:3790159‑M1/139‑2 (MQ=255)
                                                                                    cctggggaattAATTCCTGGAAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTAGCAACCTTCAAATACGCTACCGGCGGATTTGTCCTACGAAGGAGTGGGTTAACCGAAAAACAACAGATAAAACGAAAGGCCCa                                                                 <  2:99524‑M1/139‑1 (MQ=255)
                                                                                     tCCGATCCCCATTTCCGGGCTGTTTCTGGCGGGCGACCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                <  1:1329419/150‑1 (MQ=11)
                                                                                      tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTGAAATCCGCTCCCGGCGGATTTGTCCTACGGCGGGGAGGGTTGACCGACCAACAACATATAAAAACAAATGCCCAGt                                                               <  1:3320888‑M1/141‑1 (MQ=255)
                                                                                      tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTAAAATACCCGCCCCGCGGGTTTGTCCAACACAGGAGTGCGTTGCCCGACAAACAACAGATAAAAAGAAaggccccgg                                                               <  1:1009073‑M1/141‑4 (MQ=255)
                                                                                      ccGATCCCAAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                               <  1:3673546/150‑1 (MQ=255)
                                                                                       ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCAGGTGTGCGGTCACCGACCAACAACAGGTAAAACGAATGGCCCAGTc                                                              <  2:3868272‑M1/142‑1 (MQ=255)
                                                                                       ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGGGTTCACCGACAAACAACAGATAAAACGAAAGGCCCCGTc                                                              <  2:3453318‑M1/142‑1 (MQ=255)
                                                                                        gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATTTTTCCTAATCAGGAGAGCGGTAACCGACAAACAACAGATAAAACGAAAGGCCCCGTCt                                                             >  2:4201553/1‑150 (MQ=11)
                                                                                         aTCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                            <  2:552118/150‑1 (MQ=38)
                                                                                          tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                       <  1:1098011/122‑1 (MQ=38)
                                                                                            cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                  >  2:273347/1‑125 (MQ=35)
                                                                                            cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                  <  1:273347/125‑1 (MQ=35)
                                                                                             ccatgtCCTGGCATTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                        <  1:2073402/145‑1 (MQ=255)
                                                                                              tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                       <  2:4037702‑M1/149‑1 (MQ=255)
                                                                                              tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAAACCGCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAAAAGATAAAACGTAAGggccagtcagtcgcc                                                       <  2:3184296‑M1/149‑8 (MQ=255)
                                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCACCTTCAAATCAGCGCCCCGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACCACAGATAAAACGAAAGGCCCAGTCTTTCGAca                                                      >  1:663031/1‑149 (MQ=11)
                                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                     <  2:2353701/150‑1 (MQ=32)
                                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                     <  2:4072185/150‑1 (MQ=32)
                                                                                                 gtCCTGGCTGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                    <  1:3171208/149‑1 (MQ=2)
                                                                                                 gtCCTGGCAGTTTAGGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                    <  1:3217828/149‑1 (MQ=11)
                                                                                                  tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                   >  2:1824917/1‑150 (MQ=11)
                                                                                                  tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                   <  1:2616869/150‑1 (MQ=11)
                                                                                                  tCCTGGCAGTTTATGGCGGGCGGCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                   <  1:972873/150‑1 (MQ=11)
                                                                                                                       gacccgCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGcc                              <  1:1967609/145‑1 (MQ=0)
                                                                                                                         cccgCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTg                            <  1:1665325/147‑1 (MQ=1)
                                                                                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                         <  1:3084418/150‑1 (MQ=31)
                                                                                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                         <  1:3550529/150‑1 (MQ=14)
                                                                                                                            cccccccACCCGCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACGCAGGAGAGCGTTCACCGACAAACACCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                         <  1:4294804/149‑1 (MQ=255)
                                                                                                                             gatctcACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                        <  1:969619/145‑1 (MQ=12)
                                                                                                                              ccGCCACCCTCCGCGCCGCTGCCTCGCAACGGACAACTCCGCTCCCGGCCGATTTGGCCTTCTCCGGAGAGCGTTCACCGACCAACAACAGCTAAAACGAAACCCCCAGACCTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                       <  2:2904306/150‑1 (MQ=11)
                                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCTCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                      <  2:2086399/150‑1 (MQ=18)
                                                                                                                                gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                     <  2:1926054/150‑1 (MQ=34)
                                                                                                                                 ccccccGCCGGGCCGTTGCTTCGCAACGTGCAACTCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGATCACCGACAAACAACAGATAACCCGACAGGCCCAGTCTATCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatta                    <  1:3867051/147‑2 (MQ=11)
                                                                                                                                  cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTCCTCAGGAGAGCGGTCACCGACAAACAACAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                   <  1:51355/150‑1 (MQ=255)
                                                                                                                                           ggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt          <  1:3946506/150‑1 (MQ=255)
                                                                                                                                             ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt        >  2:3675345/1‑150 (MQ=255)
                                                                                                                                                  gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGACCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCGTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAg   <  1:1137981/150‑1 (MQ=255)
                                                                                                                                                   cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGCATTGGGAATTAATTCTTGAAGa  >  1:302692/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGA  >  CP009273/360009‑360305

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: