Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,512,561 Δ1 bp intergenic (‑217/+73) yhfZ ← / ← trpS putative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,512,5600T.97.9% 392.8 / 1.3 93intergenic (‑216/+74)yhfZ/trpSputative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (52/39);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 102 total aligned reads displayed.

ATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCA  >  NC_000913/3512461‑3512655
                                                                                                   |                                                                                               
aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt                                                                                                 <  2:824617/100‑1 (MQ=255)
aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt                                                                                                 >  2:34758/1‑100 (MQ=255)
  aaaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtac                                                                                               >  2:2394635/1‑98 (MQ=255)
  aaaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtac                                                                                               >  2:1105124/1‑98 (MQ=255)
   aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtact                                                                                              >  2:2372571/1‑97 (MQ=255)
   aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtact                                                                                              >  1:2948031/1‑97 (MQ=255)
    aaaaaCTGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtactt                                                                                             <  2:1973933/100‑5 (MQ=255)
    aaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtactt                                                                                             >  2:1080095/1‑96 (MQ=255)
    aaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtactt                                                                                             >  2:320205/1‑96 (MQ=255)
      aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGGTACGCTTATCAGACCTACATTTCCTTTCAATATGTGCATACTTTg                                                                                          >  2:2393295/1‑100 (MQ=37)
      aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTg                                                                                          >  1:1455025/1‑100 (MQ=255)
      aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTg                                                                                          >  2:289354/1‑100 (MQ=255)
      aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTg                                                                                          <  1:2584726/100‑1 (MQ=255)
       aaCCGGGTTTCCCCGGTTTCAGAGTGATGATTAAAGCAATATTGTCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGt                                                                                         <  1:2531869/100‑1 (MQ=37)
       aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGt                                                                                         <  2:1274960/100‑1 (MQ=255)
       aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGt                                                                                         <  2:2844609/100‑1 (MQ=255)
       aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGt                                                                                         <  2:2011007/100‑1 (MQ=255)
        aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGt                                                                                         >  1:2785276/1‑100 (MQ=255)
         ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATCCTTTGTAg                                                                                       >  1:886141/1‑100 (MQ=255)
         ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAg                                                                                       >  1:1178210/1‑100 (MQ=255)
         ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAg                                                                                       <  1:2207612/100‑1 (MQ=255)
          cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAgg                                                                                      <  1:2433549/100‑1 (MQ=255)
                  cccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATaa                                                                              >  1:1220735/1‑100 (MQ=255)
                   ccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATaa                                                                              <  2:2675606/99‑1 (MQ=255)
                   ccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAg                                                                             <  1:381360/100‑1 (MQ=255)
                    cGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAgg                                                                            <  1:1534020/100‑1 (MQ=255)
                     ggTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGc                                                                           <  1:2787110/100‑1 (MQ=255)
                     ggTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGc                                                                           <  2:1080905/100‑1 (MQ=255)
                       tttCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGt                                                                         <  1:303832/100‑1 (MQ=255)
                       tttCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGt                                                                         <  2:314171/100‑1 (MQ=255)
                       tttCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGt                                                                         >  2:1002848/1‑100 (MQ=255)
                         tCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTc                                                                       >  2:1252455/1‑100 (MQ=255)
                          cAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCa                                                                      >  1:1715448/1‑100 (MQ=255)
                          cAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCa                                                                      >  2:1157621/1‑100 (MQ=255)
                              gTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                  >  2:1680414/1‑100 (MQ=255)
                                gatgatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                >  1:779017/1‑100 (MQ=255)
                                     tAAAAGCGAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                           >  1:1595362/1‑100 (MQ=255)
                                     tAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                           >  1:2241332/1‑100 (MQ=255)
                                      aaaaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                          >  2:791839/1‑100 (MQ=255)
                                       aaaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                         <  2:489705/100‑1 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGGTAAGGCGTTCACGCCGCATCCGGCa                                                        >  1:1343769/1‑100 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGACCTACATTTCCTTGCAATATGTGCATACTTTGTAGGACGGATAAGGCGTTCACGCCGCATCCGGc                                                         >  2:666461/1‑99 (MQ=37)
                                        aaGCAAAATTGCCTGATGCGATACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                        >  2:1850025/1‑100 (MQ=37)
                                         aGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                       >  1:2210509/1‑100 (MQ=255)
                                         aGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                       <  2:1133055/100‑1 (MQ=255)
                                         aGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGCTCACGCCGCATCCGGCAt                                                       >  1:1996684/1‑100 (MQ=255)
                                          gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                        <  1:1806308/98‑2 (MQ=21)
                                           cAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                     >  1:271985/1‑100 (MQ=255)
                                            aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                    >  1:1655256/1‑100 (MQ=255)
                                             aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAc                                                   >  1:826742/1‑100 (MQ=255)
                                              aaTTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACa                                                  >  1:1324919/1‑100 (MQ=255)
                                                  gCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAgc                                              >  2:300176/1‑100 (MQ=255)
                                                  gCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAgc                                              <  1:2664928/100‑1 (MQ=255)
                                                    cTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAgcgc                                            >  1:1519583/1‑100 (MQ=255)
                                                     tGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCa                                           >  1:72068/1‑100 (MQ=255)
                                                     tGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCa                                           >  1:2660693/1‑100 (MQ=255)
                                                       aTGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCTGAACACAGCGCAAt                                         <  2:2840355/100‑1 (MQ=25)
                                                       aTGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAAt                                         <  1:441903/100‑1 (MQ=255)
                                                         gcgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAAttt                                       <  2:718028/100‑1 (MQ=255)
                                                          cgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTg                                      >  1:2522506/1‑100 (MQ=255)
                                                          cgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTg                                      >  1:2198228/1‑100 (MQ=255)
                                                           gcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGc                                     >  2:2991710/1‑100 (MQ=255)
                                                            cTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGcc                                    >  1:1691033/1‑100 (MQ=255)
                                                            cTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGcc                                    <  2:52592/100‑1 (MQ=255)
                                                              aCGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                  >  1:52593/1‑100 (MQ=255)
                                                              aCGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                  <  2:1740470/100‑1 (MQ=255)
                                                              aCGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                  <  1:2172837/100‑1 (MQ=255)
                                                                gCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCa                                >  1:1658548/1‑100 (MQ=255)
                                                                   tATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATa                             >  1:831502/1‑100 (MQ=255)
                                                                    aTCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAg                            >  1:2220164/1‑100 (MQ=255)
                                                                    aTCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAg                            <  1:251730/100‑1 (MQ=255)
                                                                       aGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGa                         <  2:1467707/100‑1 (MQ=255)
                                                                        ggCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGaa                        >  1:1409557/1‑100 (MQ=255)
                                                                        ggCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGaa                        >  1:3029721/1‑100 (MQ=255)
                                                                        ggCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGaa                        >  2:2274595/1‑100 (MQ=255)
                                                                        ggCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGaa                        >  1:1284735/1‑100 (MQ=255)
                                                                         gCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAt                       >  2:47948/1‑100 (MQ=255)
                                                                          cttACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAtt                      >  1:2059648/3‑100 (MQ=255)
                                                                            tACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTAc                    >  1:106958/1‑100 (MQ=255)
                                                                              cATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACgg                  <  2:1416215/100‑1 (MQ=255)
                                                                                 ttCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCt                <  1:47464/99‑1 (MQ=255)
                                                                                   ccTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCg             <  1:1090083/100‑1 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc            <  1:1320356/100‑1 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc            >  2:1371324/1‑100 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc            >  1:2029657/1‑100 (MQ=255)
                                                                                       gCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcac         <  1:2504212/100‑1 (MQ=255)
                                                                                        cAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcaca        <  1:1002858/100‑1 (MQ=255)
                                                                                        cAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcaca        <  2:380470/100‑1 (MQ=255)
                                                                                        cAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcaca        >  2:2993448/1‑100 (MQ=255)
                                                                                        cAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcac         <  1:2012890/99‑1 (MQ=255)
                                                                                         aaTATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACaa       >  2:2196726/1‑100 (MQ=255)
                                                                                          atatGTGCATACTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCTATAGTGAATTACGGCTTCGCCACaaa      <  1:1850041/100‑1 (MQ=21)
                                                                                            atGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAc    <  2:2897220/100‑1 (MQ=255)
                                                                                             tgtgCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAc    <  2:2785251/100‑1 (MQ=255)
                                                                                             tgtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAcc   <  1:289357/100‑1 (MQ=255)
                                                                                             tgtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAcc   >  1:2033000/1‑100 (MQ=255)
                                                                                             tgtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAcc   >  2:2184189/1‑100 (MQ=255)
                                                                                             tgtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAcc   >  1:1743322/1‑100 (MQ=255)
                                                                                             tgtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAcc   <  2:2268836/100‑1 (MQ=255)
                                                                                              gtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCa  >  1:1547337/1‑100 (MQ=255)
                                                                                                   |                                                                                               
ATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCA  >  NC_000913/3512461‑3512655

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