Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273580,7570GA31.7% 842.9 / inf 584N111N (AAC→AAT) ‡ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (199/200);  new base A (110/75);  total (309/275)
Fisher's exact test for biased strand distribution p-value = 3.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 591 total aligned reads displayed.

AAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCG  >  CP009273/580608‑580904
                                                                                                                                                     |                                                                                                                                                    
 aGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGt                                                                                                                                                     <  2:4448037/150‑1 (MQ=255)
 aGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGt                                                                                                                                                     <  2:2331244/150‑1 (MQ=255)
 aGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAAGTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGt                                                                                                                                                     >  2:618861/1‑150 (MQ=255)
  gCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGtt                                                                                                                                                    >  2:2123571/1‑150 (MQ=255)
    tataACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTAc                                                                                                                                                  <  2:3686168/150‑1 (MQ=255)
      taACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTg                                                                                                                                                >  1:2728842/1‑150 (MQ=255)
      taACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTg                                                                                                                                                <  1:2609247/150‑1 (MQ=255)
       aaCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATTATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTgg                                                                                                                                               >  1:3001224/1‑150 (MQ=255)
       aaCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTgg                                                                                                                                               <  2:1408066/150‑1 (MQ=255)
        aCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCATTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGa                                                                                                                                              <  1:1317931/150‑1 (MQ=255)
          ggCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAat                                                                                                                                            >  2:4851295/1‑150 (MQ=255)
            cactCCTGATAACCTGCCATAATAGCAAGGCGATAATCAGATTCTTTTAAAAACCAGCCTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAatcc                                                                                                                                          <  1:4288454/147‑1 (MQ=255)
               tCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAatccatc                                                                                                                                       >  1:718901/1‑150 (MQ=255)
                 cTGATATCCGGCCATGAGTCCAAGGCGGTACTTGGGTTCGTTGATGAGCCAGCCTTTGATATACAGATCAAATTCGTTGGGATAATTGTGTTGTGTTTCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGTatccatccg                                                                                                                                     >  2:768137/1‑149 (MQ=255)
                     tatCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTcc                                                                                                                                 >  2:1262769/1‑150 (MQ=255)
                        ccgattaTGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGa                                                                                                                              >  2:6483191/7‑150 (MQ=255)
                        ccGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGa                                                                                                                              >  2:4362357/1‑150 (MQ=255)
                         cGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGAt                                                                                                                             <  1:1816263/150‑1 (MQ=255)
                                  gTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCatat                                                                                                                    >  1:102151/1‑150 (MQ=255)
                                   tCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATAtt                                                                                                                   <  2:5335863/150‑1 (MQ=255)
                                   tCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATACGTTGTCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATAtt                                                                                                                   >  2:3266362/1‑150 (MQ=255)
                                             tAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCgg                                                                                                        <  2:3970151/150‑1 (MQ=255)
                                                  gggTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGcc                                                                                                   >  1:6944571/1‑150 (MQ=255)
                                                    gTTCGTTGAGGAGCCAGCCTTTGATAGTCAGATCAAATTCGTTGGCATAATTGAGGTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCCCCTCGGCTGCCga                                                                                                 >  2:6027460/1‑150 (MQ=255)
                                                       cGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGagag                                                                                              <  1:4797458/150‑1 (MQ=255)
                                                            aaaaaccagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTc                                                                                         <  1:1652271/145‑1 (MQ=255)
                                                             ggagCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑TCCACCTCGGCTGCCGAGAGTTGTcc                                                                                        <  2:1707913/150‑1 (MQ=255)
                                                             ggagCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTcc                                                                                        >  1:6634525/1‑150 (MQ=255)
                                                              aaacctgccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTcca                                                                                       <  2:1049551/144‑1 (MQ=255)
                                                                gccagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTccagc                                                                                     >  1:1532892/1‑150 (MQ=255)
                                                                accagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTccagc                                                                                     >  2:2968271/2‑150 (MQ=255)
                                                                accagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTccagc                                                                                     >  1:544803/2‑150 (MQ=255)
                                                                 ccagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTccagcc                                                                                    >  2:1541951/1‑150 (MQ=255)
                                                                  cagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTccagcca                                                                                   >  1:5698366/1‑150 (MQ=255)
                                                                   agccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTccagccag                                                                                  >  1:5009736/1‑150 (MQ=255)
                                                                   agccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGATTGCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACAATATT‑CCCCCCCCGACTATTGAGAGTGGTccagccag                                                                                  >  2:1134304/1‑150 (MQ=255)
                                                                    gccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCcagccagc                                                                                 >  1:5802320/1‑150 (MQ=255)
                                                                       tttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgc                                                                              >  2:6967713/8‑150 (MQ=255)
                                                                         tGATATTCAGATCAAATTCGATGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGccc                                                                            <  1:2460516/150‑1 (MQ=255)
                                                                           acgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcacca                                                                          >  2:3031526/4‑149 (MQ=255)
                                                                            tatTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCga                                                                         <  1:2186576/150‑1 (MQ=255)
                                                                            tatTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCga                                                                         <  1:6058090/150‑1 (MQ=255)
                                                                             gttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcaccaac                                                                        >  2:705922/2‑147 (MQ=255)
                                                                             atTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgat                                                                        <  2:5625417/150‑1 (MQ=255)
                                                                              ttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcaccaaca                                                                       <  2:5941133/150‑5 (MQ=255)
                                                                                caaaTCAAATTCGTTGGCATAATTAAGACGTGTATCTGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacaga                                                                     <  1:3305692/147‑1 (MQ=255)
                                                                                cAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgataga                                                                     <  1:3148607/150‑1 (MQ=255)
                                                                                cAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgataga                                                                     <  2:1486942/150‑1 (MQ=255)
                                                                                 aaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCCGACTAGTGAGAGTGGTCCAGCCAGCTGCACCaacaga                                                                     >  2:4591733/3‑150 (MQ=255)
                                                                                 aaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagat                                                                    <  2:629146/148‑1 (MQ=255)
                                                                                 aaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagat                                                                    >  1:6877912/3‑150 (MQ=255)
                                                                                 aaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagat                                                                    >  1:4528621/3‑150 (MQ=255)
                                                                                 aaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATAGT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagat                                                                    <  1:5555653/148‑1 (MQ=255)
                                                                                  gATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgatagata                                                                   <  1:6137392/150‑1 (MQ=255)
                                                                                    tCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATc                                                                 <  2:185515/150‑1 (MQ=255)
                                                                                      aaaTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGTATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTg                                                               >  1:4870645/1‑150 (MQ=255)
                                                                                      aaaTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTg                                                               >  1:2318820/1‑150 (MQ=255)
                                                                                      aaaTTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataacca                                                               <  1:950800/150‑5 (MQ=255)
                                                                                        aTTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTTTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCCTAGT‑CCCCCCCAGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataaccatg                                                             <  2:3433812/150‑7 (MQ=255)
                                                                                        aTTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataaccatg                                                             >  2:5449438/1‑144 (MQ=255)
                                                                                            gTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCAt                                                         >  1:5097085/1‑150 (MQ=255)
                                                                                               ggCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCaa                                                      >  2:2851292/1‑150 (MQ=255)
                                                                                                gCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCaaa                                                     <  1:6577651/150‑1 (MQ=255)
                                                                                                gCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCaaa                                                     <  1:4724219/150‑1 (MQ=255)
                                                                                                 cATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCaact                                                    <  2:3126445/150‑3 (MQ=255)
                                                                                                  aTACTTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCCCCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATc                                                   >  1:3735288/1‑150 (MQ=255)
                                                                                                  aTAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATc                                                   >  2:2583651/1‑150 (MQ=255)
                                                                                                  aTAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATc                                                   >  1:2088286/1‑150 (MQ=255)
                                                                                                     aTTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCa                                                >  2:2708799/1‑150 (MQ=255)
                                                                                                      ttGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaa                                               >  2:4786699/1‑150 (MQ=255)
                                                                                                         aGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaatta                                            <  1:4189531/150‑1 (MQ=255)
                                                                                                            tgtgTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaattaatt                                         <  1:4412127/150‑1 (MQ=255)
                                                                                                             gtgtATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTg                                        <  1:1494422/150‑1 (MQ=255)
                                                                                                                    aGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCttt                                 >  2:5760513/1‑150 (MQ=255)
                                                                                                                      ggTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTtaa                               <  2:5253578/150‑1 (MQ=255)
                                                                                                                         gtCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTtaataa                            >  2:4090125/1‑150 (MQ=255)
                                                                                                                         gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                            <  1:1225267/145‑1 (MQ=255)
                                                                                                                         gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                            >  2:5664264/6‑150 (MQ=255)
                                                                                                                         gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                            >  1:3441943/6‑150 (MQ=255)
                                                                                                                           ctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatg                          >  2:913055/4‑149 (MQ=255)
                                                                                                                           cTACTTTCATCCGTCCAGGGTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAAtt                          <  2:4258407/150‑1 (MQ=255)
                                                                                                                            tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATT‑CCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                         >  1:3939368/3‑148 (MQ=255)
                                                                                                                              cTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCa                       <  1:2024471/150‑1 (MQ=255)
                                                                                                                              cTTTCATCCGGCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCa                       <  2:6678465/150‑1 (MQ=255)
                                                                                                                                ttCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGc                     >  1:2916795/1‑150 (MQ=255)
                                                                                                                                ttCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGc                     <  2:6851767/150‑1 (MQ=255)
                                                                                                                                ttCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGc                     <  1:3601160/150‑1 (MQ=255)
                                                                                                                                ttCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGc                     >  1:5179505/1‑150 (MQ=255)
                                                                                                                                  cATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGt                   <  1:7066534/150‑1 (MQ=255)
                                                                                                                                    tCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGTGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTa                 >  1:820962/1‑150 (MQ=255)
                                                                                                                                      tgTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCAAATT‑CCCCCCCCGACTATTGAGAGTGGGCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTact               <  2:474136/149‑3 (MQ=255)
                                                                                                                                      cGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTAtt               >  1:7025866/1‑150 (MQ=255)
                                                                                                                                       gTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTg              <  1:6362530/150‑1 (MQ=255)
                                                                                                                                        tCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGa             >  1:5256143/1‑150 (MQ=255)
                                                                                                                                        tCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGa             <  2:6360746/150‑1 (MQ=255)
                                                                                                                                           aGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAAtt          <  2:5123429/150‑1 (MQ=255)
                                                                                                                                            gcttCCGGGGTTACTGGAATCCATCCAGTCCTGCTCGACCATATT‑GCCACCTAGTCAGCAGAGAGTGGTCCAGCCAGCAGCCCAGAGAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTTCAGCGTTATTGAAttt         <  2:1954139/148‑1 (MQ=255)
                                                                                                                                                ccGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAg     <  1:1387255/150‑1 (MQ=255)
                                                                                                                                                 cGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGt    >  2:4989271/1‑150 (MQ=255)
                                                                                                                                                 cGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGt    >  2:3834795/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGCTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATT‑GCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCG  >  CP009273/580608‑580904

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