Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,173 C→A 6.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 6.4% 720.9 / 6.7 253intergenic (‑107/+488)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (75/151);  new base A (16/10);  total (91/162)
Fisher's exact test for biased strand distribution p-value = 8.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

CAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAT  >  CP009273/360025‑360322
                                                                                                                                                     |                                                                                                                                                     
ttcaccgacaaacaacagataaaccgacaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                      >  2:3169126‑M1/80‑150 (MQ=255)
cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                      <  2:34781/150‑1 (MQ=255)
   gggCTTACCAACCCGCGCCACCATCCTGTGCTGGAGCTCGTATTTTATGCTTACAAACATGAGAATTCCGACCCCCATTGCGGGGCTTTTTTGGGGGGGCGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                                   <  1:912348/149‑1 (MQ=12)
    ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                                  >  2:238792/1‑150 (MQ=255)
       ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                               <  2:3238763/150‑1 (MQ=255)
       ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                               <  2:4989502/150‑1 (MQ=255)
        caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                                              <  2:3954919‑M1/79‑1 (MQ=255)
          ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                            <  1:154895/150‑1 (MQ=255)
          ccATCCAGCGCCACCATCAAGTGCTGGAGCTCCGATGATAAGCTGTAAAACATGAGAATTCCGATCCCAAATTCCTGGAAGTTTATGGCGGGCTTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                                            <  1:118353/150‑1 (MQ=25)
           cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                                           <  2:2183217/150‑1 (MQ=255)
             aacagataaaacgaaaggcccagtctttcgactgagcctttcgttgtatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGTTCGCAACGGTATAAACCCCTccccgcggct                                                                                                                                         <  1:5111563‑M1/84‑3 (MQ=255)
              acagataaaacgaaaggcccagtctttcgactgagccgttcgtgttatttgatgcctggggacttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTGCAAATCcgcgcccggggggtg                                                                                                                                        <  2:572289‑M1/85‑7 (MQ=255)
                agataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                                      <  2:4799427‑M1/87‑1 (MQ=255)
                agataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTAAAATCCGCTCCCGGCGGATTTg                                                                                                                                      <  2:463078‑M1/87‑1 (MQ=255)
                  ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                    <  2:1202572‑M1/89‑1 (MQ=255)
                     cacccTCCCGGGCCGGAGCTCCGATGATACGCGGTCACCCATGAGAAGTCCGAGCCCCAATTCCTGGCAGTTTAGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                 <  1:2471156/150‑1 (MQ=17)
                     caccaTCCAGTGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                 <  1:4358853/150‑1 (MQ=255)
                     caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACCGTAAAATCCGCGCCCGGCGGAGTTGtgcga                                                                                                                                 >  2:1819403/1‑146 (MQ=25)
                        caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                                              <  2:2034120/150‑1 (MQ=255)
                        caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                                              >  1:3114405/1‑150 (MQ=255)
                          aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGGCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                            <  2:5053077‑M1/97‑1 (MQ=255)
                           ccAGTCCAGGCGCTCCGATGATCAGCTGTAAAACATGCGAATACCGATACCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCgg                                                                                                                           <  2:2541272/150‑1 (MQ=18)
                             aGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGCCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                         <  1:5296397/150‑1 (MQ=255)
                              gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGACTTCCGATCCCCAATTCCTGGCTGTTTCTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                                        <  1:120569/150‑1 (MQ=255)
                               cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                                       <  1:682195‑M1/102‑1 (MQ=255)
                                gCAGGAGCTCCGATGATAATCTGTCAAACATGAGAATTCCTATCCCAAATGCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                      <  2:1403416/150‑1 (MQ=18)
                                 cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                     >  2:3738597/1‑150 (MQ=255)
                                 cAGGAGCGCCGAGGATAAGCTGTCAAACATGAGAATTCCGATCCCCACGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                                     <  2:2739679/150‑1 (MQ=255)
                                   ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                                   <  2:4153517/150‑1 (MQ=255)
                                          cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAAAGGTGAAATCAGCGCCCCGCGGATTTGTGCAACACAGGAGAGCGTTCACCCac                                                                                                            >  1:1992482/1‑150 (MQ=18)
                                            ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                                          <  2:1335318‑M1/115‑1 (MQ=255)
                                               atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTAGCAACGTTAAAATACGCTACCGGCGGATTTGTGCGAAACAGGGGTGGGGTAACCGacaaaac                                                                                                       >  1:1598512/1‑148 (MQ=18)
                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                                       <  2:1142065‑M1/118‑1 (MQ=255)
                                               agcctttcgttttatgtgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATACCCTCCCGGCGGATATGTCCTAATCACGATATCGTTAACCCacaaaca                                                                                                       <  1:864304‑M1/118‑1 (MQ=255)
                                                tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGGTCAAATACGCTCCCGGCGGCTTTGTCCTACTAACGAGAGCGTTGACCGacaaaaaa                                                                                                      >  1:1659660/1‑150 (MQ=18)
                                                 cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTCGCAACGTTAAAATACGCGCCCGGCGGATTTGTCCAACTCAGGAGAGCGTTAACCCACAaacaac                                                                                                     <  2:2403197‑M1/120‑1 (MQ=255)
                                                 atgCTGTCAAACCTGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacacc                                                                                                     <  2:1045089/148‑3 (MQ=255)
                                                  aGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                                    <  2:3396300/150‑1 (MQ=255)
                                                   tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                   <  1:2165190‑M1/122‑1 (MQ=255)
                                                    ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGCGTTAACCGACAAACAACAGa                                                                                                  <  2:1469367‑M1/123‑1 (MQ=255)
                                                      cgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACagaaa                                                                                                <  2:4038159‑M1/125‑3 (MQ=255)
                                                      cgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATa                                                                                                <  2:3009600‑M1/125‑1 (MQ=255)
                                                      cgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAATACGCTCCCCGCGGATTTGTCCTACTCATGTTAGCGTTGACCGACAAACAACAGATa                                                                                                <  1:110904‑M1/125‑1 (MQ=255)
                                                       tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                               <  2:1598715/150‑1 (MQ=255)
                                                        ttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                              >  1:1338207‑M1/24‑150 (MQ=255)
                                                        cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACACACAACAGATaaa                                                                                              >  2:4854329/1‑150 (MQ=255)
                                                         gttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATGCGCTCCCGGCGGAGTTGTCCGACTCAGGTGGGCGGTCACCGACCAACAACAGAtaaat                                                                                             <  1:4808848‑M1/128‑2 (MQ=255)
                                                          aaCATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                            <  1:3644010/150‑1 (MQ=255)
                                                           aaaTGTGAATGCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                           <  2:4496535/148‑1 (MQ=255)
                                                           aCATGAGAATTCCGATCCCCAATTCCTGGAAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                           <  2:1895604/150‑1 (MQ=255)
                                                             aTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTCAAATCCGCGCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACCAACAACACAAAAAACGaa                                                                                         >  1:2927538/1‑150 (MQ=14)
                                                              ttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                        >  2:4390626‑M1/18‑150 (MQ=255)
                                                              tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                                        <  2:1562703/150‑1 (MQ=255)
                                                               tgatgcgtggggaattAATTCCTGGCTGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                       >  1:1486785‑M1/17‑150 (MQ=255)
                                                               gagaATTCCGATCCCCAATTCCTGGCTTTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                       <  2:1992482/150‑1 (MQ=255)
                                                               gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTGCGCACCGGTCAAATACGCACCCGGCGGATTTGTGCGACTCAGGAGAGCGGTCACCGACAAAAAACAGATAAAACGAAAg                                                                                       >  1:4716824/1‑150 (MQ=14)
                                                                gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                   <  2:2687391‑M1/122‑1 (MQ=255)
                                                                agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTGCGCCGGAGTGCGGTGACCGACACACCACAGATAGATCGAAAgg                                                                                      >  1:3835466/1‑150 (MQ=14)
                                                                 atgcctggggaattATTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                                     >  1:3726002‑M1/15‑150 (MQ=255)
                                                                    ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                  <  2:5190703/150‑1 (MQ=255)
                                                                    cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                  >  1:1296677‑M1/12‑150 (MQ=255)
                                                                      tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                                <  1:3075426‑M1/141‑1 (MQ=255)
                                                                       cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGTTGGCAAAGGTACAATACCGAACCCGCGGTTTTGTGGTAATCAGGAGGGGGGTGACCCGCCAACACACGATAAAACGAAAGGGCCAGTc                                                                               >  1:1907454/1‑150 (MQ=9)
                                                                        gggttttATTTCCTGGTTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCGGGGCCGTTGCTTCGCAACGTTCATTTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                              >  1:4799312‑M1/8‑150 (MQ=255)
                                                                        gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTGGCAACGGTAAAATACGCTCCCGGCGGATTTGTCCAACGCAGGGGGGCGGTCACCGACAAACAACAGATAAAACGAAAGGGCCAGTca                                                                              >  1:1779075/1‑149 (MQ=11)
                                                                         aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                             <  2:539248/150‑1 (MQ=38)
                                                                         aTCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                             <  2:8017/150‑1 (MQ=38)
                                                                           ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                           <  1:3912691/150‑1 (MQ=35)
                                                                           ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACGGTAAAATACGCGCCCGGCGGCGTTGTCCAACGCCGGAGAGCGGTCACCGACCAACAACAGATAAAACGAAAGGCCCCGTCTTTc                                                                           >  1:3359928/1‑150 (MQ=11)
                                                                           ccccAATTCCGGGCTGTTTATGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                           <  1:2785839/150‑1 (MQ=255)
                                                                            cccAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                          <  1:4859220/150‑1 (MQ=25)
                                                                            cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                          <  1:5263319/150‑1 (MQ=35)
                                                                            cccAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACACTCCGGGCCGTTGATTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                                          <  2:3884653/150‑1 (MQ=255)
                                                                             ccATTTCCTGGCTGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                                         <  1:3202369/150‑1 (MQ=12)
                                                                               atttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                             <  1:5239714/126‑1 (MQ=14)
                                                                               aaTTCCTGGCAGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                         <  1:1026072/84‑1 (MQ=25)
                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                               <  1:4180340/78‑1 (MQ=35)
                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                                       >  2:5302409/1‑150 (MQ=32)
                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAAGGTCAAAACCGCTCCCGGCGGTTTTGTCCGACACAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAATGGCCAGTCattagtca                                                                       >  1:2624346/1‑145 (MQ=11)
                                                                               aaTTCCTGGCAGTTTATGGAGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                                         >  2:1026072/1‑84 (MQ=25)
                                                                                                             cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                         <  1:4739193/150‑1 (MQ=32)
                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                       <  2:2498453/150‑1 (MQ=18)
                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                       <  2:835760/150‑1 (MQ=33)
                                                                                                                gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                                      <  1:4152344/150‑1 (MQ=34)
                                                                                                                   aCCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCAACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTAGTCTGAGCCTTTCGTTTTATTTGTTGCCTGgggcattggg                                   >  2:4232590/1‑150 (MQ=11)
                                                                                                                      cTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACAAAAAACAGATAAAACGAAAGGCCCAGTCTTTAGACAGAGCCGTTCGTTTTATTTGATGCCTGgggaattaaatca                                >  2:477143/1‑144 (MQ=11)
                                                                                                                         cGGGCCGTTGCTTCGCAAAGGTCAAATCCGCTCCCGGCGGATTTGTCCTAATCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCCGTCTTTCGACTGAGCCTTTCGTTTTATTTGAGGGCTGGGGTATTGGGcattaa                             >  2:3660719/1‑150 (MQ=11)
                                                                                                                             ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt                         <  1:565081/150‑1 (MQ=255)
                                                                                                                             ccGTTGCGTCGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACACACCACAGATAAACCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt                         <  2:5111563/150‑1 (MQ=14)
                                                                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGACTTAATTCTTg                       >  1:2586129/1‑150 (MQ=255)
                                                                                                                                 tGCTTCCCAACGTACAAATCCGCACCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATACAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGaa                     <  1:4016623/150‑1 (MQ=255)
                                                                                                                                          aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAAGTGGGAATTAATTCTTGAAGACGAAAg             >  1:3268377/1‑150 (MQ=255)
                                                                                                                                           aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAgg            <  2:864304/150‑1 (MQ=255)
                                                                                                                                           aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAgg            >  2:1787234/1‑150 (MQ=255)
                                                                                                                                              ttCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGcc         <  2:682195/150‑1 (MQ=255)
                                                                                                                                              ttCAAATCCGCTCCCGGCGGAGTTGTCCTACTCAGGAGAGCGGTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGTCGTGTTCTTTGAGGCCTGGGGACTTGGGACTTAATTCTTGCAGACGACAGGGcc         >  2:4777014/1‑150 (MQ=17)
                                                                                                                                                 aaaTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGACTTGGGAATTAATTCTTGAAGACGAAAGGGCCTCg      >  2:3048320/1‑150 (MQ=255)
                                                                                                                                                   aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg    >  1:275061/1‑150 (MQ=255)
                                                                                                                                                   aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg    >  1:3439281/1‑150 (MQ=255)
                                                                                                                                                   aTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTg    >  1:4081456/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGAT  >  CP009273/360025‑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: