Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273580,7600AC39.1% 483.3 / inf 548S110R (AGT→AGG) ‡ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base A (168/165);  new base C (124/90);  total (292/256)
Fisher's exact test for biased strand distribution p-value = 9.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 552 total aligned reads displayed.

CTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTT  >  CP009273/580611‑580908
                                                                                                                                                     |                                                                                                                                                    
cTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTa                                                                                                                                                      >  1:5540191/1‑150 (MQ=255)
cTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTa                                                                                                                                                      >  1:1487792/1‑150 (MQ=255)
 tataACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTAc                                                                                                                                                     >  2:1280267/1‑150 (MQ=255)
 tataACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTAc                                                                                                                                                     >  2:1997089/1‑150 (MQ=255)
    aaCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTgg                                                                                                                                                  >  2:1391830/1‑150 (MQ=255)
    aaCGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTgg                                                                                                                                                  >  1:2579914/1‑150 (MQ=255)
      cGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGaa                                                                                                                                                <  1:1082093/150‑1 (MQ=255)
       ggCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAat                                                                                                                                               >  1:4242865/1‑150 (MQ=255)
         cTTTCCTGATAACCTGCCATAATAGCAAGGCGATAATCAGATTCTTTTAAAAACCAGCCTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAatcc                                                                                                                                             <  1:4168009/150‑1 (MQ=255)
          tttCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAatcca                                                                                                                                            <  2:3092078/150‑1 (MQ=255)
           ttCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAatccat                                                                                                                                           >  2:6407025/1‑150 (MQ=255)
           ttCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAatccat                                                                                                                                           >  2:2913911/1‑150 (MQ=255)
            tCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAatccatc                                                                                                                                          >  2:6559019/1‑150 (MQ=255)
              cTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAAtccatcca                                                                                                                                        >  1:1631300/1‑150 (MQ=255)
              cTGATAACCTGCCATAAGAGCAAGGCGATAATCAGATTCTTTTAAAAACCAGCCTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTACGGGAGTCCCGGAAtccatcca                                                                                                                                        <  2:256420/150‑1 (MQ=255)
               tGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAg                                                                                                                                       >  2:461449/1‑150 (MQ=255)
               tGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAg                                                                                                                                       <  1:2450783/150‑1 (MQ=255)
                   atCCGGCCATGAGTCCCAGGCGGTAATTGCGTTCGTTGAGGAGCCAGCCATTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCt                                                                                                                                   <  2:5934894/150‑1 (MQ=255)
                    tCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTg                                                                                                                                  >  1:1811642/1‑150 (MQ=255)
                         ccATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGa                                                                                                                             <  2:6209537/150‑1 (MQ=255)
                               gTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCatat                                                                                                                       <  2:5342350/150‑1 (MQ=255)
                                tCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATAtt                                                                                                                      <  1:2213146/150‑1 (MQ=255)
                                 cccAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTg                                                                                                                     >  2:3756777/1‑150 (MQ=255)
                                     ggCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCAc                                                                                                                 >  1:578274/1‑150 (MQ=255)
                                     ggCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCAc                                                                                                                 >  2:700440/1‑150 (MQ=255)
                                        ggTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTc                                                                                                              >  2:216561/1‑150 (MQ=255)
                                               gggTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGcc                                                                                                       >  2:4113906/1‑150 (MQ=255)
                                                 gTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCga                                                                                                     <  1:3302491/150‑1 (MQ=255)
                                                  ttCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCgag                                                                                                    >  2:6312637/1‑150 (MQ=255)
                                                     gTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGt                                                                                                 <  2:5826364/150‑1 (MQ=255)
                                                     gTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGt                                                                                                 <  1:921808/150‑1 (MQ=255)
                                                     gTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGt                                                                                                 <  2:2053677/150‑1 (MQ=255)
                                                      ttGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGtt                                                                                                >  2:6022802/1‑150 (MQ=255)
                                                         aggagCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTc                                                                                             >  1:3625753/1‑150 (MQ=255)
                                                            aaccagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTccag                                                                                          >  2:1810801/3‑150 (MQ=255)
                                                               cagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTccagcca                                                                                       >  2:6754963/1‑150 (MQ=255)
                                                               cagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTccagcca                                                                                       <  1:6574037/150‑1 (MQ=255)
                                                                agccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTccagccag                                                                                      >  1:6418375/1‑150 (MQ=255)
                                                                 gccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCcagccagc                                                                                     <  1:2927161/150‑1 (MQ=255)
                                                                    tttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgc                                                                                  >  1:4661091/8‑150 (MQ=255)
                                                                     ttaactttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCTGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgca                                                                                 <  1:3232201/144‑2 (MQ=255)
                                                                     ttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgca                                                                                 >  1:766219/7‑149 (MQ=255)
                                                                     ttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCACCCCCCGACTATTGTGAGTGGTCCAGCCagctgca                                                                                 >  2:920847/7‑149 (MQ=255)
                                                                     ttaacgttCAAATCAAATTCGGTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACGATTGAGAGTGGTCCAGCCagctgca                                                                                 >  1:6166774/7‑149 (MQ=255)
                                                                     ttGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGcc                                                                                 >  1:828534/1‑150 (MQ=255)
                                                                      tGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGccc                                                                                <  1:5644276/150‑1 (MQ=255)
                                                                         tatTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCga                                                                             >  2:3048140/1‑150 (MQ=255)
                                                                          gttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcaccaac                                                                            >  1:1719754/2‑147 (MQ=255)
                                                                          gttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcaccaac                                                                            >  2:1127923/2‑147 (MQ=255)
                                                                          atTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgat                                                                            <  1:3423089/150‑1 (MQ=255)
                                                                           ttCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgata                                                                           <  1:2763752/150‑1 (MQ=255)
                                                                           ttCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgata                                                                           <  2:5611963/150‑1 (MQ=255)
                                                                           ttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcaccaaca                                                                           <  1:1167588/150‑5 (MQ=255)
                                                                            tCATATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCAATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgatag                                                                          <  2:225983/150‑1 (MQ=255)
                                                                              aaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagat                                                                        <  2:1469736/148‑1 (MQ=255)
                                                                              aGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgatagat                                                                        >  1:5048053/1‑150 (MQ=255)
                                                                              aGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgatagat                                                                        >  1:6313231/1‑150 (MQ=255)
                                                                                aTCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGatat                                                                      <  1:4545287/150‑1 (MQ=255)
                                                                                 tCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataac                                                                     >  2:4352008/1‑148 (MQ=255)
                                                                                      ttCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTgggg                                                                <  1:4905962/150‑1 (MQ=255)
                                                                                         gTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCAt                                                             <  1:1037819/150‑1 (MQ=255)
                                                                                          ttGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCATCCCCGATAGATATCTGGGGCATc                                                            >  2:1539259/1‑150 (MQ=255)
                                                                                            ggCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCCGCAGCCCCGATAGATATCTGGGGCATCaa                                                          >  2:4794615/1‑150 (MQ=255)
                                                                                               aTAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATc                                                       <  2:6526518/150‑1 (MQ=255)
                                                                                                   ttGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaa                                                   <  2:39690/150‑1 (MQ=255)
                                                                                                            ctATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGAGATCTGGGGCATCAAATCCCAATTAATTGCa                                          <  2:2657040/149‑1 (MQ=255)
                                                                                                              aTCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCAcc                                        >  2:248245/1‑150 (MQ=255)
                                                                                                                      gtCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTtaataa                                <  1:1945342/150‑1 (MQ=255)
                                                                                                                      gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                <  1:5968780/145‑1 (MQ=255)
                                                                                                                      gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                <  1:5496627/145‑1 (MQ=255)
                                                                                                                      gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                <  1:6651978/145‑1 (MQ=255)
                                                                                                                      gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                <  1:3541314/145‑1 (MQ=255)
                                                                                                                      gcctgcTTTAATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                <  2:1053773/145‑1 (MQ=255)
                                                                                                                      gacTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTtaataa                                >  1:2121158/3‑150 (MQ=255)
                                                                                                                         tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                             <  1:2315250/148‑3 (MQ=255)
                                                                                                                         tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                             >  2:5186324/3‑148 (MQ=255)
                                                                                                                         tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCACCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                             >  1:3804329/3‑148 (MQ=255)
                                                                                                                         tACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTg                             >  1:4412714/1‑150 (MQ=255)
                                                                                                                           cTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGcg                           <  1:1277675/150‑2 (MQ=255)
                                                                                                                           cTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCa                           <  2:4047199/150‑1 (MQ=255)
                                                                                                                            tttCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGcgg                          >  1:3073557/1‑148 (MQ=255)
                                                                                                                            tttCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAg                          >  1:3946742/1‑150 (MQ=255)
                                                                                                                             ttCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATggcggc                         >  1:5280317/1‑150 (MQ=255)
                                                                                                                              tCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTAGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCg                        >  1:4339701/1‑150 (MQ=255)
                                                                                                                                aTCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGtt                      >  1:367747/1‑150 (MQ=255)
                                                                                                                                 tctgTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTa                     >  2:1582259/4‑150 (MQ=255)
                                                                                                                                 tctgTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTa                     >  2:5386970/4‑150 (MQ=255)
                                                                                                                                   cGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTAtt                   <  1:4785533/150‑1 (MQ=255)
                                                                                                                                     tCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTcc                                                                                            >  1:6346004/1‑75 (MQ=255)
                                                                                                                                     tCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGa                 >  1:1858980/1‑150 (MQ=255)
                                                                                                                                      ccAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATAAATTGCACCTTTAATACTTGCAGCGTTATTGat                >  2:3159270/1‑149 (MQ=255)
                                                                                                                                        atgTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTAtt              >  2:5732944/3‑150 (MQ=255)
                                                                                                                                        aGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAAtt              <  2:5030821/150‑1 (MQ=255)
                                                                                                                                        aGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAAtt              <  2:1625304/150‑1 (MQ=255)
                                                                                                                                           ttCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTcc           <  2:2343658/150‑1 (MQ=255)
                                                                                                                                             ccGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAg         >  2:3569720/1‑150 (MQ=255)
                                                                                                                                             ccGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAg         <  2:3349122/150‑1 (MQ=255)
                                                                                                                                                gggTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCg      >  1:510835/1‑150 (MQ=255)
                                                                                                                                                 ggTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGAGATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGa     >  2:4923078/1‑150 (MQ=255)
                                                                                                                                                   ttACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGt   >  1:5725404/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTATAACGGCTTTCCTGATATCCGGCCATGAGTCCCAGGCGGTAATTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTT  >  CP009273/580611‑580908

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: