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.99.1% 470.1 / ‑2.6 109intergenic (‑216/+74)yhfZ/trpSputative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (59/49);  total (60/49)
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.46e-01

Only 100 of 113 total aligned reads displayed.

ATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCA  >  NC_000913/3512461‑3512655
                                                                                                   |                                                                                               
   aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtact                                                                                              >  2:2381392/1‑97 (MQ=255)
   aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtact                                                                                              >  1:1602519/1‑97 (MQ=255)
   aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtact                                                                                              >  1:2429727/1‑97 (MQ=255)
      aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTg                                                                                          >  2:2737797/1‑100 (MQ=255)
         ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAg                                                                                       >  2:924576/1‑100 (MQ=255)
          cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAgg                                                                                      <  2:2331328/100‑1 (MQ=255)
          cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAgg                                                                                      <  1:389743/100‑1 (MQ=255)
           gggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGc                                                                                     <  2:192297/100‑1 (MQ=255)
            ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGcc                                                                                    <  1:481540/100‑1 (MQ=255)
            ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGcc                                                                                    <  1:2381407/100‑1 (MQ=255)
                   ccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAg                                                                             <  2:1557684/100‑1 (MQ=255)
                    cGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAgg                                                                            <  2:445565/100‑1 (MQ=255)
                     ggTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGc                                                                           <  1:513939/100‑1 (MQ=255)
                      gTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCg                                                                          <  1:1601928/100‑1 (MQ=255)
                          cAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGAtgcggcgtgaa                                                                      <  1:79359/100‑11 (MQ=25)
                          cAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGAtgcggcgtgaa                                                                      >  1:2394810/1‑90 (MQ=25)
                          cAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCa                                                                      >  1:2042375/1‑100 (MQ=255)
                             agTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgc                                                                   <  2:631115/100‑1 (MQ=255)
                              gTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                  <  2:1206000/100‑1 (MQ=255)
                              gTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                  >  2:2811454/1‑100 (MQ=255)
                               tgatgaTAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGTCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                 <  2:1548847/100‑1 (MQ=255)
                               tgatgaTAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                 >  2:2001419/1‑100 (MQ=255)
                                gatgatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGAGCATACTTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                >  1:2026408/1‑100 (MQ=255)
                                 atgatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATGCGGCGTGAAcgcctta                                                               <  1:311211/100‑4 (MQ=21)
                                 atgatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                               >  1:912538/1‑100 (MQ=255)
                                   gatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                             >  2:370260/1‑100 (MQ=255)
                                    atAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                            >  1:1794891/1‑100 (MQ=255)
                                     tAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                           >  1:1936934/1‑100 (MQ=255)
                                     tAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                           >  1:909548/1‑100 (MQ=255)
                                     tAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATAATTTGTAGGCAGGTTTAGGCGTTCACGCCGCATCCg                                                           >  2:838500/1‑100 (MQ=25)
                                      aaaaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                          >  2:248014/1‑100 (MQ=255)
                                      aaaaGCAAAATTGCCGGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                          <  1:2840629/100‑1 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                        <  1:2472066/100‑1 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                        >  2:1751146/1‑100 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                        <  2:1967994/100‑1 (MQ=255)
                                          gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                      <  1:1031353/100‑1 (MQ=255)
                                          gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                      >  2:1256585/1‑100 (MQ=255)
                                          gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                      <  1:2119692/100‑1 (MQ=255)
                                          gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                      <  1:35020/100‑1 (MQ=255)
                                            aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                    >  1:1648535/1‑100 (MQ=255)
                                              aaTTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACa                                                  <  2:2835440/100‑1 (MQ=255)
                                               aTTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACaa                                                 <  2:2681461/100‑1 (MQ=255)
                                                 tGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAg                                               <  1:2104823/100‑1 (MQ=255)
                                                   ccTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAgcg                                             >  1:2447086/1‑100 (MQ=255)
                                                    cTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAgcgc                                            <  1:705945/100‑1 (MQ=255)
                                                     tGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAgcgc                                            >  2:2454810/1‑99 (MQ=255)
                                                     tGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCa                                           >  1:1033247/1‑100 (MQ=255)
                                                      gATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCaa                                          >  2:2663867/1‑100 (MQ=255)
                                                       aTGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAAt                                         <  1:2368982/100‑1 (MQ=255)
                                                          cgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTg                                      >  2:1220996/1‑100 (MQ=255)
                                                          cgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTg                                      <  1:2142783/100‑1 (MQ=255)
                                                           gcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGc                                     >  2:87287/1‑100 (MQ=255)
                                                            cTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGcc                                    <  1:1384223/100‑1 (MQ=255)
                                                            cTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGcc                                    >  2:2324160/1‑100 (MQ=255)
                                                              aCGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                  >  1:1535791/1‑100 (MQ=255)
                                                              aCGCTTATCAGGACTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAGGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                  >  1:1901235/1‑100 (MQ=25)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                 >  2:1965264/1‑100 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                 >  2:773556/1‑100 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                 <  2:1230390/100‑1 (MQ=255)
                                                                gCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCa                                >  2:1233228/1‑100 (MQ=255)
                                                                gCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCa                                >  2:140635/1‑100 (MQ=255)
                                                                 cTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCaa                               >  1:1928905/1‑100 (MQ=255)
                                                                   tATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATa                             >  1:2100741/1‑100 (MQ=255)
                                                                    aTCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAAt                              >  1:477998/1‑98 (MQ=255)
                                                                    aTCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAg                            >  1:934390/1‑100 (MQ=255)
                                                                      cAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTg                          <  2:909547/100‑1 (MQ=255)
                                                                      cAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTg                          <  2:22999/100‑1 (MQ=255)
                                                                         gCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAt                       <  1:1965278/100‑1 (MQ=255)
                                                                         gCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAt                       <  1:2387764/100‑1 (MQ=255)
                                                                         gCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAt                       >  1:2731283/1‑100 (MQ=255)
                                                                          ccTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAtt                      <  1:2811475/100‑1 (MQ=255)
                                                                          ccTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAtt                      >  1:1876548/1‑100 (MQ=255)
                                                                          ccTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAtt                      >  1:971041/1‑100 (MQ=255)
                                                                           cTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTa                     >  2:1869897/1‑100 (MQ=255)
                                                                             aCATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACg                   <  1:2737818/100‑1 (MQ=255)
                                                                              cATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACgg                  <  1:924578/100‑1 (MQ=255)
                                                                              cATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACgg                  <  1:664350/100‑1 (MQ=255)
                                                                              cATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACgg                  >  2:1007248/1‑100 (MQ=255)
                                                                                tttCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCt                >  2:1836439/1‑100 (MQ=255)
                                                                                 ttCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCtt               >  2:2538067/1‑100 (MQ=255)
                                                                                  tCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTc              <  1:1316780/100‑1 (MQ=255)
                                                                                  tCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTc              >  2:631220/1‑100 (MQ=255)
                                                                                   ccTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCg             <  2:989965/100‑1 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc            >  2:2439925/1‑100 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc            >  1:1583862/1‑100 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc            <  1:2339894/100‑1 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc            >  2:2589793/1‑100 (MQ=255)
                                                                                     ttGCAATAGGTTCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGcc           >  2:2740636/1‑100 (MQ=21)
                                                                                      tGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCca          <  1:848502/100‑1 (MQ=255)
                                                                                      tGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCca          >  2:2782712/1‑100 (MQ=255)
                                                                                      tGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCca          >  2:2764753/1‑100 (MQ=255)
                                                                                       gCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcac         >  1:544450/1‑100 (MQ=255)
                                                                                       gCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcac         >  1:2521265/1‑100 (MQ=255)
                                                                                        cAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcaca        <  1:1829450/100‑1 (MQ=255)
                                                                                          atatGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACaaa      <  2:2023202/100‑1 (MQ=255)
                                                                                          atatGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACaaa      <  1:337835/100‑1 (MQ=255)
                                                                                           tatGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACaaaa     >  2:265175/1‑100 (MQ=255)
                                                                                             tgtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAcc   >  2:2581584/1‑100 (MQ=255)
                                                                                             tgtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAcc   <  2:961390/100‑1 (MQ=255)
                                                                                              gtgCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCa  >  1:770874/1‑100 (MQ=255)
                                                                                                   |                                                                                               
ATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCA  >  NC_000913/3512461‑3512655

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