Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,705 C→T 9.6% intergenic (+24/‑204) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7050CT 9.6% 348.2 / 19.0 116intergenic (+24/‑204)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (49/54);  new base T (5/6);  total (55/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

ACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTAT  >  NC_000913/3334559‑3334847
                                                                                                                                                  |                                                                                                                                               
  gcgTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                            >  4:13452/1‑149 (MQ=255)
     tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCggg                                                                                                                                         >  4:626353/1‑149 (MQ=255)
      cAGCGTGCCGAGGAGGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCggg                                                                                                                                         <  2:148968/148‑1 (MQ=255)
       aGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGAt                                                                                                                                       <  3:626353/149‑1 (MQ=255)
         cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGATga                                                                                                                                     <  4:78118/149‑2 (MQ=255)
          gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGt                                                                                                                                    >  3:87641/1‑149 (MQ=255)
              cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTcc                                                                                                                                 >  1:91249/1‑148 (MQ=255)
              cGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCa                                                                                                                                <  4:55084/149‑1 (MQ=255)
               gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCa                                                                                                                                <  2:383705/148‑1 (MQ=255)
               gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAg                                                                                                                               <  3:872842/149‑1 (MQ=255)
               gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAg                                                                                                                               <  4:835123/149‑1 (MQ=255)
                  gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTaa                                                                                                                            >  3:792259/1‑149 (MQ=255)
                   aagaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTaa                                                                                                                            <  2:91249/148‑1 (MQ=255)
                          aGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                    >  3:581278/1‑149 (MQ=255)
                           gACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                    >  2:889742/1‑148 (MQ=255)
                            aCAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTcc                                                                                                                                 <  1:647019/134‑1 (MQ=255)
                            aCAAAACCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTcc                                                                                                                                 >  2:647019/1‑134 (MQ=255)
                              aaaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                                                 <  2:306012/148‑1 (MQ=255)
                              aaaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaaa                                                                                                                 <  1:773637/148‑1 (MQ=255)
                               aaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTcc                                                                                                                                 <  1:429463/131‑1 (MQ=255)
                               aaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTcc                                                                                                                                 >  2:429463/1‑131 (MQ=255)
                               aaGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAac                                                                                                               <  4:792259/149‑1 (MQ=255)
                                aGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAaca                                                                                                              <  3:903165/149‑1 (MQ=255)
                                    aTCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTa                                                                                                          >  4:490233/1‑149 (MQ=255)
                                     tCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                         >  4:440617/1‑149 (MQ=255)
                                      cGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                        <  3:568193/149‑1 (MQ=255)
                                         aGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAg                                                                                                     <  3:274025/149‑1 (MQ=255)
                                                 aGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                             >  3:195722/1‑149 (MQ=255)
                                                  ggTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                             <  2:34227/148‑1 (MQ=255)
                                                  ggTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCA‑‑‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAacac                                                                                          <  4:184883/149‑1 (MQ=255)
                                                       tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGa                                                                                                                                        >  4:481627/1‑100 (MQ=255)
                                                       tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCt                                                                                        >  1:557603/1‑148 (MQ=255)
                                                       tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                       >  4:42784/1‑149 (MQ=255)
                                                       tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGC‑GCTGCGCGGGa                                                                                                                                        <  3:481627/100‑1 (MQ=255)
                                                        cccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                       >  2:373936/1‑148 (MQ=255)
                                                         ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATg                                                                                                                                      >  1:676817/1‑100 (MQ=255)
                                                         ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATg                                                                                                                                      <  2:676817/100‑1 (MQ=255)
                                                         ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGa                                                                                      >  2:157714/1‑148 (MQ=255)
                                                               acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTc                                                                                                                                  <  4:582848/98‑1 (MQ=255)
                                                               acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                               >  3:798121/1‑149 (MQ=255)
                                                                ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGATg                                                                                                                                      <  1:690784/93‑1 (MQ=255)
                                                                ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGATg                                                                                                                                      >  2:690784/1‑93 (MQ=255)
                                                                ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                               <  1:13961/148‑1 (MQ=255)
                                                                  ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                            >  4:312982/1‑149 (MQ=255)
                                                                  ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                            <  4:195722/149‑1 (MQ=255)
                                                                  ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                            >  4:845976/1‑149 (MQ=255)
                                                                   cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                            <  2:557603/148‑1 (MQ=255)
                                                                   cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                           <  4:83828/149‑1 (MQ=255)
                                                                     tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgt                                                                         <  3:523569/149‑1 (MQ=255)
                                                                     tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgt                                                                         <  3:832679/149‑1 (MQ=255)
                                                                     tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                                          <  1:462126/148‑1 (MQ=255)
                                                                       gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGc                                                                       <  3:13452/149‑1 (MQ=255)
                                                                             aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAACCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCtt                                                                  <  2:651944/148‑1 (MQ=255)
                                                                               cACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtat                                                               <  4:627646/149‑1 (MQ=255)
                                                                                aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                               <  3:416409/132‑1 (MQ=255)
                                                                                aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                               >  4:416409/1‑132 (MQ=255)
                                                                                aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGAGATGTGCATCTTtata                                                              >  3:573657/1‑149 (MQ=255)
                                                                                  tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                        <  1:618439/137‑1 (MQ=255)
                                                                                  tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                        >  2:618439/1‑137 (MQ=255)
                                                                                  tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                            >  3:403580/1‑149 (MQ=255)
                                                                                   cATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatata                                                            >  1:295683/1‑148 (MQ=255)
                                                                                    aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                          >  4:531183/1‑149 (MQ=255)
                                                                                     tCGGCGTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                         <  3:149089/149‑1 (MQ=255)
                                                                                     tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                          >  2:168541/1‑148 (MQ=255)
                                                                                       ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTg                                                       <  4:87641/149‑1 (MQ=255)
                                                                                         ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                     <  3:490233/149‑1 (MQ=255)
                                                                                         ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGa                                                      <  1:719694/148‑1 (MQ=255)
                                                                                          cTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                     <  1:889742/148‑1 (MQ=255)
                                                                                             gcgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAatat                                                  <  1:157714/148‑1 (MQ=255)
                                                                                                    tCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAACTTCCATAAAAACACTTATTGAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCGCtctc                                          <  4:573657/149‑1 (MQ=255)
                                                                                                       cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCTGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTc                                                >  1:661034/1‑140 (MQ=255)
                                                                                                       cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCTGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTc                                                <  2:661034/140‑1 (MQ=255)
                                                                                                         gTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTAc                                     <  4:581278/149‑1 (MQ=255)
                                                                                                             aaCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg                                   <  3:88329/147‑1 (MQ=255)
                                                                                                             aaCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg                                   >  4:88329/1‑147 (MQ=255)
                                                                                                              aCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGaa                                >  3:730291/1‑149 (MQ=255)
                                                                                                              aCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGaa                                <  4:670099/149‑1 (MQ=255)
                                                                                                              aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCggg                                                                                                                                         >  1:284221/1‑44 (MQ=18)
                                                                                                              aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCggg                                                                                                                                         <  2:284221/44‑1 (MQ=18)
                                                                                                                gcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTa                                      >  1:739974/1‑141 (MQ=255)
                                                                                                                    tCGTTAATACCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                          >  3:635756/1‑149 (MQ=255)
                                                                                                                     cGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                         <  4:635756/149‑1 (MQ=255)
                                                                                                                       ttAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTag                        >  1:34046/1‑148 (MQ=255)
                                                                                                                        tAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTagag                      <  3:312982/149‑1 (MQ=255)
                                                                                                                           tcccctccccTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCaa                   >  4:279176/1‑149 (MQ=255)
                                                                                                                             ccctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAAt                  >  1:616401/1‑148 (MQ=255)
                                                                                                                              cctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                 >  1:413044/1‑148 (MQ=255)
                                                                                                                               ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcac                                              <  1:649137/118‑1 (MQ=255)
                                                                                                                               ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcac                                              >  2:649137/1‑118 (MQ=255)
                                                                                                                               ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                       >  3:243461/1‑125 (MQ=255)
                                                                                                                               ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                       <  4:243461/125‑1 (MQ=255)
                                                                                                                               ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAg                >  2:820800/1‑148 (MQ=255)
                                                                                                                                   cctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCAt           >  4:431625/1‑149 (MQ=255)
                                                                                                                                      cATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                          <  2:544009/131‑1 (MQ=255)
                                                                                                                                      cATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                          >  1:544009/1‑131 (MQ=255)
                                                                                                                                      cATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa         >  2:737285/1‑148 (MQ=255)
                                                                                                                                       aTCCCGCGCGCT‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                          >  1:21542/1‑130 (MQ=255)
                                                                                                                                       aTCCCGCGC‑GCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg        >  2:599071/1‑148 (MQ=255)
                                                                                                                                         cccGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGa      >  1:76857/1‑148 (MQ=255)
                                                                                                                                             cgcAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGtat  <  1:490932/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTAT  >  NC_000913/3334559‑3334847

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