Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,715,919 T→G 6.6% intergenic (‑30/+29) pdxY ← / ← tyrS pyridoxamine kinase/tyrosyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,715,9190TG 6.6% 206.5 / 6.9 61intergenic (‑30/+29)pdxY/tyrSpyridoxamine kinase/tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base G (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

TGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAAT  >  NC_000913/1715771‑1716067
                                                                                                                                                    |                                                                                                                                                    
tGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCACATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAt                                                                                                                                                      <  4:432919/149‑1 (MQ=255)
  aaCGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAttt                                                                                                                                                    >  4:615465/1‑149 (MQ=255)
   aCGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTtc                                                                                                                                                   <  4:629925/149‑1 (MQ=255)
   aCGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTtc                                                                                                                                                   <  2:347964/149‑1 (MQ=255)
            aGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTcccc                                                                                                                                          <  2:469295/149‑1 (MQ=255)
             gCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCa                                                                                                                                                >  1:51185/1‑142 (MQ=255)
             gCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCa                                                                                                                                                <  2:51185/142‑1 (MQ=255)
                   aGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCa                                                                                                                                   >  2:718971/1‑149 (MQ=255)
                     aCGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACt                                                                                                                                 >  4:148328/1‑149 (MQ=255)
                        ttCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTaa                                                                                                                              >  3:254054/1‑149 (MQ=255)
                             gcCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCat                                                                                                                         >  1:292140/1‑149 (MQ=255)
                              cccAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTcccc                                                                                                                                          >  1:714472/1‑131 (MQ=255)
                              cccAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTcccc                                                                                                                                          <  2:714472/131‑1 (MQ=255)
                              cccAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCatt                                                                                                                        >  2:85019/1‑149 (MQ=255)
                               ccAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCatta                                                                                                                       >  1:265754/1‑149 (MQ=255)
                                        gcATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCAc                                                                                                                                  <  3:122930/129‑1 (MQ=255)
                                        gcATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCAc                                                                                                                                  >  4:122930/1‑129 (MQ=255)
                                          aTCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGAAATCACTCCCCCtgg                                                                                                                                      >  3:109493/1‑121 (MQ=255)
                                          aTCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGAGAAATCACTCCCCCtgg                                                                                                                                      <  4:109493/123‑3 (MQ=255)
                                           tCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAAt                                                                                                           <  4:757534/149‑1 (MQ=255)
                                            cGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATc                                                                                                          <  1:85019/149‑1 (MQ=255)
                                             ggAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCa                                                                                                         >  3:265054/1‑149 (MQ=255)
                                              gAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAg                                                                                                        >  1:766824/1‑149 (MQ=255)
                                                   tctGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTcccc                                                                                                                                          <  3:402553/110‑1 (MQ=255)
                                                   tctGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTcccc                                                                                                                                          >  4:402553/1‑110 (MQ=255)
                                                   tctGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGt                                                                                                   >  4:220813/1‑149 (MQ=255)
                                                          cACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtt                                                                                            <  2:292140/149‑1 (MQ=255)
                                                            cTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCtttt                                                                                          <  2:252501/149‑1 (MQ=255)
                                                               ttACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTAc                                                                                       <  3:148328/149‑1 (MQ=255)
                                                                         gACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCa                                                                                                                                   <  1:635537/95‑1 (MQ=255)
                                                                         gACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCa                                                                                                                                   >  2:635537/1‑95 (MQ=255)
                                                                           ccATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTa                                                                           >  3:84127/1‑149 (MQ=255)
                                                                             aTAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                         >  2:146670/1‑149 (MQ=255)
                                                                                  caacGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaa                                                                    >  2:809588/1‑149 (MQ=255)
                                                                                      gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                         <  4:115911/140‑1 (MQ=255)
                                                                                      gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                         >  3:115911/1‑140 (MQ=255)
                                                                                      gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                         >  3:405004/1‑140 (MQ=255)
                                                                                      gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAg                                                                         <  4:405004/140‑1 (MQ=255)
                                                                                      gTGAGACTGGAAAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGa                                                                >  1:29729/1‑149 (MQ=255)
                                                                                           aCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCaaa                                                           >  1:335951/1‑149 (MQ=255)
                                                                                                 tAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTC‑‑TTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAt                                                   <  4:641136/149‑1 (MQ=255)
                                                                                                     gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtct                                                 >  4:185859/1‑149 (MQ=255)
                                                                                                     gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtct                                                 >  2:724005/1‑149 (MQ=255)
                                                                                                        aaTATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttct                                              <  2:335951/149‑1 (MQ=255)
                                                                                                                  atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGt                                                                            <  4:541928/109‑1 (MQ=255)
                                                                                                                  atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGt                                                                            >  3:541928/1‑109 (MQ=255)
                                                                                                                   tcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                                   <  3:615465/149‑1 (MQ=255)
                                                                                                                   tcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                                   <  3:220813/149‑1 (MQ=255)
                                                                                                                   tcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa                                   <  4:254054/149‑1 (MQ=255)
                                                                                                                            cTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAgg                          >  1:86451/1‑149 (MQ=255)
                                                                                                                             tGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa                         <  1:146670/149‑1 (MQ=255)
                                                                                                                              gTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGAt                        <  3:185859/149‑1 (MQ=255)
                                                                                                                              gTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGAt                        <  2:766824/149‑1 (MQ=255)
                                                                                                                              gTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGAt                        <  1:724005/149‑1 (MQ=255)
                                                                                                                               tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATc                       <  1:809588/149‑1 (MQ=255)
                                                                                                                                       ccAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGtt               >  3:272277/1‑149 (MQ=255)
                                                                                                                                        cAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGttt              <  4:272277/149‑1 (MQ=255)
                                                                                                                                        aaGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTAAAAGAAGTATTCAGGATCGGATTAttt              >  4:691894/2‑149 (MQ=255)
                                                                                                                                          gggggAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                                          <  4:154859/119‑1 (MQ=255)
                                                                                                                                          gggggAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                                          >  3:154859/1‑119 (MQ=255)
                                                                                                                                           ggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGTCCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTc           >  4:459050/1‑149 (MQ=255)
                                                                                                                                                   aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa                            <  1:453432/124‑1 (MQ=255)
                                                                                                                                                   aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa                            >  2:453432/1‑124 (MQ=255)
                                                                                                                                                    tttCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGTCCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAAt  <  3:459050/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAAT  >  NC_000913/1715771‑1716067

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