Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 565,317 G→A 6.6% A320A (GCG→GCA ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273565,3170GA 6.6% 1216.3 / 49.3 348A320A (GCG→GCAybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base G (159/166);  new base A (8/15);  total (167/181)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 353 total aligned reads displayed.

TCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAA  >  CP009273/565169‑565466
                                                                                                                                                    |                                                                                                                                                     
tCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCgg                                                                                                                                                      <  1:2962078/150‑1 (MQ=255)
   tatCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTAt                                                                                                                                                   >  2:2384850/1‑150 (MQ=255)
    atCAGAAGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATc                                                                                                                                                  <  1:373617/150‑1 (MQ=255)
     tCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCt                                                                                                                                                 <  1:1590675/150‑1 (MQ=255)
      cAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTg                                                                                                                                                <  1:2424835/150‑1 (MQ=255)
      cAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTg                                                                                                                                                <  1:4474921/150‑1 (MQ=255)
       aGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCAGTATCTgg                                                                                                                                               >  1:1326561/1‑150 (MQ=255)
           gatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGt                                                                                                                                           <  2:1219321/150‑1 (MQ=255)
             tttGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATCGTTCATGCGGTATCTGGAAGttt                                                                                                                                         >  1:4746318/1‑150 (MQ=255)
                gTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCa                                                                                                                                      >  2:4527745/1‑150 (MQ=255)
                    tACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGc                                                                                                                                  >  2:4861814/1‑150 (MQ=255)
                      cGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCta                                                                                                                                <  1:4445187/150‑1 (MQ=255)
                             cAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGtt                                                                                                                         >  2:2643458/1‑150 (MQ=255)
                              aGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGttt                                                                                                                        >  1:623284/1‑150 (MQ=255)
                              aGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGttt                                                                                                                        >  1:2255708/1‑150 (MQ=255)
                              aGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCAGTATCTGGAAGTTTGCATGGCTATTATGttt                                                                                                                        >  2:5086941/1‑150 (MQ=255)
                               ggAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTg                                                                                                                       >  2:734297/1‑150 (MQ=255)
                               ggAAAGTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTg                                                                                                                       <  2:4955062/150‑1 (MQ=255)
                                gAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGt                                                                                                                      <  1:2609850/150‑1 (MQ=255)
                                 aaaTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTc                                                                                                                     >  2:3230614/1‑150 (MQ=255)
                                  aaTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCAGTATCTGGAAGTTTGCATGGCTATTATGTTTGTcc                                                                                                                    >  1:4709167/1‑150 (MQ=255)
                                   aTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATTTTTGTCCg                                                                                                                   >  2:4350616/1‑150 (MQ=255)
                                     tCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGAt                                                                                                                 >  2:4393172/1‑150 (MQ=255)
                                        gTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATgag                                                                                                              <  1:4624944/150‑1 (MQ=255)
                                        gTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATgag                                                                                                              >  2:1832091/1‑150 (MQ=255)
                                           ggCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATgagagg                                                                                                           <  2:3940897/150‑3 (MQ=255)
                                           ggCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGagaag                                                                                                           >  1:2156778/1‑150 (MQ=255)
                                            gCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGagaaga                                                                                                          >  2:1731273/1‑150 (MQ=255)
                                            gCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGagaaga                                                                                                          <  1:3275203/150‑1 (MQ=255)
                                             cACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGAt                                                                                                         >  1:2775465/1‑150 (MQ=255)
                                             cACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGAt                                                                                                         <  1:4586283/150‑1 (MQ=255)
                                              aCCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGAtt                                                                                                        >  2:4442555/1‑150 (MQ=255)
                                                 ttttGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACa                                                                                                     <  2:4059528/150‑1 (MQ=255)
                                                     ggTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGa                                                                                                 >  1:812226/1‑150 (MQ=255)
                                                         ttAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAtgtg                                                                                             <  1:429969/150‑1 (MQ=255)
                                                            gCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACa                                                                                          >  2:4629309/1‑150 (MQ=255)
                                                             cAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAg                                                                                         >  2:573311/1‑150 (MQ=255)
                                                              aaTAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAgg                                                                                        >  1:3560930/1‑150 (MQ=255)
                                                                tagtagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGcc                                                                                      <  2:1081447/150‑1 (MQ=255)
                                                                 agtagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAAGTGACAGGCCa                                                                                     >  1:4666326/1‑150 (MQ=255)
                                                                  gtagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCAt                                                                                    <  1:4525195/150‑1 (MQ=255)
                                                                  gtagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCAt                                                                                    <  1:2843191/150‑1 (MQ=255)
                                                                   tagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGACGATTACATCGATGTGACAGGCCATc                                                                                   >  2:1551767/1‑150 (MQ=255)
                                                                   tagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATc                                                                                   >  2:4603489/1‑150 (MQ=255)
                                                                   tagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATc                                                                                   >  2:5060752/1‑150 (MQ=255)
                                                                   tagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATc                                                                                   <  1:520362/150‑1 (MQ=255)
                                                                   tagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATc                                                                                   <  1:3226753/150‑1 (MQ=255)
                                                                     gCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCAGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCaa                                                                                 <  2:1869237/150‑1 (MQ=255)
                                                                      cAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAAt                                                                                <  1:1824925/150‑1 (MQ=255)
                                                                       aaGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATa                                                                               >  1:2089676/1‑150 (MQ=255)
                                                                       aaGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATa                                                                               <  1:956311/150‑1 (MQ=255)
                                                                          aaTCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATaaaa                                                                            >  1:3144366/1‑150 (MQ=255)
                                                                           aTCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATaaaaa                                                                           <  2:594704/150‑1 (MQ=255)
                                                                            tCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAag                                                                          <  2:4268687/150‑1 (MQ=255)
                                                                               tATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAgaga                                                                       >  2:83894/1‑150 (MQ=255)
                                                                                aTGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGAt                                                                      >  2:4727471/1‑150 (MQ=255)
                                                                                aTGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCAGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGAt                                                                      >  2:1050677/1‑150 (MQ=255)
                                                                                aTGATAATATATCTTCGTCGAACAGTGATTTAGTGTCAGGCTTGTTGTATTACAATTATGGTTCATGCGGGTTGTTGTAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGAt                                                                      <  2:2698615/150‑1 (MQ=255)
                                                                                  gCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGAttt                                                                    <  1:1368743/150‑1 (MQ=255)
                                                                                  gCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGAttt                                                                    <  1:1770274/150‑1 (MQ=255)
                                                                                  gCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGAttt                                                                    <  1:1475290/150‑1 (MQ=255)
                                                                                   cTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTg                                                                   <  1:1133986/150‑1 (MQ=255)
                                                                                    taataaATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTgg                                                                  <  1:5111051/150‑1 (MQ=255)
                                                                                       taaATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTg                                                               <  2:3509888/150‑1 (MQ=255)
                                                                                        aaaTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGa                                                              >  2:4081286/1‑150 (MQ=255)
                                                                                           tCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATta                                                           <  2:2187937/150‑1 (MQ=255)
                                                                                              aCTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTataa                                                        >  1:204313/1‑150 (MQ=255)
                                                                                              aCTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTataa                                                        <  2:1502838/150‑1 (MQ=255)
                                                                                                ttCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTataata                                                      <  1:4229861/150‑1 (MQ=255)
                                                                                                  cGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTTGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATc                                                    <  2:3395863/150‑1 (MQ=255)
                                                                                                   gAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCTTCAATAAAAAGAGATTTGGTTGATTATAATATCa                                                   <  2:4692525/150‑1 (MQ=255)
                                                                                                     actgTTATGAAGTGTGAGTCTTGTTGTATTACCATTATTGTTCATTCGGTATCTGTAAGTTTGCATGGATATTATGTTTGTCCGTTGAGAAGATTACATCGATGTGACAGGCCTTCTTTATATAGAGATTTGGTTGATTATAATATCAtt                                                 <  2:696240/147‑1 (MQ=255)
                                                                                                         ttATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATg                                             >  2:1406904/1‑150 (MQ=255)
                                                                                                         ttATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATg                                             <  2:2086313/150‑1 (MQ=255)
                                                                                                          tATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGa                                            >  1:4393351/1‑150 (MQ=255)
                                                                                                           aTGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGaa                                           >  2:1165922/1‑150 (MQ=255)
                                                                                                            tGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAAt                                          >  2:2944435/1‑150 (MQ=255)
                                                                                                             gAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAAtt                                         <  1:1932732/150‑1 (MQ=255)
                                                                                                                  gtgAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttt                                    <  2:2437479/150‑1 (MQ=255)
                                                                                                                   tgAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCtttt                                   >  1:147903/1‑150 (MQ=255)
                                                                                                                    gAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttttt                                  >  2:4172332/1‑150 (MQ=255)
                                                                                                                    gAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttttt                                  >  2:3606786/1‑150 (MQ=255)
                                                                                                                    gAGGCTTGTGGTAATACCATGATTGTTCATGCAGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttttt                                  <  2:249768/150‑1 (MQ=255)
                                                                                                                     aGGCTTGTGGTAATACAATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTa                                 <  2:1282035/150‑1 (MQ=255)
                                                                                                                      ggCTTGTGGTAATACCATGATTGTTCATGCAGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTaa                                >  1:4928550/1‑150 (MQ=255)
                                                                                                                      ggCTTGTGGTAATACCATGATTGTTCATGCAGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTaa                                >  1:2530903/1‑150 (MQ=255)
                                                                                                                       gCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAAt                               <  1:1894043/150‑1 (MQ=255)
                                                                                                                        cTTGTGGTAATACCATGATTGTACATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCTATTTGACAGGCCATCAATAAAAAGAGTGTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAAtt                              <  1:2683779/150‑1 (MQ=255)
                                                                                                                           gtgGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTa                           <  1:3880269/150‑1 (MQ=255)
                                                                                                                           gtgGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTa                           >  1:4569885/1‑150 (MQ=255)
                                                                                                                              gTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCa                        <  1:4070543/150‑1 (MQ=255)
                                                                                                                              gTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCa                        <  1:2351120/150‑1 (MQ=255)
                                                                                                                                     cATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCa                 <  1:100535/150‑1 (MQ=255)
                                                                                                                                       tGATTGTTATTTCGGTTTCTGGAAGGTGTCATGGCTTTTATGTTTGTCCGATGAGAAGTTTTCCTCGATGTGCCTGGCCATCATTAAAAAGAGATTTGGTTGATTATACTCTCATTTATGATTTGATTTTTAATTGTAGCTAAATTCAAc               <  1:3356571/150‑1 (MQ=255)
                                                                                                                                         aTTGTTCATGCAGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCa             <  2:2587721/150‑1 (MQ=255)
                                                                                                                                          ttGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAg            <  1:4847810/150‑1 (MQ=255)
                                                                                                                                                  gCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAAc    <  2:3776552/150‑1 (MQ=255)
                                                                                                                                                   cGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACa   <  2:1692476/150‑1 (MQ=255)
                                                                                                                                                    ggTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACaa  >  1:3599424/1‑150 (MQ=255)
                                                                                                                                                    ggTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACaa  <  2:2377651/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAA  >  CP009273/565169‑565466

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