Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273580,8050GA46.8% 104.7 / inf 507G95G (GGC→GGT) ‡ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (119/151);  new base A (138/99);  total (257/250)
Fisher's exact test for biased strand distribution p-value = 1.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 512 total aligned reads displayed.

TTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTT  >  CP009273/580656‑580954
                                                                                                                                                     |                                                                                                                                                     
ttGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTg                                                                                                                                                       <  1:2342197/150‑1 (MQ=255)
ttGGGTTCGTGGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTg                                                                                                                                                       >  1:1552775/1‑150 (MQ=255)
        gTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGt                                                                                                                                               <  1:5109310/150‑1 (MQ=255)
              aaaccagacTTTAACGTTCAAATCAAATTCGTTGTCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTcca                                                                                                                                         <  2:2833843/147‑1 (MQ=255)
                accagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTccagc                                                                                                                                       <  1:5988461/149‑1 (MQ=255)
                accagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTccagc                                                                                                                                       >  1:434731/2‑150 (MQ=255)
                accagccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTccagc                                                                                                                                       >  1:4976400/2‑150 (MQ=255)
                 ccagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTccagcc                                                                                                                                      >  1:1342230/1‑150 (MQ=255)
                  cagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTccagcca                                                                                                                                     >  2:665249/1‑150 (MQ=255)
                  cagccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTccagcca                                                                                                                                     <  2:1320905/150‑1 (MQ=255)
                    gccTTTAACGTTCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCcagccagc                                                                                                                                   >  2:3040094/1‑150 (MQ=255)
                     ccTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCcagca                                                                                                                                  >  2:322957/1‑150 (MQ=255)
                       tttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgc                                                                                                                                >  2:3914377/8‑150 (MQ=255)
                       tttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgc                                                                                                                                >  1:620435/8‑150 (MQ=255)
                        ttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgca                                                                                                                               <  1:3836030/144‑2 (MQ=255)
                        ttaacgttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCagctgca                                                                                                                               >  2:4903939/7‑149 (MQ=255)
                          gATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGcccc                                                                                                                             >  1:1030280/1‑150 (MQ=255)
                           atatTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCg                                                                                                                            >  2:5685445/1‑150 (MQ=255)
                             atTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgat                                                                                                                          >  1:87587/1‑150 (MQ=255)
                              ttCAAATCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcaccaaca                                                                                                                         <  1:2807430/150‑5 (MQ=255)
                               tcaaatcaaaTTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGcaccaacag                                                                                                                        >  1:3747303/1‑145 (MQ=255)
                                caaaTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacaga                                                                                                                       >  2:5073802/4‑150 (MQ=255)
                                 aGATCAAATTCGTTGGCATAATTGAGGTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgatagat                                                                                                                      >  2:1475206/1‑150 (MQ=255)
                                  gATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCgatagata                                                                                                                     <  2:2970808/150‑1 (MQ=255)
                                   aTCAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataa                                                                                                                    >  1:1593502/1‑149 (MQ=255)
                                     cAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCt                                                                                                                  <  1:4135224/150‑1 (MQ=255)
                                     cAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCt                                                                                                                  <  2:6061198/150‑1 (MQ=255)
                                     cAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCt                                                                                                                  <  1:3241474/150‑1 (MQ=255)
                                     cAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataacc                                                                                                                  >  1:2580867/1‑147 (MQ=255)
                                     cAAATTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataacc                                                                                                                  >  1:4836972/1‑147 (MQ=255)
                                       aaTTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataaccat                                                                                                                >  1:3382837/1‑145 (MQ=255)
                                        aTTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTggg                                                                                                               <  1:2139047/150‑1 (MQ=255)
                                        aTTCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataaccatg                                                                                                               >  2:3380410/1‑144 (MQ=255)
                                         ttCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTgggg                                                                                                              <  2:2437128/150‑1 (MQ=255)
                                         ttCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTgggg                                                                                                              >  2:30975/1‑150 (MQ=255)
                                         ttCGTTGGCATAATTAAGACGTGTATCAGGATGCCTGCTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCaacagataaccatgg                                                                                                              <  2:4432039/150‑8 (MQ=255)
                                             ttGGCAAAAGTGAGTAGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTCCTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGGCCAGCCAGCAGCCCCGATAGATATCTGGGGCATc                                                                                                          <  2:1880586/150‑1 (MQ=255)
                                                 cATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAAt                                                                                                      <  2:2613439/150‑1 (MQ=255)
                                                   tAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATcc                                                                                                    >  1:5331552/1‑150 (MQ=255)
                                                    aaTTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATccc                                                                                                   >  1:3001670/1‑150 (MQ=255)
                                                     aTTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCa                                                                                                  <  1:5726561/150‑1 (MQ=255)
                                                       tGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaat                                                                                                <  1:753809/150‑1 (MQ=255)
                                                       tGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaat                                                                                                >  1:2397025/1‑150 (MQ=255)
                                                           ttGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCaattaat                                                                                            <  1:1433976/150‑1 (MQ=255)
                                                                   cAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCtt                                                                                    <  2:1806001/150‑1 (MQ=255)
                                                                         gtCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTtaataa                                                                              <  2:427849/150‑1 (MQ=255)
                                                                         gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                                                              <  1:2923376/145‑1 (MQ=255)
                                                                         gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                                                              <  2:1711925/145‑1 (MQ=255)
                                                                         gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                                                              <  1:4475237/145‑1 (MQ=255)
                                                                         gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                                                              <  1:938710/145‑1 (MQ=255)
                                                                         gcctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTtaagaa                                                                              >  2:5166507/6‑150 (MQ=255)
                                                                          tCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTaataat                                                                             >  1:5486386/1‑150 (MQ=255)
                                                                          cctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaat                                                                             <  1:2051023/146‑1 (MQ=255)
                                                                          cctgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaat                                                                             <  2:4531266/146‑1 (MQ=255)
                                                                            tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                                                                           <  1:3573118/148‑3 (MQ=255)
                                                                            tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                                                                           >  1:3410407/3‑148 (MQ=255)
                                                                            tgcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTaagaatgg                                                                           <  2:2702650/148‑3 (MQ=255)
                                                                             gcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGc                                                                          >  2:5796926/2‑150 (MQ=255)
                                                                             gcTTTCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGc                                                                          >  2:529589/2‑150 (MQ=255)
                                                                               tttCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGcgg                                                                        >  2:5974921/1‑148 (MQ=255)
                                                                                ttCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATggcggc                                                                       >  2:4149883/1‑150 (MQ=255)
                                                                                 tCATCTGTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGgcggcg                                                                      <  1:2252672/150‑1 (MQ=255)
                                                                                  cATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGt                                                                     <  1:6072856/150‑1 (MQ=255)
                                                                                    tctgTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTa                                                                   >  2:3287957/4‑150 (MQ=255)
                                                                                    tctgTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTa                                                                   <  1:4668776/147‑1 (MQ=255)
                                                                                       gTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTg                                                                >  1:1970/1‑150 (MQ=255)
                                                                                       gTCCATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTg                                                                >  1:3882817/1‑150 (MQ=255)
                                                                                        tCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGa                                                               >  2:780358/1‑150 (MQ=255)
                                                                                        tCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGa                                                               >  1:6122985/1‑150 (MQ=255)
                                                                                        tCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTCAATAATTGCAGCGTTATTGa                                                               >  1:3170836/1‑150 (MQ=255)
                                                                                         ccATGTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTgta                                                              <  1:5635858/150‑3 (MQ=255)
                                                                                         ccAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGaa                                                              >  2:5538458/1‑150 (MQ=255)
                                                                                           aGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAAtt                                                            >  1:4518486/1‑150 (MQ=255)
                                                                                            tgTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTCACGGCGCCTTTAAGAATGGCGGCGTTACTGTAttt                                                           >  2:6136256/2‑150 (MQ=255)
                                                                                            tgTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTAttt                                                           <  2:1169338/149‑1 (MQ=255)
                                                                                             gTTCCGGGAGTCCCGGAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTc                                                          <  1:4572914/150‑1 (MQ=255)
                                                                                              ttCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTcc                                                         >  2:402096/1‑150 (MQ=255)
                                                                                               tCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCa                                                        <  1:2693662/150‑1 (MQ=255)
                                                                                                 cGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGt                                                      <  1:3316864/150‑1 (MQ=255)
                                                                                                     attcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTtccagtccag                                                  <  1:5801411/144‑1 (MQ=255)
                                                                                                     agtcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCACCAGATAACCATGGATTCAGCTCCCAGTTATCGGCGCCTTTAAGAATGGCGGCGTTACTGTATTtccagtccag                                                  <  2:3779492/144‑1 (MQ=255)
                                                                                                      ttACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGt                                                 >  1:1747125/1‑150 (MQ=255)
                                                                                                       tcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTCCAGTCCAGtt                                                >  2:3613927/5‑150 (MQ=255)
                                                                                                       tcccggAATCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACTATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTCCAGTCCAGtt                                                >  2:227279/5‑150 (MQ=255)
                                                                                                          tGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTtgac                                             <  1:1824280/150‑1 (MQ=255)
                                                                                                          tGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTtgac                                             <  2:4352355/150‑1 (MQ=255)
                                                                                                             aaTCCATCCAGTCCTGATCAACCATATTCCCCCCCCGACGATTGAGAGTGGTCCAGCCAGCTGCACCAACAGATAACCATGGATTCAGCTCCCAGTTAACGGCGCCTTTAAGAATGGCGGCGTTACTGTATTTCCAGTCCAGTtggctga                                          >  1:2014530/1‑150 (MQ=255)
                                                                                                               tccatccaGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTgactgact                                        <  2:2319283/150‑1 (MQ=255)
                                                                                                                    ccaGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCg                                   <  1:4231311/150‑1 (MQ=255)
                                                                                                                       gTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGcc                                >  2:2558180/1‑150 (MQ=255)
                                                                                                                        tCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGcct                               >  2:1080108/1‑150 (MQ=255)
                                                                                                                             aTCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCcttc                          >  1:1238847/1‑150 (MQ=255)
                                                                                                                                    atatTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCt                   >  1:1911391/1‑150 (MQ=255)
                                                                                                                                       ttGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGa                <  2:3534498/150‑1 (MQ=255)
                                                                                                                                          ccACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATaa             <  1:2941440/150‑1 (MQ=255)
                                                                                                                                           cACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATaaa            <  1:5236693/150‑1 (MQ=255)
                                                                                                                                           cACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATaaa            <  2:2096721/150‑1 (MQ=255)
                                                                                                                                             ccTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAaca          <  1:1141820/150‑1 (MQ=255)
                                                                                                                                              cTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAacac         >  1:3105984/1‑150 (MQ=255)
                                                                                                                                                    tGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCtt   >  1:3644589/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGGGTTCGTTGAGGAGCCAGCCTTTGATATTCAGATCAAATTCGTTGGCATAATTGAGTTGTGTATCAGGGTGTCTACTTTCATCCGTCCAGGTTCCGGGGTTACTGGAATCCATCCAGTCCTGATCGACCATATTGCCACCTCGGCTGCCGAGAGTTGTCCAGCCAGCAGCCCCGATAGATATCTGGGGCATCAAATCCCAATTAATTGCACCTTTAATAATTGCAGCGTTATTGAATTTCCAGTCGAGTTGACTGACTTTTCGGCCTCCTTCTTCGGCTAGATAAACACGCTCTTT  >  CP009273/580656‑580954

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