Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,024,149 T→A 5.9% L199Q (CTG→CAG)  tatD → quality control of Tat‑exported FeS proteins; Mg‑dependent cytoplasmic DNase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,024,1490TA 5.9% 212.0 / 3.2 68L199Q (CTG→CAG) tatDquality control of Tat‑exported FeS proteins; Mg‑dependent cytoplasmic DNase
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base A (2/2);  total (35/33)
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.36e-01

CATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCAC  >  NC_000913/4024007‑4024296
                                                                                                                                              |                                                                                                                                                   
cATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCAATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt                                                                                                                                                  <  3:128057/146‑1 (MQ=255)
    gCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAAc                                                                                                                                            <  7:265806/148‑1 (MQ=255)
            ggCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGAGGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGcc                                                                                                                                   <  5:93447/149‑1 (MQ=255)
                 acGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc                                                                                                                              <  1:36424/149‑1 (MQ=255)
                 acGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACAGGGTGGGGTTGCGATGATCGACGCTGACTGGAGCTGCGGGAACTTTTGCCTTTGATTCCGGCGGAAAAATTACTGAGAGATACTGATGAGCCGTATc                                                                                                                              >  6:51698/1‑149 (MQ=255)
                     gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCt                                                                                                                          >  3:324790/1‑149 (MQ=255)
                       agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGcc                                                                                                                                   <  5:18887/138‑1 (MQ=255)
                       agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGcc                                                                                                                                   >  6:18887/1‑138 (MQ=255)
                           aTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATGGAAACTGATGCGCCGTATCTGCTCCCTc                                                                                                                     <  6:313950/148‑1 (MQ=255)
                              cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGa                                                                                                                 <  6:400422/149‑1 (MQ=255)
                               aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGAt                                                                                                                <  2:370752/149‑1 (MQ=255)
                                 gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGAtct                                                                                                              >  8:13308/1‑149 (MQ=255)
                                     gcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACg                                                                                                          >  6:79764/1‑149 (MQ=255)
                                         ggCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGGCGTTGATTCCGGCGGATAAATTACAGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGcaa                                                                                                       >  5:129416/1‑146 (MQ=255)
                                            gcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGc                                                                                                   >  1:358157/1‑149 (MQ=255)
                                               tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCaaa                                                                                                     >  1:24411/1‑144 (MQ=255)
                                               tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCaaa                                                                                                     <  2:24411/144‑1 (MQ=255)
                                               tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCcat                                                                                                >  2:177638/1‑149 (MQ=255)
                                                ggAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTTCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCcat                                                                                                >  5:162609/1‑148 (MQ=255)
                                                 gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCcatc                                                                                               <  2:358157/148‑1 (MQ=255)
                                                  aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTc                                                                                                                     <  5:385445/125‑1 (MQ=255)
                                                  aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTc                                                                                                                     >  6:385445/1‑125 (MQ=255)
                                                          tCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGgcgc                                                                                     <  2:263917/149‑1 (MQ=255)
                                                            ggCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCaa                                                                                   <  7:13308/149‑1 (MQ=255)
                                                              cATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACg                                                                                 <  3:344212/149‑1 (MQ=255)
                                                                ttACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGag                                                                               >  7:48759/1‑149 (MQ=255)
                                                                ttACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGag                                                                               >  7:157209/1‑149 (MQ=255)
                                                                ttACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGag                                                                               >  1:371413/1‑149 (MQ=255)
                                                                 tACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGagc                                                                              >  7:154004/1‑149 (MQ=255)
                                                                    cGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACAGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGagccag                                                                           >  1:233672/1‑149 (MQ=255)
                                                                      gTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGagccagcc                                                                         >  8:23130/1‑149 (MQ=255)
                                                                               gCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTTACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGcc                                                                 <  3:371927/148‑1 (MQ=255)
                                                                                  aTGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt                                                                                                                                                  <  2:40749/64‑1 (MQ=255)
                                                                                  aTGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt                                                                                                                                                  >  1:40749/1‑64 (MQ=255)
                                                                                  aTGAACGACGCGGACTGGAGCTGCGGGAACCTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCa                                                              >  2:127142/1‑148 (MQ=255)
                                                                                      acgacgCGGACGGGAGCTGCGGGAACGTTTGCCGTTGATTCCGGCGGAAAAAGTACTGATCGACACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCCTCTGCCCCAGAttt                                                         >  6:77711/1‑149 (MQ=255)
                                                                                       cgacgCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATAtttt                                                        <  8:157209/149‑1 (MQ=255)
                                                                                       cgacgCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATAtttt                                                        <  5:79764/149‑1 (MQ=255)
                                                                                           gcgGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCaa                                                    >  8:167610/1‑149 (MQ=255)
                                                                                           gcgGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCaa                                                    <  8:48759/149‑1 (MQ=255)
                                                                                                cTGGAGGTGCGGGAACTATTGCCGTTGATTGCGGCGTAAAAATTACAGATCGAAACTGATGCGCCGTATGTGTTCCGTCGCGAGCTCACGCCAAAGCCATGATCCCGGCGCAAGGTGCCAGGCCATCTGCCCCATATTTTGCACCGttt                                               <  8:149504/149‑3 (MQ=255)
                                                                                                   gAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTAtt                                              <  1:177638/147‑1 (MQ=255)
                                                                                                      cTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCa                                         <  8:154004/149‑1 (MQ=255)
                                                                                                        gCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACt                                       >  5:149072/1‑149 (MQ=255)
                                                                                                         cGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTg                                      >  3:155774/1‑149 (MQ=255)
                                                                                                          gggAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTgg                                     >  8:93723/1‑149 (MQ=255)
                                                                                                           ggAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGc                                    >  1:37172/1‑149 (MQ=255)
                                                                                                            gAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTACCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCg                                   >  3:239214/1‑149 (MQ=255)
                                                                                                              aCTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGTCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTg                                 >  1:114069/1‑149 (MQ=255)
                                                                                                              aCTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTg                                 <  7:23130/149‑1 (MQ=255)
                                                                                                              aCTTTTGCAGTTGATTCCGGCGGAAAAATTACAGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGAGCTCCCGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTg                                 <  6:129416/149‑1 (MQ=255)
                                                                                                                     ccGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGAGCTCCCGCCAAAGCCATCATCACGTCGCCACGCGCCAGACCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGagaagc                          >  5:355356/1‑148 (MQ=255)
                                                                                                                      cGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTAtt                                              <  1:182721/128‑1 (MQ=255)
                                                                                                                      cGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACCCCAAAGCCCACATCCCGCCGCAACGAGCCAGCCCCTCTGCCCCATAATTTGCAACgaaat                                              >  2:182721/1‑124 (MQ=255)
                                                                                                                           aTTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCg                     <  6:149072/148‑1 (MQ=255)
                                                                                                                            ttCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCCACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCa                   >  2:302077/1‑149 (MQ=255)
                                                                                                                              ccGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCAt                  <  1:11790/148‑1 (MQ=255)
                                                                                                                               cggcggAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCAtgg                <  2:37172/149‑1 (MQ=255)
                                                                                                                                ggcggAAAAATTACTGATAGAAACTGATGCGACGTATATGCTAACTCGCGATCTCACGCCAAAGACATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCAATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCAtggc               >  2:234440/1‑149 (MQ=255)
                                                                                                                                 gcggAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCAtggct              <  3:134509/149‑1 (MQ=255)
                                                                                                                                 gcggAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCAtggct              <  3:128795/149‑1 (MQ=255)
                                                                                                                                   ggAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCAtggctgg            <  4:324790/149‑1 (MQ=255)
                                                                                                                                    gAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCAtggctggc           >  4:97448/1‑149 (MQ=255)
                                                                                                                                    gAAAAATTACTGATCGAAACTGATGCGCCGTATCAGCTCCCTCGCGATCTCACGCCAAAGCCATCATACCGGAGCAACGAGCCAGCCCATCTGCCCCAAATTTTCCAACGTATTTCGCCCTTGCGTGGAGAAGATGCCGCAtggctggc           >  4:63924/1‑149 (MQ=255)
                                                                                                                                     aaaaaTTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGTCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATggctggct          <  2:114069/149‑1 (MQ=255)
                                                                                                                                       aaaTTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGc        <  7:167610/149‑1 (MQ=255)
                                                                                                                                        aaTTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGcc       >  8:177775/1‑149 (MQ=255)
                                                                                                                                             cTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGACCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCcaccac  >  1:218879/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCTCGCGATCTCACGCCAAAGCCATCATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCAC  >  NC_000913/4024007‑4024296

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