Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 280,968 G→T 60.6% intergenic (+89/‑138) yagF → / → yagG CP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273280,9680GT60.6% 75.7 / inf 286intergenic (+89/‑138)yagF/yagGCP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base G (60/52);  new base T (109/65);  total (169/117)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

Only 100 of 295 total aligned reads displayed.

ATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTA  >  CP009273/280821‑281117
                                                                                                                                                   |                                                                                                                                                     
aTGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGat                                                                                                                                                     >  1:93036/1‑150 (MQ=255)
aTGACACCGATAAAATAATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAAGTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGatat                                                                                                                                                     <  1:3842326/150‑1 (MQ=255)
  gACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatat                                                                                                                                                   >  2:1196665/1‑150 (MQ=255)
    cacCGATATAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCTGAATTTAATTATTGTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatatat                                                                                                                                                 <  2:3355842/150‑1 (MQ=255)
      ccGATAAAATTATCGAGGTAATTAACGCCGGTCAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCg                                                                                                                                               >  1:1100707/1‑150 (MQ=255)
      ccGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCg                                                                                                                                               <  1:2830429/150‑1 (MQ=255)
       cGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGt                                                                                                                                              >  1:3938131/1‑150 (MQ=255)
        taTAAAATTATCGTGTTATTTAACGCCGTTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGtt                                                                                                                                             <  1:768354/149‑1 (MQ=255)
        gATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAATCAATTAGTTGATATATATCGtt                                                                                                                                             >  2:3949589/1‑150 (MQ=255)
        gATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGtt                                                                                                                                             >  2:3391125/1‑150 (MQ=255)
         aTAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTg                                                                                                                                            >  2:2151530/1‑150 (MQ=255)
           aaaaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                                                                                                                                          <  2:5300346/150‑1 (MQ=255)
            aaaTTATCGAGGTAATTAACGCCGGTCAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                                                                                                                                         >  1:3017105/1‑150 (MQ=255)
            aaaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                                                                                                                                         >  2:1340502/1‑150 (MQ=255)
            aaaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                                                                                                                                         <  1:4893072/150‑1 (MQ=255)
                tATCGAGGTAATTAACGCCGGTAAAACAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTcac                                                                                                                                     >  1:2412592/1‑150 (MQ=255)
                tATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTcac                                                                                                                                     >  1:2610604/1‑150 (MQ=255)
                tATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTcac                                                                                                                                     <  2:1070080/150‑1 (MQ=255)
                 aTCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTcaca                                                                                                                                    <  1:4181959/150‑1 (MQ=255)
                      ggTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGc                                                                                                                               >  2:3993430/1‑150 (MQ=255)
                         aaTTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCaaa                                                                                                                            >  1:516071/1‑150 (MQ=255)
                          aTTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCaaaa                                                                                                                           <  2:408367/150‑1 (MQ=255)
                          aTTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCaaaa                                                                                                                           >  1:4657659/1‑150 (MQ=255)
                           ttAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAAt                                                                                                                          <  1:4740886/150‑1 (MQ=255)
                           ttAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAAt                                                                                                                          <  1:3686260/150‑1 (MQ=255)
                           ttAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAAt                                                                                                                          <  2:4785510/150‑1 (MQ=255)
                           ttAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAAt                                                                                                                          <  2:2106811/150‑1 (MQ=255)
                            tAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATa                                                                                                                         >  2:72753/1‑150 (MQ=255)
                             aaCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATaa                                                                                                                        >  2:3181104/1‑150 (MQ=255)
                              aCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAACAGCAAAATAAc                                                                                                                       <  1:2733362/150‑1 (MQ=255)
                               cGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACg                                                                                                                      <  2:1112926/150‑1 (MQ=255)
                                 ccGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTa                                                                                                                    <  2:5011211/150‑1 (MQ=255)
                                 ccGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATAGTTGGCGTCACAAAAGCAAAATAACGTa                                                                                                                    >  1:4335830/1‑150 (MQ=255)
                                  cGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTaa                                                                                                                   >  2:2788758/1‑150 (MQ=255)
                                   ggaaaaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAAt                                                                                                                  <  2:3858364/147‑1 (MQ=255)
                                   ggTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAAt                                                                                                                  >  2:4626505/1‑150 (MQ=255)
                                   ggTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAAt                                                                                                                  <  2:4836589/150‑1 (MQ=255)
                                    gTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAAtt                                                                                                                 <  2:2477016/150‑1 (MQ=255)
                                     tAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTc                                                                                                                <  1:2221273/150‑1 (MQ=255)
                                      aaaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCg                                                                                                               >  2:3612188/1‑150 (MQ=255)
                                          aaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGaaa                                                                                                           <  2:1500/150‑1 (MQ=255)
                                          aaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGaaa                                                                                                           >  1:1955086/1‑150 (MQ=255)
                                                tCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGat                                                                                                     <  1:4045940/150‑1 (MQ=255)
                                                tCGGAATTTAACTATTTTAAGAGATCAAACCGTCTGCGGAATATTTCCCGCAGACGGCATTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGat                                                                                                     <  1:5219577/150‑1 (MQ=255)
                                                 cGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAtata                                                                                                    >  1:692508/1‑150 (MQ=255)
                                                 cGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGata                                                                                                    >  2:5233648/1‑150 (MQ=255)
                                                  ggAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGatat                                                                                                   >  1:4605679/1‑150 (MQ=255)
                                                  ggAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGatat                                                                                                   <  2:1796993/150‑1 (MQ=255)
                                                   gAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATg                                                                                                  <  2:4289004/150‑1 (MQ=255)
                                                   gAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATg                                                                                                  >  1:2034381/1‑150 (MQ=255)
                                                         aaTTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCTCAGAAGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCAtt                                                                                            <  2:5326897/150‑1 (MQ=255)
                                                         aaTTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGTTATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACAAtt                                                                                            >  1:2946375/1‑150 (MQ=255)
                                                          attcttTTCAGCGATCAAACCGTCAGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATAGCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTg                                                                                           <  1:4507844/150‑1 (MQ=255)
                                                              ttttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTgg                                                                                       >  2:497895/1‑150 (MQ=255)
                                                                   agagATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAAt                                                                                  >  1:212439/1‑150 (MQ=255)
                                                                   agagATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAAt                                                                                  >  1:18884/1‑150 (MQ=255)
                                                                   agagATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAAt                                                                                  >  1:5065325/1‑150 (MQ=255)
                                                                    gagaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAAtt                                                                                 <  1:3264128/150‑1 (MQ=255)
                                                                    gagaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCATTGTTGTTATTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAAtt                                                                                 <  1:2086973/150‑1 (MQ=255)
                                                                     agaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTg                                                                                >  1:177334/1‑150 (MQ=255)
                                                                        tAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAAt                                                                             <  2:4482963/150‑1 (MQ=255)
                                                                            aCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCt                                                                         >  1:1437330/1‑150 (MQ=255)
                                                                              cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAATCAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCa                                                                       >  1:1256084/1‑150 (MQ=255)
                                                                              cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCa                                                                       <  1:3610071/150‑1 (MQ=255)
                                                                               gTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCAt                                                                      <  1:4144619/150‑1 (MQ=255)
                                                                                tCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCAtt                                                                     <  2:2584675/150‑1 (MQ=255)
                                                                                  tGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTac                                                                   >  1:558592/1‑150 (MQ=255)
                                                                                   gCGGACTATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTaca                                                                  <  1:5319660/150‑1 (MQ=255)
                                                                                     ggAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACt                                                                >  1:3590092/1‑150 (MQ=255)
                                                                                                cGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGATGAGATATATCGTTGGCGTCACAAAAGCAAAATAAAGTAATTCGGAAATAAGATATGACCAATGCTGGTTAATTGAAAAGCTCATAACACTCCATAAACACg                                                     >  2:3556583/1‑150 (MQ=255)
                                                                                                        gCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGAACATTGCTGGTTAATTGAATATCTCATTACACTCCATTAACACGATGTTGTa                                             >  2:2303387/1‑150 (MQ=255)
                                                                                                           ttgttgTTTTTGAAATTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAAtt                                          >  1:1028428/1‑150 (MQ=255)
                                                                                                            tgttgtTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTc                                         >  2:4673604/1‑150 (MQ=255)
                                                                                                                        aaTTTATTAATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACAt                             >  2:1701120/1‑150 (MQ=255)
                                                                                                                        aaTTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACAt                             >  2:1442803/1‑150 (MQ=255)
                                                                                                                            tattaATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATTACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAgg                         >  2:2556921/1‑150 (MQ=255)
                                                                                                                             attaATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAggg                        <  2:357726/150‑1 (MQ=255)
                                                                                                                              ttaATTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGgt                       <  2:138601/150‑1 (MQ=255)
                                                                                                                               taATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGACACACTACATAAGGgtg                      >  1:1537875/1‑150 (MQ=255)
                                                                                                                                aaTTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGgtgt                     >  2:3061343/1‑150 (MQ=255)
                                                                                                                                aaTTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGgtgt                     >  2:856983/1‑150 (MQ=255)
                                                                                                                                aaTTTAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGgtgt                     >  1:4968986/1‑150 (MQ=255)
                                                                                                                                  tttAAACAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTaa                   <  2:1682256/150‑1 (MQ=255)
                                                                                                                                     aaaCAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTTATTCGGCACACTACATAAGGGTGTAATTc                <  1:2350837/150‑1 (MQ=255)
                                                                                                                                     aaaCAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGATGTAATTCGGCACACTACATAAGGGTGTAATTc                >  1:3055798/1‑150 (MQ=255)
                                                                                                                                     aaaCAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTc                <  2:5065325/150‑1 (MQ=255)
                                                                                                                                      aaCAATTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCt               >  1:2255734/1‑150 (MQ=255)
                                                                                                                                       aCAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCtg              >  2:3132308/1‑150 (MQ=255)
                                                                                                                                         aaTTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCtgat            >  1:336458/1‑150 (MQ=255)
                                                                                                                                         aaTTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCtgat            >  1:3106957/1‑150 (MQ=255)
                                                                                                                                         aaTTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCtgat            >  2:1714410/1‑150 (MQ=255)
                                                                                                                                         aaTTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCtgat            >  2:673928/1‑150 (MQ=255)
                                                                                                                                         aaTTAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCtgat            >  1:1008900/1‑150 (MQ=255)
                                                                                                                                          aTTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCtgatg           <  1:5308408/150‑1 (MQ=255)
                                                                                                                                            tAGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGAc         <  1:268325/150‑1 (MQ=255)
                                                                                                                                             aGTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACg        >  2:777643/1‑150 (MQ=255)
                                                                                                                                              gTTGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGc       >  2:2879632/1‑150 (MQ=255)
                                                                                                                                                tGATATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCaa     >  2:2690169/1‑150 (MQ=255)
                                                                                                                                                  agaTATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAAtt   >  1:379354/1‑150 (MQ=255)
                                                                                                                                                  agaTATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAAtt   <  1:48296/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATAAGATATGACCATTGCTGGTTAATTGAATAGCTCATTACACTCCATTAACACGATGTTGTAATTCGGCACACTACATAAGGGTGTAATTCTGATGACGCAATTA  >  CP009273/280821‑281117

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