Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,173 C→A 5.4% 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 5.4% 1036.2 / 10.8 337intergenic (‑107/+488)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (105/201);  new base A (15/16);  total (120/217)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 353 total aligned reads displayed.

CCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAT  >  CP009273/360024‑360322
                                                                                                                                                     |                                                                                                                                                     
ccAGCGGCTGCCCCTCCAGCCCCACCCGCCAGTGCAGGAGCACCGAGCTTACGCTGTCCAACACGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGCGCGGCCGACCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                                       <  2:2396772/150‑1 (MQ=18)
 ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                      <  2:5525865‑M1/71‑1 (MQ=255)
 cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                      <  2:103285/150‑1 (MQ=255)
 cAGCGCCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                      <  2:428379/150‑1 (MQ=255)
  aGCGGCTTACCATCCAGCGCCACCATCCTGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                                     <  1:4095249/150‑1 (MQ=255)
     ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGAGAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                                  >  2:5351973/1‑150 (MQ=255)
        acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                               >  2:2996764‑M1/73‑150 (MQ=255)
         tACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGGCGCCACGGTCAAAACCCCTccccg                                                                                                                                              >  1:2516789/1‑148 (MQ=37)
            acaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                                           <  2:4989925‑M1/82‑1 (MQ=255)
              tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcggct                                                                                                                                         >  1:3040660/1‑148 (MQ=255)
              aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                         <  2:2673247‑M1/84‑1 (MQ=255)
                 agataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                      <  2:2030977‑M1/87‑1 (MQ=255)
                 agataaaacgaaaggcccagtctttcgactgagcctttcgttttagttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                      <  2:4982182‑M1/87‑1 (MQ=255)
                  gcgcCCCCCTCCAGTGCAGGAGCTCCGATGATAAGCGGTCAACCATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                                     <  2:3877037/150‑1 (MQ=21)
                     ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                                  <  2:4758490/150‑1 (MQ=255)
                     ccaccaTCCAGTGCAGGAGATCCGATGATATGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                                  <  2:710876/150‑1 (MQ=255)
                     aaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                                  >  2:3981398‑M1/60‑150 (MQ=255)
                      caacaTCCCGTGCAGTAGCTCCGCTGATAAGCTGTCAAAAATGAGAATACCGAGCCCAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                 <  2:1260903/147‑1 (MQ=18)
                       accaTCCAGTGCAGGAGCTCCTATGATAAGCTGACAAACATGATAATTCCCACCCCAAATTCCTGGAAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                                <  1:2378957/150‑1 (MQ=18)
                       accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                                <  2:5451859/150‑1 (MQ=255)
                           aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGGTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                            <  2:1351892‑M1/97‑1 (MQ=255)
                             cggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                                          >  1:2657982‑M1/52‑150 (MQ=255)
                                tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTGTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                                       >  1:4986015/1‑150 (MQ=255)
                                cccagtctttcgactgagcctttcgttttatttgatgcctggggacttAATACCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCGTCGCAACGTTACAATACGCACCCGGCGGAGGTGTCCTACTCAGGagag                                                                                                                       <  2:1058266‑M1/102‑1 (MQ=255)
                                   agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTTGCAGTTTATGGCGGGCGTCCTGCCAGCCACGCTCCGGGCAGTTGCTGCGCAACGTTCAAATCCGCTCCCGGCGGATGTGTCCTACTCAGGAGAGcgg                                                                                                                    <  2:932474‑M1/105‑2 (MQ=255)
                                    gccttacgactgcgcctttcgttttatttgagccctggggaataAATTCCGGGCAGTTTATGGCGGGCGTCCTCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                   >  1:4822036‑M1/45‑150 (MQ=255)
                                     gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                                  <  2:1811720/150‑1 (MQ=255)
                                      aGCTCCGATGATAAGCTGTCAAACATGAGAAGCCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                                 <  1:4344083/150‑1 (MQ=255)
                                       gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                                <  2:4116273/150‑1 (MQ=255)
                                           gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATCCGCTACCGGCGGATTTGTCCTACTCAGGTGAGCGGTCACCGac                                                                                                            <  1:1728400‑M1/113‑1 (MQ=255)
                                           cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                                            <  2:1616082/150‑1 (MQ=255)
                                               gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                        <  2:2563898/150‑1 (MQ=255)
                                               gagcctttcgttttatttgatgcctggggacttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAccgacacac                                                                                                        <  2:3513124‑M1/117‑4 (MQ=255)
                                                agcctttcgtattatttgttgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                                       >  1:4862748‑M1/33‑150 (MQ=255)
                                                   ctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCAACACAGGCGAGCGTTCCCCGACAAacaaca                                                                                                    <  2:5183986‑M1/121‑1 (MQ=255)
                                                   aGCTGTAAAAATTTATAATTCCGAGACCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                    <  1:1699831/150‑1 (MQ=14)
                                                    tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                   <  1:3164816‑M1/122‑1 (MQ=255)
                                                     cTGACAACCATGAGAATTCCGATCCCCAATGCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                  <  1:717694/150‑1 (MQ=255)
                                                       gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                <  2:4668318/150‑1 (MQ=255)
                                                       gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                                <  2:3100794/150‑1 (MQ=255)
                                                       gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATACGCTCCCGGCGGATTTGTCCTACTCATGATAGCGTTCACCGACAAACAAAAGaaa                                                                                                >  1:417768/1‑148 (MQ=17)
                                                        tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                               <  2:2134381/150‑1 (MQ=255)
                                                        tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                               <  2:417850/150‑1 (MQ=255)
                                                           ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                            >  2:1701210‑M1/22‑150 (MQ=255)
                                                            aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGGCCTACTCAGGTGAGCGTTCACCGACAAACAACAGATaaaaag                                                                                           >  2:2474362/1‑148 (MQ=255)
                                                             atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTAGCAACGTTAAAATACGCTCCCGGCGGATTTGTACTACACAGGAGAGCGTTAACCGACAAAAAAaagaaaaatcca                                                                                          <  1:2238187‑M1/131‑5 (MQ=255)
                                                               tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                        <  2:5232579/150‑1 (MQ=255)
                                                                tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAacgatag                                                                                       <  2:5063625‑M1/134‑4 (MQ=255)
                                                                gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCCACGTTCAAATACGCGCCCGGCGGAGTTGTCCTACTCAGGAGGGCGGTCACCGACAAACAACACATAAAACGAAAg                                                                                       >  1:4907285/1‑150 (MQ=14)
                                                                 agaTTTCCGAGCCCCAAGTCCTGGCTGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                      <  1:3805957/150‑1 (MQ=11)
                                                                  atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGCCGCCGGAGAGCGTTCACCGCCCAACACCAGATAAAACGAAAGGc                                                                                     <  2:3349563‑M1/136‑1 (MQ=255)
                                                                       tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTACCCACCGTCAAATCCGCTCCCCGCGGAATTGTCCAACTCACGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGCCCAGt                                                                                <  1:4148274‑M1/141‑1 (MQ=255)
                                                                           tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                                            <  1:2767578/150‑1 (MQ=37)
                                                                            ccgaaagtCCTGGCAGTGGCGGGCGGGCGGCCGGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTCCTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGACCGGCCCAGTCTTTc                                                                           <  1:4680976/143‑1 (MQ=11)
                                                                            ccccAATTCCTGGAAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                           <  1:1140780/150‑1 (MQ=25)
                                                                            aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGCGAGCGTGCACCGACCAACAACAGGTAAAACGCACGGCCCAGtcttgc                                                                           <  2:3402729‑M1/146‑3 (MQ=255)
                                                                             attaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCCACGGTGAAAAACCCTCCCCGCGGGTTTTGGCTACACAGGTGAGCGTTAACCCTCAAAAAACAGATTAAACCGAAGGCCCAGTCTTTCg                                                                          >  1:1069120/4‑150 (MQ=2)
                                                                             attATTTCCTGGCAGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                          >  2:31716‑M1/4‑150 (MQ=255)
                                                                             attAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                          <  2:4508002‑M1/147‑1 (MQ=255)
                                                                              ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                    <  2:224951‑M1/121‑1 (MQ=255)
                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                                      <  1:4998493/135‑1 (MQ=33)
                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                       >  2:1198144/1‑150 (MQ=32)
                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                       <  2:3593679/150‑1 (MQ=32)
                                                                                 aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                                      <  1:4556028/150‑1 (MQ=32)
                                                                                  ttCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCGTACTCAGGAGAGCGTTCAACGACAAAAACCAGATAAAACGATAGTCCCAGTCTTTGGACTGa                                                                     <  2:1984384/150‑1 (MQ=11)
                                                                                   tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCAACGTTCAAAAACGCTACCCGCGGATTTTTACAACTAAGGAGAGCGTTGACCGAAAAAAAACAGATAAAAAGAAAGTCCCCGTCTTTCGaaagag                                                                    >  1:4408439/1‑145 (MQ=11)
                                                                                   tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACGTTAAAATACGCTCCCCGCGGATTTTTGCTACTCAGGTGAGCGTTCACCGACAAACAACAGATAAAAAGCAAGGCCCAGTCGTTCGagagag                                                                    >  1:4524148/1‑145 (MQ=11)
                                                                                   tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTggg                                                                    >  2:4495239/1‑148 (MQ=11)
                                                                                              tagggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGtttt                                                         <  2:3033791/147‑1 (MQ=0)
                                                                                                     gcgctcgtGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGAt                                                  <  1:2578881/143‑1 (MQ=0)
                                                                                                        gcccggCCCCCCACCCTCCGGTCCGTTGCTTCGCACCGTTCAAATCCGCTCCCGGCGGATTTGTCCTCCTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACTAAAGGCCCAGTCTTTCGACTGAGCCTTTCGATTTATTTGATGcc                                               <  1:1001778/145‑1 (MQ=1)
                                                                                                            tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                           <  2:146497/150‑1 (MQ=11)
                                                                                                            tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                           <  1:5525865/150‑1 (MQ=11)
                                                                                                             gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                                          <  1:3175077/150‑1 (MQ=31)
                                                                                                              cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                         <  1:2899181/150‑1 (MQ=32)
                                                                                                                cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGTGCCTTTCGTTTTATTTGTTGGCTGgggaat                                       >  2:2945286/1‑150 (MQ=16)
                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                      <  1:3464901/150‑1 (MQ=34)
                                                                                                                          cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                                    >  2:3632565/1‑143 (MQ=35)
                                                                                                                          cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                                    <  1:3632565/143‑1 (MQ=35)
                                                                                                                           gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat                            >  2:5378180/1‑150 (MQ=255)
                                                                                                                           gggCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat                            <  1:828354/150‑1 (MQ=255)
                                                                                                                             gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                         >  1:369692/1‑135 (MQ=32)
                                                                                                                             gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                         <  2:369692/135‑1 (MQ=32)
                                                                                                                             gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGTATTAATTc                          >  1:5121852/1‑150 (MQ=255)
                                                                                                                               cGTTGCTTCGCACCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                        <  1:3097148/150‑1 (MQ=255)
                                                                                                                               cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt                        >  2:5000551/1‑150 (MQ=255)
                                                                                                                                gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTg                       <  1:5346876/150‑1 (MQ=255)
                                                                                                                                 ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTGTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa                      >  2:2279214/1‑150 (MQ=255)
                                                                                                                                  tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACTGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa                     <  1:394700/150‑1 (MQ=255)
                                                                                                                                  tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGCATTGGGAATTAATTCTTGaa                     >  1:2113377/1‑150 (MQ=255)
                                                                                                                                    cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAAAACTTGatgt                   >  2:1300711/1‑147 (MQ=255)
                                                                                                                                         cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGaaa              <  1:1285700/150‑1 (MQ=255)
                                                                                                                                          aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAg             <  1:1609200/150‑1 (MQ=255)
                                                                                                                                             gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGc          <  1:4508002/150‑1 (MQ=255)
                                                                                                                                               tCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAATGCCCAGTCTTTCGAATGAGCCTTTCGTTTTAGTTGATGACTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCt        >  1:4483652/1‑150 (MQ=17)
                                                                                                                                               tCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCt        >  2:2866590/1‑150 (MQ=255)
                                                                                                                                                cAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTc       >  1:3309234/1‑150 (MQ=255)
                                                                                                                                                 aaaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCg      >  1:5507894/1‑150 (MQ=255)
                                                                                                                                                   aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg    >  2:856267/1‑150 (MQ=255)
                                                                                                                                                    tCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGa   >  2:2644978/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAT  >  CP009273/360024‑360322

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: