Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273580,8140TG44.3% 260.3 / inf 535T92T (ACA→ACC) ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base T (138/160);  new base G (138/99);  total (276/259)
Fisher's exact test for biased strand distribution p-value = 6.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 538 total aligned reads displayed.

GAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCC  >  CP009273/580667‑580963
                                                                                                                                                   |                                                                                                                                                     
gaggagCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGt                                                                                                                                                     <  2:3501508/150‑1 (MQ=255)
 aaaaaccagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTc                                                                                                                                                    <  1:487741/145‑1 (MQ=255)
  aaaaccagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGACCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTcc                                                                                                                                                   <  1:2056613/146‑1 (MQ=255)
    agccagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCCGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTccag                                                                                                                                                 >  1:6160579/1‑150 (MQ=255)
      ccagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTccagcc                                                                                                                                               >  2:2871366/1‑150 (MQ=255)
        agccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTccagccag                                                                                                                                             >  2:789447/1‑150 (MQ=255)
         gccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCcagccagc                                                                                                                                            <  2:5292874/150‑1 (MQ=255)
          ccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCcagca                                                                                                                                           >  1:2425056/1‑150 (MQ=255)
          ccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCcagca                                                                                                                                           <  2:5385029/150‑1 (MQ=255)
           cTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCcagcag                                                                                                                                          >  2:2109541/1‑150 (MQ=255)
            tttGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCagcagc                                                                                                                                         >  2:1112699/1‑150 (MQ=255)
             ttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTTAGAGTGGTCCAGCCagctgca                                                                                                                                        >  2:578425/7‑149 (MQ=255)
             ttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgca                                                                                                                                        <  1:3836030/144‑2 (MQ=255)
             ttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgca                                                                                                                                        >  2:4903939/7‑149 (MQ=255)
             ttGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGcc                                                                                                                                        <  2:3089442/150‑1 (MQ=255)
              taacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgcac                                                                                                                                       >  1:3929118/6‑148 (MQ=255)
               gATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGcccc                                                                                                                                      >  1:1030280/1‑150 (MQ=255)
                 tatTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCga                                                                                                                                    >  2:2455713/1‑150 (MQ=255)
                    tcaaatcaaaTTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcaccaacag                                                                                                                                 >  1:3747303/1‑145 (MQ=255)
                     caaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacaga                                                                                                                                <  1:4794034/147‑1 (MQ=255)
                     caaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacaga                                                                                                                                >  2:5073802/4‑150 (MQ=255)
                     caaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacaga                                                                                                                                <  2:2731886/147‑1 (MQ=255)
                      aaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagat                                                                                                                               >  1:779364/3‑150 (MQ=255)
                      aGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgatagat                                                                                                                               >  1:5134891/1‑150 (MQ=255)
                        aTCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGatat                                                                                                                             <  1:2757828/150‑1 (MQ=255)
                          cAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCt                                                                                                                           <  1:4135224/150‑1 (MQ=255)
                          cAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCt                                                                                                                           <  1:3241443/150‑1 (MQ=255)
                          cAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCt                                                                                                                           <  2:6061198/150‑1 (MQ=255)
                           aaaTTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataacca                                                                                                                          <  1:3042969/150‑5 (MQ=255)
                             aTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTggg                                                                                                                        <  1:5753959/150‑1 (MQ=255)
                             aTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTggg                                                                                                                        >  1:225576/1‑150 (MQ=255)
                             aTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTggg                                                                                                                        <  1:2139047/150‑1 (MQ=255)
                              ttCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataaccatgg                                                                                                                       <  2:4432039/150‑8 (MQ=255)
                                   tGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCa                                                                                                                  >  2:2283674/1‑150 (MQ=255)
                                    ggCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCaa                                                                                                                 <  1:296839/150‑1 (MQ=255)
                                        tAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATcc                                                                                                             >  1:5332763/1‑150 (MQ=255)
                                        tAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATcc                                                                                                             >  1:656537/1‑150 (MQ=255)
                                        tAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATcc                                                                                                             >  2:1837709/1‑150 (MQ=255)
                                         aaTTGAGTTGTGTATCAGGGTGACTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATccc                                                                                                            >  2:3253278/1‑150 (MQ=255)
                                           ttGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaa                                                                                                          >  1:1562733/1‑150 (MQ=255)
                                                   tgtATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGc                                                                                                  <  1:6061937/150‑1 (MQ=255)
                                                    gtATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCa                                                                                                 <  1:3556898/150‑1 (MQ=255)
                                                              gtCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAACTCCCAATTAATTGCACCTTtaataa                                                                                       <  2:4273586/150‑1 (MQ=255)
                                                              gcctgcTTTCATCTGTCCATGGTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATCGTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                                                                       <  2:2662337/145‑1 (MQ=255)
                                                               tCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTaataat                                                                                      >  1:5486386/1‑150 (MQ=255)
                                                               tCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTaataat                                                                                      <  2:1562421/150‑1 (MQ=255)
                                                               actgcTTTCATCTGTCCATGTTCCGGTAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGAATATTGAGAGTGTTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaat                                                                                      <  2:4531419/146‑1 (MQ=255)
                                                                 tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                                                                                    >  2:717596/3‑148 (MQ=255)
                                                                 tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                                                                                    <  1:3403442/148‑3 (MQ=255)
                                                                  gcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGc                                                                                   <  1:2016544/149‑1 (MQ=255)
                                                                  gcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGc                                                                                   <  1:3166019/149‑1 (MQ=255)
                                                                  gcTTTAATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGc                                                                                   <  2:3020848/149‑1 (MQ=255)
                                                                  aCTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGc                                                                                   <  1:4459077/150‑1 (MQ=255)
                                                                   cTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGcg                                                                                  <  2:1907269/150‑2 (MQ=255)
                                                                    tttCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGcgg                                                                                 >  2:1021526/1‑148 (MQ=255)
                                                                    tttCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAg                                                                                 <  1:1347951/150‑1 (MQ=255)
                                                                     ttCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATggcggc                                                                                >  2:4149883/1‑150 (MQ=255)
                                                                       cATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGt                                                                              <  1:93809/150‑1 (MQ=255)
                                                                        aTCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGtt                                                                             >  1:2381882/1‑150 (MQ=255)
                                                                          ctgTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTac                                                                           >  2:1738310/3‑149 (MQ=255)
                                                                           tgTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTact                                                                          <  2:3225391/149‑3 (MQ=255)
                                                                            gTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTg                                                                         >  1:1970/1‑150 (MQ=255)
                                                                            gTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTg                                                                         >  1:1889366/1‑150 (MQ=255)
                                                                              ccATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTgta                                                                       <  1:5635858/150‑3 (MQ=255)
                                                                               cAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAAt                                                                      >  2:5921474/1‑150 (MQ=255)
                                                                                atgTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTAtt                                                                     >  1:5965925/3‑150 (MQ=255)
                                                                                 tgTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTCACGGCGCCTTTAAGAATGGCGGCGTTACTGTAttt                                                                    >  2:6136256/2‑150 (MQ=255)
                                                                                  gTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGTCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTc                                                                   <  1:2645319/150‑1 (MQ=255)
                                                                                   ttCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTcc                                                                  >  2:402096/1‑150 (MQ=255)
                                                                                    tCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCa                                                                 >  2:3868713/1‑150 (MQ=255)
                                                                                          agtcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCACCAGATAACCATGGATTCAGCTCCCAGTTATCGGCGCCTTTAAGAATGGCGGCGTTACTGTATTtccagtccag                                                           <  2:3779492/144‑1 (MQ=255)
                                                                                          agtcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTtccagtccag                                                           <  1:2925242/144‑1 (MQ=255)
                                                                                           ttACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGt                                                          >  1:1747125/1‑150 (MQ=255)
                                                                                           gtcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTccagtccagt                                                          <  1:6095526/145‑1 (MQ=255)
                                                                                            tcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTCCAGTCCAGtt                                                         >  2:2355264/5‑150 (MQ=255)
                                                                                            tcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTCCAGTCCAGtt                                                         >  2:3613927/5‑150 (MQ=255)
                                                                                              ccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTCCAGTCCAGTtgg                                                       >  1:3342485/3‑149 (MQ=255)
                                                                                               tGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTtgac                                                      <  2:4352355/150‑1 (MQ=255)
                                                                                               tGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTtgac                                                      <  1:1824280/150‑1 (MQ=255)
                                                                                               cggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTCCAGTCCAGTtggc                                                      <  2:255937/149‑3 (MQ=255)
                                                                                                 gAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGGGGCGTTACTGTATTTCCAGTCCAGTtggctg                                                    >  1:2872415/1‑150 (MQ=255)
                                                                                                        tccaGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTc                                             >  2:3376911/1‑150 (MQ=255)
                                                                                                         ccaGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCg                                            <  1:4231311/150‑1 (MQ=255)
                                                                                                           aGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGc                                          >  2:1347659/1‑150 (MQ=255)
                                                                                                            gTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGcc                                         >  2:5620567/1‑150 (MQ=255)
                                                                                                                  aTCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCcttc                                   >  1:1238847/1‑150 (MQ=255)
                                                                                                                    cGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCcttctt                                 >  1:4506132/1‑150 (MQ=255)
                                                                                                                            ttGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGa                         <  2:1273172/150‑1 (MQ=255)
                                                                                                                                  ccTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAaca                   <  1:1141820/150‑1 (MQ=255)
                                                                                                                                      ggCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGct               >  1:4565271/1‑150 (MQ=255)
                                                                                                                                         tGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCtt            >  1:3644589/1‑150 (MQ=255)
                                                                                                                                           ccGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCtttt          <  1:5032554/150‑1 (MQ=255)
                                                                                                                                           ccGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCtttt          >  1:1647415/1‑150 (MQ=255)
                                                                                                                                            cGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTg         >  1:5829796/1‑150 (MQ=255)
                                                                                                                                             gagaGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGt        >  1:4641231/1‑150 (MQ=255)
                                                                                                                                             gagaGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGt        <  1:4315543/150‑1 (MQ=255)
                                                                                                                                              agagTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGtt       <  1:1443836/150‑1 (MQ=255)
                                                                                                                                                 gTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATACTTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGttttt    <  1:3302037/150‑1 (MQ=255)
                                                                                                                                                   tGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTcc  <  1:5247451/150‑1 (MQ=255)
                                                                                                                                                   tGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTcc  <  2:5207377/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTTTGTTTTTCC  >  CP009273/580667‑580963

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: ‑