Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,181 G→C 7.0% intergenic (‑115/+480) 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,1810GC 7.0% 1029.6 / 13.3 329intergenic (‑115/+480)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (95/206);  new base C (16/11);  total (111/218)
Fisher's exact test for biased strand distribution p-value = 5.32e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

TACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTA  >  CP009273/360033‑360329
                                                                                                                                                      |                                                                                                                                                    
tACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTCAAATCCGCGCCcg                                                                                                                                                      >  1:3751677/1‑149 (MQ=255)
 aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTAAAATCCGCTCCcgg                                                                                                                                                     <  1:2583058‑M1/79‑1 (MQ=255)
  ccATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                                                    <  2:3601596/149‑1 (MQ=255)
  ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                                                    >  1:3447169/1‑149 (MQ=255)
    aTCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTCAAAACCGCGCCccgcgg                                                                                                                                                  >  1:4308210/1‑149 (MQ=25)
     aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCGCCGGGCCGTTGGGTCGCCACGTGCCAATACGCTCccggcggc                                                                                                                                                 <  2:2210156‑M1/83‑2 (MQ=255)
      ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAAGCCGCGCCCGGCGGAt                                                                                                                                                >  1:225573/1‑149 (MQ=255)
      acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgccgggggaattAATTCCTGGCAGTTGATGGCGGGCGTCCTGCGCGCCACCCTCCGGGCCGGTGCTGCGCAACGTGCACAGGCGCTcccggcgggt                                                                                                                                                <  2:2862181‑M1/84‑3 (MQ=255)
       ccgcgcCACCATCCTGTGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGAGCCCCAATGCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                               <  1:2865555/147‑1 (MQ=21)
       ccccgcCACCAGCCCGTGAAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATACCGACCCCCATTTCCTGGCTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                               <  1:2769780/146‑1 (MQ=21)
       cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                               <  2:4580174/149‑1 (MQ=255)
         gcgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                             <  2:5061131/149‑1 (MQ=255)
         gataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGTTAAAATCCCCACCCCGCGGatttt                                                                                                                                             <  1:2675662‑M1/87‑2 (MQ=255)
          agacaacgaacggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                                            >  2:2590753‑M1/62‑149 (MQ=255)
           taaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAATACGCTCCCGGCGGATTTgta                                                                                                                                           <  2:3132111‑M1/89‑2 (MQ=255)
            ccaccaGCCCGGGCCGGAGCTCCGATTATAAGCTGAAAAACATGAGAAGTCCGAGCCCCACTTCCTGGCAGTTTTGGGGGGGCGTCCTCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                          <  1:1455651/149‑1 (MQ=17)
             aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                                         <  2:2137179‑M1/91‑1 (MQ=255)
               ccaTCCAGTGCAGGAGCTCCGATGATAAGGTGACAAACCTGCCAATTCCGATCCCCAATTCCTGGCATTTTTTGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                                       <  1:4493418/149‑1 (MQ=18)
               acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCAACGTTGAAAAAAGCTCCCGGCGGATTTGTCCTAc                                                                                                                                       <  1:5136209‑M1/93‑1 (MQ=255)
                caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                                      >  2:2724770/1‑149 (MQ=255)
                 gaacggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                                                     >  2:3724026‑M1/55‑149 (MQ=255)
                 gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAAACGCAACCGGCGGAtttgtccaacaa                                                                                                                                     <  1:1585698‑M1/95‑6 (MQ=255)
                 gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGGTCAAATCCGCGCCCGGCGGCTTTGTCCGActg                                                                                                                                     <  2:3990759‑M1/95‑2 (MQ=255)
                  aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCctacaca                                                                                                                                    <  1:839747‑M1/96‑4 (MQ=255)
                      gcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGTTTTGTCCTACTCAGgag                                                                                                                                >  1:1043970‑M1/50‑149 (MQ=255)
                       cccagtctttcgactgagcctttcttattattggatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                               >  1:1114263‑M1/49‑149 (MQ=255)
                       cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                               <  2:2932817‑M1/101‑1 (MQ=255)
                       cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATACGCACCCGGCGGATTTTTGCTACTCAGgaga                                                                                                                               <  1:2533904‑M1/101‑1 (MQ=255)
                        gCAGGAGCTCCTATTATAAGCTGTCAAACATTAGATTTCCGATCCCCAATTCCTGGCTTTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                                              <  1:1976849/149‑1 (MQ=18)
                              tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                                        <  2:280929‑M1/108‑1 (MQ=255)
                              gCTCCGATGATATGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTTGGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTAAGGAGAGCGTTCa                                                                                                                        <  1:4908785/149‑1 (MQ=18)
                              gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                                        <  2:2780923/149‑1 (MQ=255)
                                tcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                                      <  2:1531229‑M1/110‑1 (MQ=255)
                                tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                                      <  1:1308176/149‑1 (MQ=255)
                                acctatgatatgatGTCAAACTTTTCAATTCCTATCCCCAATTCCTGGATGTTTATGGCGGGCTACCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                                                      <  1:948879/145‑1 (MQ=14)
                                 cgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                                                     <  2:1258547‑M1/111‑1 (MQ=255)
                                    ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                                                  <  2:3179548‑M1/114‑1 (MQ=255)
                                    ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCAACCGGCGGATTTGTCCTACTCAGGAGAGCGTTAACCGaca                                                                                                                  <  1:1244391‑M1/114‑1 (MQ=255)
                                    ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCCACGTTCAAATACGCTCCCGGCGGGTTTTTGCTAAACAGGAGGGCGTTCACCGaca                                                                                                                  <  1:4258516‑M1/114‑1 (MQ=255)
                                      gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                                                <  2:1180146/149‑1 (MQ=255)
                                      gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATACGCTCCCGGCGGATTGGTCCGACTCAGGAGAGCGTTCACCGaccaa                                                                                                                >  1:4647458/1‑149 (MQ=18)
                                      gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCACCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGacaaa                                                                                                                <  2:4429714‑M1/116‑1 (MQ=255)
                                        gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCACCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGccaaaca                                                                                                              <  2:231949‑M1/118‑1 (MQ=255)
                                           tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                           <  2:2392814‑M1/121‑1 (MQ=255)
                                           tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCCCTCCCCGCGGATTTGTCCGAATAACGAGAGCGTTCACCGACAAacaaca                                                                                                           <  2:3572979‑M1/121‑1 (MQ=255)
                                           tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCCACGTTCAAAAACCCTCCCCGCCGATTTGTACAAATCAGGATAGCCTTAAACGACAAacaaca                                                                                                           <  1:1640768‑M1/121‑1 (MQ=255)
                                             tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                         <  2:720114/149‑1 (MQ=255)
                                                 aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATACGCACCCGGCGGATTTTTCCTAAACAGGTGAGGCTTGACCGACAaacaaccccaaaa                                                                                                     >  1:2872334/1‑142 (MQ=17)
                                                 aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACAAACAACAGAaaaa                                                                                                     >  1:247658/1‑149 (MQ=255)
                                                  ttatttgatgcctgggtaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                                    >  1:3233249‑M1/22‑149 (MQ=255)
                                                  ttatttgatgcctggggacttAATTCCTGTCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCAACTAAGGCGAGCGTTCACCGACAAACAAAAGATaaaa                                                                                                    <  2:1472563‑M1/128‑1 (MQ=255)
                                                  ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                                    >  1:2597153‑M1/22‑149 (MQ=255)
                                                  aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAAACGTCAAAAACCCGACCGGCGGAATTGTCCAACACAGGAGAGCCGTCACCCACAAACAAAACATaaaa                                                                                                    >  1:976421/1‑149 (MQ=17)
                                                   tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                                   <  2:208735‑M1/129‑1 (MQ=255)
                                                      tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAAACGTGAAATACGCACCCGGCGGGGTGTGCCGACGCCGGAGAGCGGGGCCCGACAAACACCACATAAAAccta                                                                                                >  1:2342080/1‑146 (MQ=9)
                                                      tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGGGCGTTCACCGACAAACCACAGATAAAACGaa                                                                                                >  2:4995999/1‑149 (MQ=255)
                                                      tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAAAACGCACCCCGGCGAATTTTCCCAAAAAGGTTAGCGTTCCACGACCAACAACAAATAAAAAGaa                                                                                                >  1:5040027/1‑149 (MQ=9)
                                                      tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGCGCGTTCACCGACAAACAACCGATAAATCGac                                                                                                >  2:3857650/1‑148 (MQ=255)
                                                       gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCCAATCCGCGCCCGGCGGAGTTGTCCGACGCAGGAGAGCGGTCACCGACACACCACAGAGAAAACGaaa                                                                                               >  2:1924048/1‑149 (MQ=14)
                                                          tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTCAAAACCGCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGGTCACCCACCAACAAAAGATAAAAACAAAGGc                                                                                            <  1:1917218‑M1/136‑1 (MQ=255)
                                                          tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAAAACGCTACCCGCGGATTTGTACAACACAGGAGAGCGGTCACCGACAAACAACACATAAAaagaaatgc                                                                                            <  2:1344470‑M1/136‑4 (MQ=255)
                                                          aaTTCCGATCCCCAATTCCTGGCTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                                            <  1:4190157/149‑1 (MQ=255)
                                                           aTTCCGATCCCATTTTCTTGGCTGTTTTTGGCGGGCTTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGAAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                           <  1:4064639/149‑1 (MQ=11)
                                                            ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                          <  2:871378/149‑1 (MQ=255)
                                                             ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACACACAACAGATAAAACGAAAGGCCCa                                                                                         <  2:2541358‑M1/139‑1 (MQ=255)
                                                               cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGCCCAACGTCAAATCCGCCCCCCGCCGATTTTTCCAAATCAGGAGTGCGGTACCCGACAAACAAAAGATAAAAAGAAAAGcaccct                                                                                       >  1:4124241/1‑144 (MQ=9)
                                                               cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCACCCCGCGGCTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAAAAGATAAAACGAAAGGCCCAga                                                                                       >  1:3857485/1‑148 (MQ=14)
                                                                gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                                      <  2:4806325‑M1/142‑1 (MQ=255)
                                                                gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                                      <  2:4537264‑M1/142‑1 (MQ=255)
                                                                 gtaattATTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                                     >  1:4995999‑M1/7‑149 (MQ=255)
                                                                  gaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                                    >  1:3056433‑M1/6‑149 (MQ=255)
                                                                   aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGGGTTGTCCTACTCAGGAGGGCGGTCACCGACAAACAACAGATAAAAAGAAAGGCCCCGTCttt                                                                                   <  2:492666‑M1/145‑1 (MQ=255)
                                                                    cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                                  <  2:3342923/149‑1 (MQ=35)
                                                                    attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCAACAAACAACAGATAAAACGa                                                                                                 <  2:2966302‑M1/131‑1 (MQ=255)
                                                                    attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAaaggcccactcgtta                                                                                  <  2:2303102‑M1/146‑8 (MQ=255)
                                                                     ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                          <  2:5202985‑M1/106‑1 (MQ=255)
                                                                     ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                          >  1:5203119‑M1/3‑108 (MQ=255)
                                                                      cacttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                                <  1:3700865/146‑1 (MQ=14)
                                                                       aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                           >  2:3909470/1‑105 (MQ=33)
                                                                       aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                                           <  1:3909470/105‑1 (MQ=33)
                                                                       aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                               <  2:4455340/149‑1 (MQ=32)
                                                                        aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                              >  2:4744193/1‑149 (MQ=17)
                                                                        aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                              <  2:3420501/149‑1 (MQ=32)
                                                                         ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                             <  1:3470223/149‑1 (MQ=31)
                                                                          tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                            <  2:5097607/149‑1 (MQ=11)
                                                                                        ggggggcgGCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAt                                                              <  1:468216/146‑1 (MQ=0)
                                                                                                     cccGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCGTTCGACTGAGCCGTTCGGGTTATTTGATggcgggca                                                 >  2:565593/1‑142 (MQ=11)
                                                                                                          cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                                            <  1:3325508/149‑1 (MQ=35)
                                                                                                           aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                                           >  2:788004/1‑149 (MQ=35)
                                                                                                                    gCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt                                  <  2:3525488/149‑1 (MQ=255)
                                                                                                                      cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAGCGACAAACAACAGATAAAACGAAAGGCCCAGTCGTGCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGAAGTGGGacttattgct                                >  2:1963661/1‑144 (MQ=11)
                                                                                                                      cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt                                <  1:4756391/149‑1 (MQ=255)
                                                                                                                          gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCGTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa                            >  2:2706713/1‑149 (MQ=255)
                                                                                                                            ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa                          <  1:3008585/149‑1 (MQ=255)
                                                                                                                             tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACCGATAAAACGACAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAc                         <  1:1525242/149‑1 (MQ=255)
                                                                                                                               cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACg                        <  2:1640768/149‑1 (MQ=255)
                                                                                                                                           aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                                 <  2:2817586/112‑1 (MQ=31)
                                                                                                                                            tCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg           <  1:2098084/149‑1 (MQ=255)
                                                                                                                                               gCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATa        <  1:2137179/149‑1 (MQ=255)
                                                                                                                                               gCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTGCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGACGACGACAGGGCCACGTGATa        >  2:4321424/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGATACGCCTA  >  CP009273/360033‑360329

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: