Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,600,990 C→T 5.7% intergenic (‑207/+38) 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,9900CT 5.7% 435.6 / 9.5 142intergenic (‑207/+38)rpoH/ftsXRNA polymerase, sigma 32 (sigma H) factor/putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (58/73);  new base T (3/5);  total (63/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 142 total aligned reads displayed.

GCTTAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC  >  NC_000913/3600842‑3601138
                                                                                                                                                    |                                                                                                                                                    
gCTTAAGAGCAACTAGGTTATCATTCACTGTTTTATCAGACCGTGATATTATCCACAAGTTCAATTCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACg                                                                                                                                                     <  1:111305/150‑1 (MQ=255)
 ctaaaGAGCAACGAGGTTAGCATTCACTGTTTTATCAGACCGGGATTTTATCCACAACTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGATACAACGGGGTCTCTTTCCCTGCTACgg                                                                                                                                                    <  2:570939/147‑1 (MQ=255)
  ttAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATATAGATGAAAAACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGa                                                                                                                                                   <  2:1608437/150‑1 (MQ=255)
   tAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGTGATGAGACAACGGGGTATCTTTAACTACTAAGGaa                                                                                                                                                  >  1:85620/1‑150 (MQ=255)
    aaGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATTAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAc                                                                                                                                                 <  2:71756/150‑1 (MQ=255)
    aaGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAc                                                                                                                                                 >  2:405199/1‑150 (MQ=255)
     agagCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATAACGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAAcc                                                                                                                                                >  2:5723/1‑150 (MQ=255)
         cAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCAtt                                                                                                                                            >  2:1330468/1‑150 (MQ=255)
          aaCGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTg                                                                                                                                           >  1:1573413/1‑150 (MQ=255)
           aCGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTAAGGAACCCATTGc                                                                                                                                          >  1:1573018/1‑150 (MQ=255)
             gAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAg                                                                                                                                        <  2:424474/150‑1 (MQ=255)
               ggTTATCATTCACTGTTTTCTCAGACCGTGATTTTATCCACAAGGTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAggg                                                                                                                                      <  2:131011/150‑1 (MQ=255)
               ggTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAggg                                                                                                                                      <  1:1705492/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)
                     cATTCACTGTTTTATCAGACCGAGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCAATGGGTTCCGTAGCAGGGAAagag                                                                                                                                <  2:60791/150‑1 (MQ=38)
                       ttCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGta                                                                                                                              >  2:631947/1‑150 (MQ=255)
                        tCACTGTTTTATCAGACCTTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCTTTGCAGGGAAAGAGtgt                                                                                                                             <  2:772481/150‑3 (MQ=255)
                        tCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtgt                                                                                                                             <  2:616920/150‑3 (MQ=255)
                        tCACTGTTTTACCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtaa                                                                                                                             >  2:367145/1‑149 (MQ=255)
                         cACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGtata                                                                                                                            >  2:1569413/1‑150 (MQ=255)
                           cTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAac                                                                                                                          <  2:1201703/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)
                               tttATCAGACCGTGATTTTATCAACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGc                                                                                                                      <  2:1154742/150‑1 (MQ=255)
                                  aTCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCttt                                                                                                                   <  1:374829/150‑1 (MQ=255)
                                     aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                <  2:1666609/150‑1 (MQ=255)
                                     aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                <  2:1040451/150‑1 (MQ=255)
                                     aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                <  1:1539155/150‑1 (MQ=255)
                                     aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCACTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                >  1:1476542/1‑150 (MQ=255)
                                     aGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttat                                                                                                                <  1:599051/150‑1 (MQ=255)
                                       aCCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTttatta                                                                                                              >  2:1139772/1‑150 (MQ=255)
                                         cGTGATTTTATCCACAAGTTGAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTTCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                            >  1:1315119/1‑150 (MQ=255)
                                         cGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTattatt                                                                                                            >  1:559798/1‑150 (MQ=255)
                                           tGATTTTATCCACAAGTTCAATGCAATCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCa                                                                                                          <  2:1021998/150‑1 (MQ=255)
                                             aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                        <  1:804665/150‑1 (MQ=255)
                                             aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                        <  2:256909/150‑1 (MQ=255)
                                             aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                        <  1:410877/150‑1 (MQ=255)
                                             aTTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAgg                                                                                                        <  1:299300/150‑1 (MQ=255)
                                                ttATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGt                                                                                                     >  2:1284341/1‑150 (MQ=255)
                                                 tataCACAATTTCAATGAAATCTGTTGAATAAATTACGGACAATATGTGACATAGAGATGTAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAATCCATTTCAGGGAATGAGTATAACAGGCTTTTATTATTCAGGCGTa                                                                                                    <  1:1416576/150‑1 (MQ=255)
                                                   tCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTaaa                                                                                                  >  1:1063133/1‑150 (MQ=255)
                                                     cacaAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTACCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGt                                                                                                >  1:186444/1‑150 (MQ=255)
                                                         aGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAAAACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCg                                                                                            >  1:723745/1‑150 (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)
                                                                    aGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTa                                                                                 >  2:1548267/1‑150 (MQ=255)
                                                                        tgtgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGt                                                                             >  2:300448/1‑150 (MQ=255)
                                                                          tgAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCATTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTgg                                                                           >  1:1226355/1‑150 (MQ=255)
                                                                           gAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGc                                                                          <  1:1411852/150‑1 (MQ=255)
                                                                             aTAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCaa                                                                        >  1:1178083/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:1501035/1‑150 (MQ=255)
                                                                                      gCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTg                                                               <  1:1569413/150‑1 (MQ=255)
                                                                                       cacaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTAAGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGc                                                              <  2:1573018/150‑1 (MQ=255)
                                                                                       cacaAAATGTGACATAGAGAGGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGCCCCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTCCCGTGGCAAGCCCCGCTGc                                                              >  1:206794/1‑150 (MQ=255)
                                                                                        acaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGcc                                                             <  2:556039/150‑1 (MQ=255)
                                                                                        acaAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGcc                                                             <  2:821109/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)
                                                                                                        agaTGAAATACCGGGAAGAGACAACGGTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCg                                             <  1:405199/150‑1 (MQ=255)
                                                                                                        agaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCg                                             <  2:196148/150‑1 (MQ=255)
                                                                                                         gaTGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGa                                            >  2:860455/1‑150 (MQ=255)
                                                                                                          aTGAAATACCGGGAAGAGACAACGGGGGCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAg                                           <  1:1076935/150‑1 (MQ=255)
                                                                                                               aTACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATAc                                      >  1:1032317/1‑150 (MQ=255)
                                                                                                                 aCCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACca                                    >  1:717698/1‑150 (MQ=255)
                                                                                                                   gggaaaGAGACCCCGTTGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCCGCCAATCATCGAGCATACcagc                                  >  2:518504/5‑150 (MQ=255)
                                                                                                                   cGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACcagc                                  >  2:40773/1‑150 (MQ=255)
                                                                                                                           gaCAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCCCCCAGCCAATCATCGAGCATACCAGCAGCAATAg                          >  2:147152/1‑150 (MQ=255)
                                                                                                                             cAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCa                        >  2:624131/1‑150 (MQ=255)
                                                                                                                             cAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGGTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCa                        <  2:501013/150‑1 (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)
                                                                                                                                      ctctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACAGTGGCAAGCCCCGCTGCCACCCAGCCAATAATCGAGCATACCAGCAGCAATAGCAGGCATTCAt               >  2:116864/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)
                                                                                                                                        ctTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCCAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCg             >  1:1375983/1‑150 (MQ=255)
                                                                                                                                         tttCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGTGTAAATTGGCGTAAATGTTGTACCGTGGCAAGAAACGATGCAACCCAGCCAATCATCGAGCATACCAGAAGCAATAGCAGGCATTCATCGa            >  2:342164/1‑150 (MQ=255)
                                                                                                                                           tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt          <  2:1063133/150‑1 (MQ=255)
                                                                                                                                           tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt          <  2:723745/150‑1 (MQ=255)
                                                                                                                                           tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt          <  2:225575/150‑1 (MQ=255)
                                                                                                                                           tCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAAt          <  2:669204/150‑1 (MQ=255)
                                                                                                                                              cTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGAt       >  2:1287089/1‑150 (MQ=255)
                                                                                                                                                   aCGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGcc  <  1:367145/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTTAAGAGCAACGAGGTTATCATTCACTGTTTTATCAGACCGTGATTTTATCCACAAGTTCAATGCAAGCTTGTGAATAAATTACGCACAAAATGTGACATAGAGATGAAATACCGGGAAGAGACAACGGGGTCTCTTTCCCTGCTACGGAACCCATTGCAGGGAAAGAGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC  >  NC_000913/3600842‑3601138

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