Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,643,394 A→G 5.4% 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.4% 1099.2 / 36.6 315intergenic (+96/‑245)ydfC/dicFuncharacterized protein, Qin prophage/sRNA antisense regulator, inhibits ftsZ, Hfq‑dependent; Qin prophage
Reads supporting (aligned to +/- strand):  ref base A (178/119);  new base G (10/7);  total (189/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 318 total aligned reads displayed.

TGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAG  >  CP009273/1643246‑1643541
                                                                                                                                                    |                                                                                                                                                   
 ggAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAg                                                                                                                                                    >  2:826669/1‑149 (MQ=255)
  gAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGt                                                                                                                                                   >  1:3252855/1‑149 (MQ=255)
     aaaGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTa                                                                                                                                                <  1:2031212/149‑1 (MQ=255)
       agcagcGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACt                                                                                                                                              >  2:5104168/1‑149 (MQ=255)
        gcagcGGCATTAAGCATTGCGAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTg                                                                                                                                             >  2:24411/1‑149 (MQ=255)
           gcGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGaaa                                                                                                                                          >  1:1204241/1‑149 (MQ=255)
           gcGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGaaa                                                                                                                                          >  1:2721633/1‑149 (MQ=255)
              gCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGa                                                                                                                                       >  2:4482449/1‑149 (MQ=255)
              gCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGa                                                                                                                                       >  2:1889303/1‑149 (MQ=255)
                aTTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGata                                                                                                                                     <  2:1254957/149‑1 (MQ=255)
                aTTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGata                                                                                                                                     <  2:663925/149‑1 (MQ=255)
                 ttAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGatat                                                                                                                                    >  1:1179859/1‑149 (MQ=255)
                 ttAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGatat                                                                                                                                    >  1:1643643/1‑149 (MQ=255)
                 ttAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGatat                                                                                                                                    <  1:2539992/149‑1 (MQ=255)
                   aaGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATAttt                                                                                                                                  <  2:3758988/149‑1 (MQ=255)
                    aGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTc                                                                                                                                 <  1:2813616/149‑1 (MQ=255)
                    aGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTc                                                                                                                                 >  1:1485827/1‑149 (MQ=255)
                      cATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTccc                                                                                                                               >  2:109209/1‑149 (MQ=255)
                            cAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAAGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAGAGGATATTTCCCTGaaaa                                                                                                                         <  1:1425968/149‑1 (MQ=255)
                            cAGAATGATTACAGACCCGCGTTAGGAGCTAATCAACAAAGCTAAGGTTAGTAATTCAGGAGTGCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGaaaa                                                                                                                         <  2:3942536/149‑1 (MQ=255)
                                aTGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCa                                                                                                                     >  1:4712530/1‑149 (MQ=255)
                                    ttACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGAcc                                                                                                                 >  2:3349299/1‑149 (MQ=255)
                                     tACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCa                                                                                                                >  1:387525/1‑149 (MQ=255)
                                       cATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATa                                                                                                              >  2:4040457/1‑149 (MQ=255)
                                         tACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAc                                                                                                            >  1:1566507/1‑149 (MQ=255)
                                           cccGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAcgc                                                                                                          >  1:3717105/1‑149 (MQ=255)
                                            ccGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAAcgcg                                                                                                         <  1:2902594/149‑1 (MQ=255)
                                              gcgTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGaa                                                                                                       >  1:1527140/1‑149 (MQ=255)
                                                gTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAg                                                                                                     <  2:1853742/149‑1 (MQ=255)
                                                    gtagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACCCGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCa                                                                                                 >  2:1508784/1‑149 (MQ=255)
                                                    gtagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCa                                                                                                 <  1:3507518/149‑1 (MQ=255)
                                                     tagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCAt                                                                                                <  1:1014123/149‑1 (MQ=255)
                                                     tagCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCAt                                                                                                >  2:4178396/1‑149 (MQ=255)
                                                       gCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATgg                                                                                              >  2:3094742/1‑149 (MQ=255)
                                                         tAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCg                                                                                            >  1:4263829/1‑149 (MQ=255)
                                                         tAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCg                                                                                            <  1:4394482/149‑1 (MQ=255)
                                                         tAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCg                                                                                            >  1:655203/1‑149 (MQ=255)
                                                         tAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCg                                                                                            >  2:2836051/1‑149 (MQ=255)
                                                         tAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCg                                                                                            >  2:1507043/1‑149 (MQ=255)
                                                           aTCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAg                                                                                          >  2:3491649/1‑149 (MQ=255)
                                                            tCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAgg                                                                                         <  1:4905502/149‑1 (MQ=255)
                                                            tCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAgg                                                                                         >  2:2786968/1‑149 (MQ=255)
                                                             caacaaAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGt                                                                                        <  1:3577082/149‑1 (MQ=255)
                                                             caacaaAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGCGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCCTGGCGAGGt                                                                                        >  2:4516865/1‑149 (MQ=255)
                                                              aacaaAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTa                                                                                       <  1:4466330/149‑1 (MQ=255)
                                                                 aaaGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCt                                                                                    <  1:3772002/149‑1 (MQ=255)
                                                                     cTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGtt                                                                                >  1:1276204/1‑149 (MQ=255)
                                                                       aaGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCa                                                                              >  2:1611073/1‑149 (MQ=255)
                                                                       aaGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCa                                                                              >  1:4775566/1‑149 (MQ=255)
                                                                       aaGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCa                                                                              <  1:3084304/149‑1 (MQ=255)
                                                                        aGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCat                                                                             >  2:5031403/1‑149 (MQ=255)
                                                                                 aTTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCt                                                                    >  1:693455/1‑149 (MQ=255)
                                                                                 aTTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCt                                                                    >  1:5014388/1‑149 (MQ=255)
                                                                                 aTTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCt                                                                    >  2:1406946/1‑149 (MQ=255)
                                                                                  ttAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTg                                                                   >  1:4979389/1‑149 (MQ=255)
                                                                                    aaGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTc                                                                 >  2:227394/1‑149 (MQ=255)
                                                                                    aaGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTc                                                                 <  2:1628556/149‑1 (MQ=255)
                                                                                    aaGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTc                                                                 <  2:2628159/149‑1 (MQ=255)
                                                                                    aaGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTc                                                                 <  2:4866304/149‑1 (MQ=255)
                                                                                     aGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCg                                                                >  1:4625306/1‑149 (MQ=255)
                                                                                     aGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCg                                                                >  1:190190/1‑149 (MQ=255)
                                                                                       gAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGtt                                                              >  2:2375692/1‑149 (MQ=255)
                                                                                        aGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTa                                                             >  1:2795403/1‑149 (MQ=255)
                                                                                         gTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTaa                                                            >  2:1997488/1‑149 (MQ=255)
                                                                                         gTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTaa                                                            >  1:2635726/1‑149 (MQ=255)
                                                                                          ttCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTaaa                                                           >  2:4818943/1‑149 (MQ=255)
                                                                                               cACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtc                                                      <  2:3216590/149‑1 (MQ=255)
                                                                                                 cGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtcgt                                                    <  2:1393571/149‑1 (MQ=255)
                                                                                                 cGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtcgt                                                    <  1:5199329/149‑1 (MQ=255)
                                                                                                  gggTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTtcgtc                                                   >  2:716725/1‑149 (MQ=255)
                                                                                                      gAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGAcc                                               <  2:3057442/149‑1 (MQ=255)
                                                                                                        ggTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGt                                             >  2:2888618/1‑149 (MQ=255)
                                                                                                         gTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTg                                            >  1:3311331/1‑149 (MQ=255)
                                                                                                         gTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTg                                            >  1:1093447/1‑149 (MQ=255)
                                                                                                          tGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTgc                                           <  2:903660/149‑1 (MQ=255)
                                                                                                          tGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTgc                                           <  2:2958790/149‑1 (MQ=255)
                                                                                                           ggAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTgcg                                          >  1:4574429/1‑149 (MQ=255)
                                                                                                                gcgtgcgCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTGAAATTTTCGGCGTCCGTGCGCGTcc                                     >  1:550270/1‑149 (MQ=255)
                                                                                                                 cgtgcgCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCg                                    >  2:2536147/1‑149 (MQ=255)
                                                                                                                  gtgcgCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTAGTCGACCGTGCGCTTCCgg                                   >  1:4860118/1‑149 (MQ=255)
                                                                                                                   tgcgCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGt                                  <  1:4013315/149‑1 (MQ=255)
                                                                                                                   tgcgCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGt                                  <  1:5031403/149‑1 (MQ=255)
                                                                                                                    gcgcCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGtt                                 <  2:2078348/149‑1 (MQ=255)
                                                                                                                      gcCGGACACGGGTGCGCATCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTtgt                               <  2:4346973/149‑1 (MQ=255)
                                                                                                                          gACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCa                           >  2:1202417/1‑149 (MQ=255)
                                                                                                                                ggTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCcgc                     <  2:533178/149‑1 (MQ=255)
                                                                                                                                  tGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTACGGTTGTGGCAACCCGCGa                   >  1:4139949/1‑149 (MQ=255)
                                                                                                                                    cgcATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGaaa                 <  1:3027076/149‑1 (MQ=255)
                                                                                                                                      cATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATg               <  2:1544939/149‑1 (MQ=255)
                                                                                                                                       aTCCGGAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTGAAATTTTAGTCGACCCTGCGCGTACGGTTGTGGCAACCCCCGCAATgg              >  1:2324431/1‑149 (MQ=255)
                                                                                                                                       aTCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATgg              <  2:394947/149‑1 (MQ=255)
                                                                                                                                           ggAACTGACGGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGcgcg          >  2:61457/1‑149 (MQ=255)
                                                                                                                                            gAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCgcgg         >  2:1216407/1‑149 (MQ=255)
                                                                                                                                             aaCTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTAGTCGACCGTGCGCGTCCGGGTGTGGCAACCCCCGAAATGGCgcggc        >  1:2210969/1‑149 (MQ=255)
                                                                                                                                                 gACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTa    >  2:1277533/1‑149 (MQ=255)
                                                                                                                                                 gACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTa    <  1:3339393/149‑1 (MQ=255)
                                                                                                                                                  aCAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGGTGTGGCAACCCGCGAAATGGCGCGGCGGTaa   >  2:1156282/1‑149 (MQ=255)
                                                                                                                                                   cAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAg  <  2:355205/149‑1 (MQ=255)
                                                                                                                                                   cAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAg  <  2:1689713/149‑1 (MQ=255)
                                                                                                                                                   cAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCGTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAg  <  2:4733811/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAG  >  CP009273/1643246‑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: