Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,430,873 G→T 100% S93* (TCA→TAA)  ompN ← outer membrane pore protein N, non‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,430,8730GT98.5% 863.9 / 3.0 268S93* (TCA→TAA) ompNouter membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (137/126);  total (141/127)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 268 total aligned reads displayed.

TTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAG  >  CP009273/1430726‑1431017
                                                                                                                                                   |                                                                                                                                                
ttGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTAc                                                                                                                                                >  1:2387298/1‑150 (MQ=255)
 tGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACt                                                                                                                                               >  1:3652864/1‑150 (MQ=255)
  gTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTg                                                                                                                                              >  2:4174385/1‑150 (MQ=255)
     tGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGtt                                                                                                                                           <  2:905391/150‑1 (MQ=255)
     tGCATGGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGtt                                                                                                                                           <  2:135457/150‑1 (MQ=255)
         ttGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGtttttt                                                                                                                                       <  1:1030752/150‑1 (MQ=255)
           ggTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGa                                                                                                                                     <  1:3402231/150‑1 (MQ=255)
            gTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGaa                                                                                                                                    >  1:4632398/1‑150 (MQ=255)
                aaGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGAtt                                                                                                                                <  2:975587/150‑1 (MQ=255)
                aaGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGAtt                                                                                                                                <  2:1780643/150‑1 (MQ=255)
                aaGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGAtt                                                                                                                                >  1:4231839/1‑150 (MQ=255)
                    gTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAgt                                                                                                                            <  1:4123806/150‑1 (MQ=255)
                     tCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAgtg                                                                                                                           <  1:2191695/150‑1 (MQ=255)
                       aCCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGtgtt                                                                                                                         >  2:2885691/1‑150 (MQ=255)
                       aCCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGtgtt                                                                                                                         <  1:5018045/150‑1 (MQ=255)
                       aCCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGtgtt                                                                                                                         <  2:3222533/150‑1 (MQ=255)
                        ccGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGtgttg                                                                                                                        >  2:4383708/1‑150 (MQ=255)
                          gCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTgttgtt                                                                                                                      <  1:3194617/150‑1 (MQ=255)
                           ccAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGttt                                                                                                                     <  2:3844799/150‑1 (MQ=255)
                             aaaTTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGc                                                                                                                   >  2:1100521/1‑150 (MQ=255)
                              aaTTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGcc                                                                                                                  <  2:1757936/150‑1 (MQ=255)
                                   aGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAat                                                                                                             <  2:4911870/150‑1 (MQ=255)
                                    ggCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACTCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAata                                                                                                            <  2:313677/150‑1 (MQ=255)
                                      cagcaTATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATAtt                                                                                                          >  1:1829019/1‑150 (MQ=255)
                                      cagcaTATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATAtt                                                                                                          <  1:3558340/150‑1 (MQ=255)
                                        gcaTATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGt                                                                                                        >  2:1621751/1‑150 (MQ=255)
                                          atatCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTAt                                                                                                      >  2:569734/1‑150 (MQ=255)
                                            atCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTc                                                                                                    >  2:3196005/1‑150 (MQ=255)
                                             tCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTcc                                                                                                   <  1:291251/150‑1 (MQ=255)
                                             tCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTcc                                                                                                   <  1:2582330/150‑1 (MQ=255)
                                             tCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTcc                                                                                                   >  2:214970/1‑150 (MQ=255)
                                                  ccAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTg                                                                                              <  2:1510503/150‑1 (MQ=255)
                                                    aGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGAc                                                                                            >  2:2314189/1‑150 (MQ=255)
                                                    aGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGAc                                                                                            >  2:1032876/1‑150 (MQ=255)
                                                    aGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGAc                                                                                            <  2:1929267/150‑1 (MQ=255)
                                                    aGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGAc                                                                                            >  2:4487128/1‑150 (MQ=255)
                                                     gCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGAcc                                                                                           >  1:1165270/1‑150 (MQ=255)
                                                     gCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGAcc                                                                                           <  1:3142283/150‑1 (MQ=255)
                                                        ttCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTa                                                                                        <  1:4379600/150‑1 (MQ=255)
                                                        ttCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTa                                                                                        >  2:4813719/1‑150 (MQ=255)
                                                        ttCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTa                                                                                        <  2:1231295/150‑1 (MQ=255)
                                                        ttCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTa                                                                                        >  1:628358/1‑150 (MQ=255)
                                                           gATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGcc                                                                                     <  2:413443/150‑1 (MQ=255)
                                                            aTGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCa                                                                                    <  2:5033986/150‑1 (MQ=255)
                                                            aTGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGATTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTATCCa                                                                                    <  1:907045/150‑1 (MQ=255)
                                                             tGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAg                                                                                   >  1:1515434/1‑150 (MQ=255)
                                                                   aCATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGtt                                                                             >  2:485658/1‑150 (MQ=255)
                                                                   aCATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGtt                                                                             <  1:4212770/150‑1 (MQ=255)
                                                                      tgaCGCCATAATTACGTCCGTAATCGAAAGAACCGGAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGAt                                                                          <  2:2931784/148‑1 (MQ=255)
                                                                             aTAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAAt                                                                   >  2:3652470/1‑150 (MQ=255)
                                                                              tAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAAtt                                                                  <  1:5085799/150‑1 (MQ=255)
                                                                               aaTTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAAttt                                                                 >  2:1943977/1‑150 (MQ=255)
                                                                                    cGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGtt                                                            <  2:3325875/150‑1 (MQ=255)
                                                                                     gTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGttt                                                           >  1:1646004/1‑150 (MQ=255)
                                                                                      tCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTc                                                          <  2:4073094/150‑1 (MQ=255)
                                                                                      tCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTc                                                          <  2:2426577/150‑1 (MQ=255)
                                                                                         gTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGcc                                                       <  1:1702165/150‑1 (MQ=255)
                                                                                              cGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTna                                                  <  1:2703121/150‑1 (MQ=255)
                                                                                                aaaGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTaaaa                                                >  1:3448078/1‑150 (MQ=255)
                                                                                                 aaGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAAc                                               >  2:4297762/1‑150 (MQ=255)
                                                                                                      ccGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGa                                          >  1:3700644/1‑150 (MQ=255)
                                                                                                       cGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGAc                                         <  1:2264229/150‑1 (MQ=255)
                                                                                                        gTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGAcg                                        <  1:4928533/150‑1 (MQ=255)
                                                                                                        gTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGAcg                                        <  2:3533026/150‑1 (MQ=255)
                                                                                                         tAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGAcgc                                       >  2:3097719/1‑150 (MQ=255)
                                                                                                              tGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATa                                  <  2:2803978/150‑1 (MQ=255)
                                                                                                               gCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAg                                 >  2:3495374/1‑150 (MQ=255)
                                                                                                                cAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGc                                >  1:101237/1‑150 (MQ=255)
                                                                                                                 aaaTTTCAGCACGGCAAATGCCAGACGGGTCCATTACTGGTTTTTAGAAGATACAGTCTTGTTTGCCTCAATAATGTATTCCCATAGACCGAAGCCAGTGAGTTGATAGTTAATTTGGGTTTCGCCTATAAAACCCAGACGCGCATAGct                               >  2:2932661/1‑150 (MQ=255)
                                                                                                                  aaTTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGctc                              >  2:301807/1‑150 (MQ=255)
                                                                                                                  aaTTTCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGctc                              >  2:5156154/1‑150 (MQ=255)
                                                                                                                     ttCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTgg                           <  2:4953836/150‑1 (MQ=255)
                                                                                                                     ttCAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTgg                           >  2:706472/1‑150 (MQ=255)
                                                                                                                       cAGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTc                         >  1:3815824/1‑150 (MQ=255)
                                                                                                                        aGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                        >  1:3366751/1‑150 (MQ=255)
                                                                                                                        aGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                        <  1:4755546/150‑1 (MQ=255)
                                                                                                                        aGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                        <  1:2885691/150‑1 (MQ=255)
                                                                                                                        aGCCCGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                        >  1:1911643/1‑150 (MQ=255)
                                                                                                                          cccGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGcc                      >  1:3478303/1‑150 (MQ=255)
                                                                                                                           ccGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCa                     >  1:774902/1‑150 (MQ=255)
                                                                                                                           ccGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCa                     <  2:3258261/150‑1 (MQ=255)
                                                                                                                           ccGGCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCa                     >  2:1386058/1‑150 (MQ=255)
                                                                                                                             ggCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATc                   <  1:1481889/150‑1 (MQ=255)
                                                                                                                             ggCAAATGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATc                   >  1:3199805/1‑150 (MQ=255)
                                                                                                                                aaaTGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCttt                >  2:1270495/1‑150 (MQ=255)
                                                                                                                                aaaTGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGACTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCttt                <  1:3928173/150‑1 (MQ=255)
                                                                                                                                 aaTGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTc               <  2:3770757/150‑1 (MQ=255)
                                                                                                                                 aaTGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTc               <  1:3684379/150‑1 (MQ=255)
                                                                                                                                  ctGCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCg              <  2:2681694/149‑1 (MQ=255)
                                                                                                                                    gCCAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCt            >  2:3502077/1‑150 (MQ=255)
                                                                                                                                     ccAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTg           >  1:4467451/1‑150 (MQ=255)
                                                                                                                                     ccAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTg           <  1:1232966/150‑1 (MQ=255)
                                                                                                                                      cAGACGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGa          >  2:2519513/1‑150 (MQ=255)
                                                                                                                                         aCGGTTCCATTACTGGTTTTTTGAAGATTCAGAGTTGTTTGCCTGAATATTTTATTCACATTTACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCAAAGCTCTGGTCGCCATCTTTCGCTGAAtt       >  1:129701/1‑150 (MQ=255)
                                                                                                                                          cGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTa      <  2:3611744/150‑1 (MQ=255)
                                                                                                                                          cGGGTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTa      >  2:1262625/1‑150 (MQ=255)
                                                                                                                                          cGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTa      >  1:2958704/1‑150 (MQ=255)
                                                                                                                                           gggTCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTAt     >  2:1441765/1‑150 (MQ=255)
                                                                                                                                            ggTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATc    >  1:264313/1‑150 (MQ=255)
                                                                                                                                              tCCATTACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAg  <  2:2794748/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTGTCTGCATTGGTATAAGAGTCACCGCCAAATTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGAACCGTAATCTGCAAATTTCAGCCCGGCAAATGCCAGACGGGTCCATGACTGGTTTTTTGAAGATTCAGTGTTGTTTGCCTGAATATTGTATTCCCATTGACCGTAGCCAGTGAGTTGATCGTTAATTTGGGTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAG  >  CP009273/1430726‑1431017

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: