Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,131,162 G→T 5.3% intergenic (+30/+82) PP_3634 → / ← PP_3635 GNAT family acetyltransferase/sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,131,1620GT 5.3% 460.1 / 6.0 149intergenic (+30/+82)PP_3634/PP_3635GNAT family acetyltransferase/sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (62/77);  new base T (2/7);  total (65/84)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 153 total aligned reads displayed.

TGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACG  >  NC_002947/4131014‑4131307
                                                                                                                                                    |                                                                                                                                                 
tGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGC                                                                                                                                                  >  2:203579‑M1/1‑150 (MQ=255)
      aCCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAG                                                                                                                                            <  2:243927‑M1/150‑1 (MQ=255)
       ccAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGG                                                                                                                                           >  2:1415451‑M1/1‑150 (MQ=255)
        cagcagGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGt                                                                                                                                          >  1:1647288‑M1/1‑149 (MQ=255)
        cagcagGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGa                                                                                                                                          >  1:1278358/1‑150 (MQ=255)
                tGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTa                                                                                                                                  <  2:372982/150‑1 (MQ=255)
                tGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTCGCCCGCTCCCACAGGtacatcacc                                                                                                                                  >  1:2391201‑M1/1‑141 (MQ=255)
                     aTGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcgaa                                                                                                                                                 <  2:953284/130‑3 (MQ=255)
                     aTGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcgaa                                                                                                                                                 >  1:953284/1‑128 (MQ=255)
                           tGATCGAAGAGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCGCAGGAATGTGTGCGCCAGCGCCCGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGc                                                                                                                       <  1:2348755/150‑1 (MQ=255)
                            gATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGcc                                                                                                                      >  2:1567534/1‑150 (MQ=255)
                              tCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCgg                                                                                                                    <  1:1518322/150‑1 (MQ=255)
                              tCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCgg                                                                                                                    >  2:2292765/1‑150 (MQ=255)
                              tCGAAGCGGAGCAGATGCACGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCCCAGGAATACCGTAAGGCCCAAGGCCgg                                                                                                                    <  1:1777846/150‑1 (MQ=255)
                               cggaGCGGAGCAGTGGGGCGATTGGCGGCATGATGAGCCGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCTGATGGCCTGTTCGCGGGCTGGCGCGGGCCCATAGGAATACCGTAAGGCCCACGGCCGGt                                                                                                                   <  2:784504/147‑1 (MQ=255)
                               cGAAGCGGAGCAGATGCGCGATGAGCGGCTGGCTGAGCCCTGTTCGTAATGCCAGGAATGTGTGCGCCAGAGTTCAGACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                                   <  1:1739081/150‑1 (MQ=255)
                                 aaGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCCTAGTAATACCGTAAGGCCCAAGGCCGGTGa                                                                                                                 <  2:760629/150‑1 (MQ=255)
                                  aGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCTCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                                <  1:1789604/150‑1 (MQ=255)
                                  aGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                                <  1:1524920/150‑1 (MQ=255)
                                  aGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                                <  1:2179627/150‑1 (MQ=255)
                                  aGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                                <  2:1243293/150‑1 (MQ=255)
                                           gATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGtacatcaccggccttgggccttacggtattcctgtg                                                                                                       <  2:1917358‑M1/150‑37 (MQ=255)
                                           gATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgtg                                                                                                       >  2:1371187‑M2/1‑144 (MQ=255)
                                           gATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTCGCCCGCTCCCCCAGGtacatcaccggccttgggccttacggtattcctgtg                                                                                                       <  2:2271148‑M1/150‑37 (MQ=255)
                                            aTGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCNCNCGCCCANAGGAATACCGTANGGCCCAAGGCCGGTGATGTAcctgngg                                                                                                      >  1:1576422‑M2/1‑143 (MQ=255)
                                              gcgcGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGa                                                                                                    <  1:1882425/150‑1 (MQ=255)
                                                 cgATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCg                                                                                                 <  2:185610/150‑1 (MQ=255)
                                                  gATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAAGAATACCGTAAGGCCCAAGGCCGGAGATGTACCTGTGGTAGCgg                                                                                                >  1:2393026/1‑150 (MQ=255)
                                                      aGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGGAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGa                                                                                            <  2:2052229/150‑1 (MQ=255)
                                                      aGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGa                                                                                            <  1:985938/150‑1 (MQ=255)
                                                       gCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAg                                                                                           >  2:1266145/1‑150 (MQ=255)
                                                       gCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgcgggcaag                                                                                           <  1:704180‑M1/150‑49 (MQ=255)
                                                              gATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                                <  2:1118432/122‑1 (MQ=255)
                                                              gATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                                >  1:1118432/1‑122 (MQ=255)
                                                               aTGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTa                                                                                                                                  >  1:2237481/1‑103 (MQ=255)
                                                               aTGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTa                                                                                                                                  <  2:2237481/103‑1 (MQ=255)
                                                                 gAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgcgggcaagcccgcgaaca                                                                                 <  2:1100783‑M1/150‑59 (MQ=255)
                                                                  aGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCgaagaa                                                                                <  2:776613/150‑1 (MQ=255)
                                                                   gCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCggg                                                                                               <  1:2088784/134‑1 (MQ=255)
                                                                   gCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGaagaag                                                                               >  2:1302302/1‑150 (MQ=255)
                                                                      gACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGccc                                                                            >  1:281129/1‑150 (MQ=255)
                                                                          gCTCAAGCCCATGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCNGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                                <  1:724102/110‑1 (MQ=255)
                                                                          gCTCAAGCCCATGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                                >  2:724102/1‑110 (MQ=255)
                                                                          gCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCAtt                                                                        <  2:403059/150‑1 (MQ=255)
                                                                            tCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcc                                                                             >  1:854090/1‑143 (MQ=255)
                                                                            tCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcc                                                                             <  2:854090/143‑1 (MQ=255)
                                                                            tCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTgg                                                                      >  1:2409849/1‑150 (MQ=255)
                                                                                 cccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCTTCTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgc                                                                                                  >  2:661482‑M1/1‑76 (MQ=255)
                                                                                 cccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCTTCTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgc                                                                                                  <  1:661482‑M1/117‑42 (MQ=255)
                                                                                  ccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCTGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTc                                                                >  2:2237988/1‑150 (MQ=255)
                                                                                   cAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                               <  2:505482/150‑1 (MQ=255)
                                                                                    aGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAg                                                                                                   <  2:2391281/113‑1 (MQ=255)
                                                                                    aGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAg                                                              <  1:1266145/150‑1 (MQ=255)
                                                                                    aGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATAATATGTGTTACAAGTCCTGGGATGTACCTGTGGGAGCGGGCGACACCCTGGAGGATGGCCATTTGCCTTTCAg                                                              >  1:44490/1‑150 (MQ=255)
                                                                                               cgcCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGNGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCa                                                   <  1:270279/150‑1 (MQ=255)
                                                                                                gcCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCTTCCGAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAg                                                  >  1:738530/1‑150 (MQ=255)
                                                                                                   aGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCTTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGCGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAgcgc                                               <  2:1278358/150‑1 (MQ=255)
                                                                                                      gcAAGAACTTGCATTGATCCGAGGCCCTGTTCGCAGCCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGa                                            <  2:1083728/150‑1 (MQ=255)
                                                                                                        aataaCGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATg                                          >  1:1667724/4‑150 (MQ=255)
                                                                                                        aaGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATg                                          >  2:1578360/1‑150 (MQ=255)
                                                                                                          gAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                               >  1:2028790/1‑127 (MQ=255)
                                                                                                          gAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCa                                        <  1:370054/150‑1 (MQ=255)
                                                                                                              gTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgtgggagcgggcgagcacccggagaagccatcggatcaatgcacgttcttgcgctggcgcacacattcct                                    <  1:760055‑M2/150‑74 (MQ=255)
                                                                                                              gTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCAtt                                                                        >  1:1225530/1‑114 (MQ=255)
                                                                                                              gTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCAtt                                                                        <  2:1225530/114‑1 (MQ=255)
                                                                                                               tGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCCAGGCCGGTGATGTAcctgtgggagcgggcgagcccgcgaagaagccatcggatcaatgcacgttcttgcgctggcgcacacattcctg                                   <  1:1519054‑M2/150‑75 (MQ=255)
                                                                                                                     gATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCgg                                                                                                <  2:560803/83‑1 (MQ=255)
                                                                                                                        ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                          <  2:370413/150‑1 (MQ=255)
                                                                                                                        ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                          <  1:282183/150‑1 (MQ=255)
                                                                                                                        ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                          <  1:1415451/150‑1 (MQ=255)
                                                                                                                        ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                          <  1:1567534/150‑1 (MQ=255)
                                                                                                                         cGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCagag                         <  1:1302302/150‑1 (MQ=255)
                                                                                                                         cGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCagag                         >  1:2116946/1‑150 (MQ=255)
                                                                                                                         cGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCagag                         >  2:2028306/1‑150 (MQ=255)
                                                                                                                           aTGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCTCCCAGCGCCGATGCAGCCCACgcgc                              >  1:89906/1‑143 (MQ=255)
                                                                                                                           aTGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGTCCACGCGCCAGAGCa                       <  2:1415035/150‑1 (MQ=255)
                                                                                                                           aTGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACgcgc                              <  2:89906/143‑1 (MQ=255)
                                                                                                                            tGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                                   <  1:1421920/57‑1 (MQ=255)
                                                                                                                            tGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCAt                      >  2:1159962/1‑150 (MQ=255)
                                                                                                                              gCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATcc                    <  1:820810/150‑1 (MQ=255)
                                                                                                                               ccTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCa                   >  1:1343479/1‑150 (MQ=255)
                                                                                                                                cTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgt                                                                                                        <  1:1839347‑M2/64‑6 (MQ=255)
                                                                                                                                cTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgt                                                                                                        >  2:1839347‑M2/1‑59 (MQ=255)
                                                                                                                                cTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAg                  <  2:281129/150‑1 (MQ=255)
                                                                                                                                 tGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGc                 >  2:191475/1‑150 (MQ=255)
                                                                                                                                 tGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGc                 >  2:178991/1‑150 (MQ=255)
                                                                                                                                    gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGcc                                     <  2:1767060/126‑1 (MQ=255)
                                                                                                                                    gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGca                >  1:2014707/2‑148 (MQ=255)
                                                                                                                                    gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGca                <  2:2014707/147‑1 (MQ=255)
                                                                                                                                      gcgGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCAGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCACCCCCGATGCAGCCCACGCGCCAGACCACCCACCACgaa            >  2:444228/1‑150 (MQ=255)
                                                                                                                                        gggCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCttt                                                                 >  1:653044/1‑95 (MQ=255)
                                                                                                                                           cTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgtgggagcgggcgagcccgcgaacaggccatcggatcaatgcacgttcttgcgctggcgcacacattcctgggcttgagcggtcagctcatccagccgc       <  1:704180‑M2/150‑103 (MQ=255)
                                                                                                                                           cTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCa       >  1:77364/1‑150 (MQ=255)
                                                                                                                                           cTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTCGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGACGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCa       <  1:2292765/150‑1 (MQ=255)
                                                                                                                                                ccGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgtgggcgcgggcaagcccgcgaacaggccatcggatcaatgcacgttcttgcgctggcgcacacattcctgggcttgagcggtcagctcatccagccgctcatc  <  1:2167945‑M2/150‑108 (MQ=255)
                                                                                                                                                ccGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGaagaag                                                                               <  2:1732821/73‑1 (MQ=255)
                                                                                                                                                ccGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCAcgc                                <  1:1514766/120‑1 (MQ=255)
                                                                                                                                                ccGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCAcgc                                >  2:1514667/1‑120 (MQ=255)
                                                                                                                                                ccGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCAcgc                                >  2:1514766/1‑120 (MQ=255)
                                                                                                                                                  gcgcCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgtgggagcgggcgagcccgcgaagaagccatcggatcaatgcacgttcttgcgctggcgcacacattcctggg                                 >  2:661482‑M2/1‑41 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACG  >  NC_002947/4131014‑4131307

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