Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,600,994 A→T 5.4% intergenic (‑211/+34) rpoH ← / ← ftsX RNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,600,9940AT 5.4% 401.7 / 6.8 138intergenic (‑211/+34)rpoH/ftsXRNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (59/70);  new base T (3/5);  total (63/75)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

TAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCAT  >  NC_000913/3600845‑3601140
                                                                                                                                                     |                                                                                                                                                  
 aaGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATTAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAc                                                                                                                                                   <  2:71756/150‑1 (MQ=255)
 aaGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAc                                                                                                                                                   >  2:405199/1‑150 (MQ=255)
       aaCGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTg                                                                                                                                             >  2:1125703/1‑150 (MQ=255)
       aaCGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTg                                                                                                                                             >  1:1573413/1‑150 (MQ=255)
          gAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAg                                                                                                                                          <  2:424474/150‑1 (MQ=255)
          gAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAg                                                                                                                                          <  2:1571651/150‑1 (MQ=255)
            ggTTATCATTCACTGTTTTCTCAGACCGTGATTTTATCCACAAGGTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAggg                                                                                                                                        <  2:131011/150‑1 (MQ=255)
             gTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGa                                                                                                                                       <  2:1477682/150‑1 (MQ=255)
                aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGATGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAag                                                                                                                                    <  2:543237/150‑1 (MQ=255)
                aTCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAag                                                                                                                                    <  2:1650562/150‑1 (MQ=255)
                  gaTTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAATATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCCCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAagag                                                                                                                                  <  1:514540/149‑1 (MQ=255)
                  cATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAagag                                                                                                                                  >  2:1236988/1‑150 (MQ=255)
                  cATTCACTGTTTTATCAGACCGAGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAagag                                                                                                                                  <  2:60791/150‑1 (MQ=38)
                     tCACTGTTTTACCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtaa                                                                                                                               >  2:367145/1‑149 (MQ=255)
                       aCTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATaa                                                                                                                             >  1:843912/1‑150 (MQ=255)
                        cTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAac                                                                                                                            <  2:1201703/150‑1 (MQ=255)
                        cTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGAAATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAac                                                                                                                            <  2:621343/150‑1 (MQ=255)
                          gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAGATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                          >  1:980532/1‑150 (MQ=255)
                          gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                          >  1:196148/1‑150 (MQ=255)
                          gTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAacac                                                                                                                          >  2:1424605/1‑150 (MQ=255)
                           ttttATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACg                                                                                                                         <  2:1412967/150‑1 (MQ=255)
                            tttATCAGACCGTGATTTTATCAACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGc                                                                                                                        <  2:1154742/150‑1 (MQ=255)
                               aTCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCttt                                                                                                                     <  1:374829/150‑1 (MQ=255)
                                tCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCtttt                                                                                                                    <  2:53427/150‑1 (MQ=255)
                                  aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                  <  2:1666609/150‑1 (MQ=255)
                                  aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                  <  2:1040451/150‑1 (MQ=255)
                                  aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCACTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                  >  1:1476542/1‑150 (MQ=255)
                                  aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                  <  1:599051/150‑1 (MQ=255)
                                   gACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttatt                                                                                                                 <  2:604351/150‑1 (MQ=255)
                                      cGTGATTTTATCCACAAGTTGAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTTCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                              >  1:1315119/1‑150 (MQ=255)
                                      cGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                              >  1:1208642/1‑150 (MQ=255)
                                      cGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                              >  1:559798/1‑150 (MQ=255)
                                          aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                          <  1:804665/150‑1 (MQ=255)
                                          aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                          <  1:410877/150‑1 (MQ=255)
                                          aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                          <  1:299300/150‑1 (MQ=255)
                                          aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                          <  2:256909/150‑1 (MQ=255)
                                             ttATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGt                                                                                                       >  2:1284341/1‑150 (MQ=255)
                                              tATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTa                                                                                                      <  2:791522/150‑1 (MQ=255)
                                                tCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTaaa                                                                                                    >  1:1063133/1‑150 (MQ=255)
                                                  cacaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGt                                                                                                  >  1:1378500/1‑150 (MQ=255)
                                                  cacaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTACCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGt                                                                                                  >  1:186444/1‑150 (MQ=255)
                                                    caAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTgg                                                                                                <  2:1327856/150‑1 (MQ=255)
                                                      aGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAAAACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCg                                                                                              >  1:723745/1‑150 (MQ=255)
                                                      aGTTAAATGTAATCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTTCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGTCg                                                                                              <  2:1164498/150‑1 (MQ=255)
                                                         tCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTaa                                                                                           <  1:1724358/150‑1 (MQ=255)
                                                           aaTGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAt                                                                                         >  1:7908/1‑150 (MQ=255)
                                                            aTGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtg                                                                                        <  1:185702/150‑1 (MQ=255)
                                                              gCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgtt                                                                                      <  2:1573413/150‑1 (MQ=255)
                                                               cAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttg                                                                                     <  1:1160368/150‑1 (MQ=255)
                                                                aaGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAAtgttgt                                                                                    <  2:660654/150‑1 (MQ=255)
                                                                 aGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTa                                                                                   >  2:1548267/1‑150 (MQ=255)
                                                                     tgtgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGt                                                                               >  2:300448/1‑150 (MQ=255)
                                                                        gAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGc                                                                            <  1:1411852/150‑1 (MQ=255)
                                                                         aaTAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAGGTGGCGTAAATGTTGTACCGTGGCa                                                                           <  2:568625/150‑1 (MQ=255)
                                                                          aTAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGGTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGTCaa                                                                          >  2:435033/1‑150 (MQ=255)
                                                                              aTTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCa                                                                      >  1:204197/1‑150 (MQ=255)
                                                                                 aCGCACAAAATGTTACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                   <  2:881468/150‑1 (MQ=255)
                                                                                 aCGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                   <  2:843912/150‑1 (MQ=255)
                                                                                 aCGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGc                                                                   <  2:1713966/150‑1 (MQ=255)
                                                                                  cGCACAAAATGTGACATAGATATGAAATACCGGGAAGAGACAACGGGGTCTATTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCTTGCCACggt                                                                  <  2:699902/150‑3 (MQ=255)
                                                                                   gCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTg                                                                 <  1:1569413/150‑1 (MQ=255)
                                                                                   gCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTg                                                                 >  1:1501035/1‑150 (MQ=255)
                                                                                    cacaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTAAGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGc                                                                <  2:1573018/150‑1 (MQ=255)
                                                                                     acaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGcc                                                               <  2:821109/150‑1 (MQ=255)
                                                                                       aaaaTGTGACATAGAGATGAAATACCTGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCAc                                                             <  2:1678095/150‑1 (MQ=255)
                                                                                       aaaaTGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCAc                                                             <  2:1723093/150‑1 (MQ=255)
                                                                                       aaaaTGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCAc                                                             <  2:221029/150‑1 (MQ=255)
                                                                                             tgACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAGCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGcc                                                       <  1:288728/150‑1 (MQ=255)
                                                                                              gACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCa                                                      >  1:41050/1‑150 (MQ=255)
                                                                                                  tAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAatca                                                  >  2:605046/1‑150 (MQ=255)
                                                                                                    gagaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAatcatc                                                <  1:1260605/150‑1 (MQ=255)
                                                                                                     agaTGAAATACCGGGAAGAGACAACGGTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCg                                               <  1:405199/150‑1 (MQ=255)
                                                                                                     agaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCg                                               <  2:196148/150‑1 (MQ=255)
                                                                                                      gaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGa                                              <  1:1548267/150‑1 (MQ=255)
                                                                                                      gaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGa                                              >  2:860455/1‑150 (MQ=255)
                                                                                                       aTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAg                                             <  2:300650/150‑1 (MQ=255)
                                                                                                       aTGAAATACCGGGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAg                                             <  1:1076935/150‑1 (MQ=255)
                                                                                                            aTACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATAc                                        >  1:1032317/1‑150 (MQ=255)
                                                                                                              aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACca                                      >  1:717698/1‑150 (MQ=255)
                                                                                                               ccGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCGTCGAGCATACcag                                     >  1:1651405/1‑150 (MQ=255)
                                                                                                               ccGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcag                                     <  1:989549/150‑1 (MQ=255)
                                                                                                                gggaaaGAGACCCCGTTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCCGCCAATCATCGAGCATACcagc                                    >  2:518504/5‑150 (MQ=255)
                                                                                                                cGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcagc                                    >  2:40773/1‑150 (MQ=255)
                                                                                                                 ggaaaGAGACCCCGTTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcagca                                   >  1:211575/4‑150 (MQ=255)
                                                                                                                          cAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCa                          >  2:624131/1‑150 (MQ=255)
                                                                                                                            aCGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAgg                        <  2:559798/150‑1 (MQ=255)
                                                                                                                             cGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGc                       <  1:1236988/150‑1 (MQ=255)
                                                                                                                               gggTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCAt                     >  2:1348841/1‑150 (MQ=255)
                                                                                                                                  actctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCa                  >  2:211589/2‑150 (MQ=38)
                                                                                                                                   ctctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCAt                 >  1:1654584/1‑150 (MQ=255)
                                                                                                                                   ctctTTCCCTGCAATGGGTTCCGTAGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCAt                 >  1:1116083/1‑150 (MQ=39)
                                                                                                                                    tctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACAAGCTTTTATTATTCAGGCGTAAAGTGTCGTAAATGTTTTACCGTGGCAAGCAACGCTGCCACACAACAAATCATAGAGCAAACCAGCAGAAAAAGCAGGCATTCATc                >  2:262168/1‑150 (MQ=255)
                                                                                                                                     ctTTCCCTGCTACGGAACCCATTGCAGTGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCg               <  1:118445/150‑1 (MQ=255)
                                                                                                                                        tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt            <  2:1063133/150‑1 (MQ=255)
                                                                                                                                        tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt            <  2:225575/150‑1 (MQ=255)
                                                                                                                                        tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt            <  2:723745/150‑1 (MQ=255)
                                                                                                                                        tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt            <  1:1476805/150‑1 (MQ=255)
                                                                                                                                        tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt            <  2:669204/150‑1 (MQ=255)
                                                                                                                                            tGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATa        <  1:1139772/150‑1 (MQ=255)
                                                                                                                                                aCGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGcc    <  1:367145/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCAT  >  NC_000913/3600845‑3601140

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