Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 89,754 T→C Y41H (TAT→CAT)  mraZ → DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091389,7540TC97.3% 115.6 / ‑4.1 37Y41H (TAT→CAT) mraZDNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (18/18);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

AACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCTGTGGGA  >  NC_000913/89514‑89995
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
aaCATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGGAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTAAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGCACAGCTGCTTGAGAAAGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTg                                                                                                                                                                                                                                         >  1:247676/1‑251 (MQ=255)
       ggTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCctgctg                                                                                                                                                                                                                                  >  1:314680/1‑251 (MQ=255)
                     gTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAACGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGATTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCCTTTACCTTCATCACCCGTGCCTGCTGCTTTAccccattcc                                                                                                                                                                                                                    >  2:140609/1‑246 (MQ=255)
                      tATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCcctgcct                                                                                                                                                                                                                   >  1:203728/1‑251 (MQ=255)
                                           cAAGCTTTGCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCaa                                                                                                                                                                                              <  1:142907/251‑1 (MQ=255)
                                                   tcctccGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGGTAGTCAATCTCGACAGCAAAGGGCGCTTAGCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTAt                                                                                                                                                                                       <  1:35983/250‑1 (MQ=255)
                                                                aaaCTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCcctgcct                                                                                                                                                                                                                   >  2:26884/1‑209 (MQ=255)
                                                                aaaCTCCGTGCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAg                                                                                                                                                                          <  2:55659/250‑1 (MQ=255)
                                                                     ccTTTCAGGGGGAAATTGGGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATg                                                                                                                                                                      <  2:314680/249‑1 (MQ=255)
                                                                                      gtgGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCTTGTTCCGGGGATCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTACCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGTTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTAGCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtg                                                                                                                                                   <  2:365864/251‑1 (MQ=255)
                                                                                      gtgGGGAAAAGTGTGAATAGGGGGTGAGGCTTGCATGTGCCGGGGAGTACCTTTTGACAATCTCGACATCAAAGGGCGCTTATCAGTGCCTGCCCGTGGCCGGGAACAGCTGCTGGAGAACGCTGCCGGTCAAATGGTTTGCAGCGTTGACAGTCATCATCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAATTTAGCGAGTAGAAATTATCGCGTCTGGTGAGCATTAACCCGGGTGAGCGCtgtg                                                                                                                                                   <  2:59063/251‑1 (MQ=255)
                                                                                             aaaGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgcg                                                                                                                                            >  1:90362/1‑251 (MQ=255)
                                                                                                   ggAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACt                                                                                                                                      >  1:386207/1‑251 (MQ=255)
                                                                                                     aaTAAGGGGTGAGGCTTGCATGTTTCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGt                                                                                                                                    <  2:173645/251‑1 (MQ=255)
                                                                                                                      gCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGa                                                                                                                    >  1:187635/1‑250 (MQ=255)
                                                                                                                        aTGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg                                                                                                                 >  1:186590/1‑251 (MQ=255)
                                                                                                                                 ggAGCAACGTTAGTCAATCTCGACAGCACAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAACTGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATAGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGACTGTCAgctgg                                                                                                         >  1:281279/1‑250 (MQ=255)
                                                                                                                                 ggAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTGTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCAGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatgga                                                                                                        <  2:364381/251‑1 (MQ=255)
                                                                                                                                            cgTCAATCTCGACAGCATAGGTCGCTTATCAGTTCCTACCCGTTACCCGGAACAGCTTCTTGAGAACGCTGCCGGTCAAATGGTGTTCACCATTGACATGCATCTCCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATTAACCCGGTTGAGCGGCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTc                                                                                             <  2:4279/250‑1 (MQ=255)
                                                                                                                                               cAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAggaggat                                                                                          >  2:68207/1‑246 (MQ=255)
                                                                                                                                                     cGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGATTGTCAGATGGATGgc                                                                                                    >  2:315380/1‑235 (MQ=255)
                                                                                                                                                     cGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgc                                                                                                    <  1:315380/235‑1 (MQ=255)
                                                                                                                                                       aCAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGAATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATc                                                                                  >  2:310723/1‑251 (MQ=255)
                                                                                                                                                          gCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATcgcg                                                                               <  1:241963/251‑1 (MQ=255)
                                                                                                                                                                  cgcTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACt                                                                       >  1:222910/1‑251 (MQ=255)
                                                                                                                                                                          aGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAAc                                                               >  1:117508/1‑251 (MQ=255)
                                                                                                                                                                                        cccGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACa                                                 <  2:201376/249‑1 (MQ=255)
                                                                                                                                                                                         ccGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACa                                                 >  1:327997/2‑250 (MQ=255)
                                                                                                                                                                                                     gctTGAGAACGCGGCCGGTCAAATGGTGTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGAt                                      <  2:90362/249‑1 (MQ=255)
                                                                                                                                                                                                      ctTGAGAACGCTGACGGCCAAGTGGTTTGCACCATTGACATTCGTCACCCCTGCCTGCTTCCTTACCCCCGGCATGAATGGGAAATTATCGAGCCAAAATTATCGCGTCCGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCTACCACATGACGGGCTGAGAAAAGAAGTGATGCt                                   >  2:197871/1‑251 (MQ=255)
                                                                                                                                                                                                         gagaACGCTTCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCTGGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTgg                                 <  1:321930/250‑1 (MQ=255)
                                                                                                                                                                                                                 tGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCTCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGAcggt                        >  2:220278/1‑248 (MQ=255)
                                                                                                                                                                                                                         aaaTGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAg                <  1:310723/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggTTTGCACCATTGACATTCAGCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCGACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCTGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGttt             <  1:18535/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gccccATTGACATTCAACCCCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCtt       <  2:333860/248‑2 (MQ=255)
                                                                                                                                                                                                                                    aCCACTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCCCGTCTGTCTAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGTCGCAGGTCGATTGTTAATCGCGCCAGTACTTCGGAAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGAAAGTTCAACAAGATTGAGCtgcg     >  2:78181/1‑249 (MQ=255)
                                                                                                                                                                                                                                       aTTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCTGTGGGa  <  2:186590/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
AACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGATTGTTAATCGCGCCAGTACTGCGGCAACATGCCGGGCTGACAAAAGAAGTGATGCTGGTTGGACAGTTCAACAAGTTTGAGCTGTGGGA  >  NC_000913/89514‑89995

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