Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,643,394 A→G 5.7% intergenic (+96/‑245) ydfC → / → dicF uncharacterized protein, Qin prophage/sRNA antisense regulator, inhibits ftsZ, Hfq‑dependent; Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,643,3940AG 5.7% 1138.5 / 40.9 331intergenic (+96/‑245)ydfC/dicFuncharacterized protein, Qin prophage/sRNA antisense regulator, inhibits ftsZ, Hfq‑dependent; Qin prophage
Reads supporting (aligned to +/- strand):  ref base A (176/134);  new base G (9/10);  total (186/145)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 331 total aligned reads displayed.

GGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAG  >  CP009273/1643247‑1643541
                                                                                                                                                   |                                                                                                                                                   
ggAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAGGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGt                                                                                                                                                   >  2:2975505/1‑150 (MQ=255)
 gAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGtt                                                                                                                                                  <  1:2061396/150‑1 (MQ=255)
     aaGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACt                                                                                                                                              >  2:868698/1‑150 (MQ=255)
     aaGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACt                                                                                                                                              <  2:552183/150‑1 (MQ=255)
         agcGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGaaa                                                                                                                                          >  1:1203373/1‑150 (MQ=255)
         agcGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGaaa                                                                                                                                          >  1:2693864/1‑150 (MQ=255)
           cGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAgg                                                                                                                                        <  2:667611/150‑1 (MQ=255)
              cATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGata                                                                                                                                     >  1:4565586/1‑150 (MQ=255)
              cATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGata                                                                                                                                     <  1:2177403/150‑1 (MQ=255)
               aTTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGatat                                                                                                                                    <  1:4945234/150‑1 (MQ=255)
                 tAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATAttt                                                                                                                                  >  1:5261318/1‑150 (MQ=255)
                 tAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATAttt                                                                                                                                  <  1:5338384/150‑1 (MQ=255)
                         gCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGaaa                                                                                                                          >  2:4517586/1‑150 (MQ=255)
                          ccAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGaaaa                                                                                                                         >  1:3474133/1‑150 (MQ=255)
                           cAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAg                                                                                                                        >  2:2542289/1‑150 (MQ=255)
                             gAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTc                                                                                                                      >  1:4723496/1‑150 (MQ=255)
                             gAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTc                                                                                                                      <  2:232434/150‑1 (MQ=255)
                              aaTGATTACATACCAGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCa                                                                                                                     >  1:3695601/1‑150 (MQ=255)
                                 gATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGAc                                                                                                                  <  1:3606154/150‑1 (MQ=255)
                                  aTTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGAcc                                                                                                                 <  2:3985158/150‑1 (MQ=255)
                                     aCATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATa                                                                                                              >  1:2218430/1‑150 (MQ=255)
                                       aTACCCGCTTTTGTTGCTAATCAACAAATCTATGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACTCGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAc                                                                                                            <  2:2991239/150‑1 (MQ=255)
                                       aTACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAc                                                                                                            >  1:75389/1‑150 (MQ=255)
                                        tACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAcg                                                                                                           <  2:1623830/150‑1 (MQ=255)
                                        tACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAcg                                                                                                           >  2:3811862/1‑150 (MQ=255)
                                        tACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAcg                                                                                                           <  2:5636689/150‑1 (MQ=255)
                                        tACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAcg                                                                                                           >  2:882615/1‑150 (MQ=255)
                                          cccGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAcgcg                                                                                                         >  1:140033/1‑150 (MQ=255)
                                            cgcgTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCAAAACGCGaa                                                                                                       >  2:5146447/1‑150 (MQ=255)
                                              cgTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAg                                                                                                     >  1:545710/1‑150 (MQ=255)
                                              cgTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAg                                                                                                     >  2:1015068/1‑150 (MQ=255)
                                              cgTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAg                                                                                                     >  2:3859127/1‑150 (MQ=255)
                                               gTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAgc                                                                                                    >  1:938563/1‑150 (MQ=255)
                                                 tagtagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAgcgc                                                                                                  <  1:3018423/150‑1 (MQ=255)
                                                 tagtagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAgcgc                                                                                                  <  1:1904056/150‑1 (MQ=255)
                                                  agtagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCa                                                                                                 >  1:1054552/1‑150 (MQ=255)
                                                  agtagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCa                                                                                                 >  2:4417969/1‑150 (MQ=255)
                                                   gtagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCAt                                                                                                >  2:1072388/1‑150 (MQ=255)
                                                      gCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGc                                                                                             <  1:2795398/150‑1 (MQ=255)
                                                      gCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGc                                                                                             >  1:1432036/1‑150 (MQ=255)
                                                       cTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCg                                                                                            <  2:2417098/150‑1 (MQ=255)
                                                         aaTCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAg                                                                                          <  2:3172584/150‑1 (MQ=255)
                                                         aaTCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTACACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAg                                                                                          >  2:1332421/1‑150 (MQ=255)
                                                          aTCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAgg                                                                                         <  1:3925885/150‑1 (MQ=255)
                                                          aTCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAgg                                                                                         <  1:1083054/150‑1 (MQ=255)
                                                            caacaaAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTa                                                                                       >  2:1670632/1‑150 (MQ=255)
                                                             aacaaAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAg                                                                                      >  1:5053338/1‑150 (MQ=255)
                                                                aaaGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTg                                                                                   >  2:2769955/1‑150 (MQ=255)
                                                                aaaGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTg                                                                                   >  2:1795933/1‑150 (MQ=255)
                                                                aaaGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTg                                                                                   >  2:4626804/1‑150 (MQ=255)
                                                                   gCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGtt                                                                                <  2:2738052/150‑1 (MQ=255)
                                                                   gCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGtt                                                                                <  1:1625189/150‑1 (MQ=255)
                                                                   gCTAAGGTTAGGAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGtt                                                                                <  2:4945597/150‑1 (MQ=255)
                                                                      aaGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCat                                                                             <  1:1772664/150‑1 (MQ=255)
                                                                        gtttAGAAATTAAGGAGTTCTACACGGGTGAGGTGTATTGCGTGAGCCGTACCCGGGTGCGCATCCGCAACTCACAGATTACTGAAAGGATATATCCCTGAACAGTCAGACCATAACGCGAAATCGCATGGCGAGGTAGCTGGTTCatag                                                                           >  2:2051849/3‑150 (MQ=255)
                                                                        ggatAGTAATTAAGGAGGTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACGGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCatag                                                                           <  2:3972065/147‑1 (MQ=255)
                                                                        ggTTAGTAATTAAGGAGGTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCatag                                                                           <  1:4273725/150‑1 (MQ=255)
                                                                           tAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCatagata                                                                        >  1:1654562/1‑150 (MQ=255)
                                                                            aGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCatagatag                                                                       >  1:1478885/1‑150 (MQ=255)
                                                                             gTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGc                                                                      >  1:2472830/1‑150 (MQ=255)
                                                                              tAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGcc                                                                     >  2:4130185/1‑150 (MQ=255)
                                                                               aaTTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCt                                                                    >  2:813868/1‑150 (MQ=255)
                                                                                  tAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTc                                                                 >  1:2565845/1‑150 (MQ=255)
                                                                                  tAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTc                                                                 <  1:3390743/150‑1 (MQ=255)
                                                                                   aaGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCg                                                                >  2:4175105/1‑150 (MQ=255)
                                                                                    aGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGt                                                               <  2:2173736/150‑1 (MQ=255)
                                                                                          tctcCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAAtt                                                         >  1:3284032/1‑150 (MQ=255)
                                                                                            tcCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCCTCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAAtttt                                                       >  2:3072119/1‑150 (MQ=255)
                                                                                            tcCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAAtttt                                                       <  2:5546959/150‑1 (MQ=255)
                                                                                            tcCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAAtttt                                                       <  1:4189980/150‑1 (MQ=255)
                                                                                            tcCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAAtttt                                                       <  2:741121/150‑1 (MQ=255)
                                                                                             ccACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtc                                                      >  1:5067568/1‑150 (MQ=255)
                                                                                              cACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtcg                                                     >  1:4845462/1‑150 (MQ=255)
                                                                                                cGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtcgtc                                                   <  2:5402570/150‑1 (MQ=255)
                                                                                                 gggTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtcgtcg                                                  <  1:1247139/150‑1 (MQ=255)
                                                                                                 gggTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtcgtcg                                                  <  1:4086820/150‑1 (MQ=255)
                                                                                                  ggTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTCGCTGGTTCATAGATAGCCTGTCGTTAAAGTGTCGTCGa                                                 >  2:2879751/1‑150 (MQ=255)
                                                                                                  ggTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGa                                                 >  2:3476/1‑150 (MQ=255)
                                                                                                    tGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGAcc                                               <  2:4164744/150‑1 (MQ=255)
                                                                                                         tGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTgcg                                          <  1:2813179/150‑1 (MQ=255)
                                                                                                           gAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCt                                        <  2:1907140/150‑1 (MQ=255)
                                                                                                            aGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCtt                                       <  2:3627818/150‑1 (MQ=255)
                                                                                                            aGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCtt                                       <  1:3171981/150‑1 (MQ=255)
                                                                                                               gcgtgcgCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCg                                    <  1:5634906/150‑1 (MQ=255)
                                                                                                                cgtgcgCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCgg                                   >  1:3493748/1‑150 (MQ=255)
                                                                                                                cgtgcgCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCgg                                   >  1:4412466/1‑150 (MQ=255)
                                                                                                                       cGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATCGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGc                            >  2:1457173/1‑150 (MQ=255)
                                                                                                                       cGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGc                            <  2:1406807/150‑1 (MQ=255)
                                                                                                                               ggTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCcgcg                    >  2:3210282/1‑150 (MQ=255)
                                                                                                                                gTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGa                   >  2:3278985/1‑150 (MQ=255)
                                                                                                                                     cATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATgg              >  1:3721453/1‑150 (MQ=255)
                                                                                                                                      aTCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGgc             >  1:4014656/1‑150 (MQ=255)
                                                                                                                                          ggAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCgcgg         <  1:2941021/150‑1 (MQ=255)
                                                                                                                                           gAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCgcggc        <  2:4892192/150‑1 (MQ=255)
                                                                                                                                              cTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGt     >  1:599081/1‑150 (MQ=255)
                                                                                                                                               tGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTa    >  2:1753708/1‑150 (MQ=255)
                                                                                                                                               tGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTa    <  2:12468/150‑1 (MQ=255)
                                                                                                                                                gACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTaa   >  2:4290581/1‑150 (MQ=255)
                                                                                                                                                 aCAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAg  <  1:3902811/150‑1 (MQ=255)
                                                                                                                                                 aCAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAg  >  1:2312975/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAG  >  CP009273/1643247‑1643541

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: