Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,081,526 G→A 6.3% A37V (GCC→GTC)  tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,081,5260GA 6.3% 203.0 / 3.4 63A37V (GCC→GTC) tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base A (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTT  >  NC_000913/3081388‑3081671
                                                                                                                                          |                                                                                                                                                 
ggCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAg                                                                                                                                          <  3:87305/148‑1 (MQ=255)
  cAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGc                                                                                                                                         >  7:56950/1‑147 (MQ=255)
     aTCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCa                                                                                                                                                 <  7:157257/136‑1 (MQ=255)
     aTCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCa                                                                                                                                                 >  8:157257/1‑136 (MQ=255)
     aTCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGcc                                                                                                                                        >  3:333721/1‑145 (MQ=255)
     aTCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGcc                                                                                                                                        <  4:333721/145‑1 (MQ=255)
        gTAACCGGGGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATGCTGCGGGGTGTGTGTCAGGAACTCACGCCACAGGACTTCGGCAATGTCAGCCATAccc                                                                                                                                  <  2:156735/148‑1 (MQ=255)
         tAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGcc                                                                                                                                        >  7:267733/1‑141 (MQ=255)
         tAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGcc                                                                                                                                        <  8:267733/141‑1 (MQ=255)
           aCCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATAcc                                                                                                                                   >  8:317374/1‑144 (MQ=255)
           aCCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATAcc                                                                                                                                   <  7:317374/144‑1 (MQ=255)
           aCCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCa                                                                                                                                 >  6:175028/1‑146 (MQ=255)
           aCCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCAt                                                                                                                                <  7:366847/147‑1 (MQ=255)
           aCCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATa                                                                                                                               <  6:155205/148‑1 (MQ=255)
           aCCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATa                                                                                                                               <  5:175028/148‑1 (MQ=255)
              ggTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAggg                                                                                                                            >  3:307173/1‑148 (MQ=255)
                      gcagcaGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCgg                                                                                                                    >  8:162838/1‑148 (MQ=255)
                                 tAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGAtt                                                                                                         <  4:307173/148‑1 (MQ=255)
                                     tCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGc                                                                                                                          <  7:280047/127‑1 (MQ=255)
                                     tCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGc                                                                                                                          >  8:280047/1‑127 (MQ=255)
                                     tCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGc                                                                                                     >  8:242911/1‑148 (MQ=255)
                                     tCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGc                                                                                                     >  2:105153/1‑148 (MQ=255)
                                      cAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCt                                                                                                    >  7:303669/1‑148 (MQ=255)
                                      cAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCt                                                                                                    >  4:299338/1‑148 (MQ=255)
                                        gCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCttt                                                                                                  >  1:381819/1‑148 (MQ=255)
                                         cATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTc                                                                                                 <  7:242911/148‑1 (MQ=255)
                                         cATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTc                                                                                                 <  3:179884/148‑1 (MQ=255)
                                            ggAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATAc                                                                                                                                    <  7:346868/110‑1 (MQ=255)
                                            ggAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATAc                                                                                                                                    >  8:346868/1‑110 (MQ=255)
                                            ggAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGAt                                                                                                          <  7:237336/136‑1 (MQ=255)
                                            ggAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGAt                                                                                                          >  8:237336/1‑136 (MQ=255)
                                             gAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGt                                                                                              <  8:56950/147‑1 (MQ=255)
                                             gAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTa                                                                                             <  1:105153/148‑1 (MQ=255)
                                                  gTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCAt                                                                                                                                      <  6:81564/102‑1 (MQ=255)
                                                  gTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCAt                                                                                                                                      >  5:81564/1‑102 (MQ=255)
                                                    ggCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTa                                                                                             <  4:306199/141‑1 (MQ=255)
                                                    ggCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTa                                                                                             >  3:306199/1‑141 (MQ=255)
                                                    ggCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGACa                                                                                                                                                 >  2:131855/1‑89 (MQ=255)
                                                    ggCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGACa                                                                                                                                                 <  1:131855/89‑1 (MQ=255)
                                                     gCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTcc                                                                                     <  5:181327/148‑1 (MQ=255)
                                                          tGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCt                                                                                >  5:165572/1‑148 (MQ=255)
                                                                gCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAgcg                                                                           >  1:114993/1‑147 (MQ=255)
                                                                 cACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGCGGCACCCGGGTGACCGGATTTGGCTTTc                                                                                                 >  1:62221/1‑124 (MQ=255)
                                                                 cACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGCGGCACCCGGGTGACCGGATTTGGCTTTc                                                                                                 <  2:62221/124‑1 (MQ=255)
                                                                  aCGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCgg                                                                                                                    >  8:5869/1‑104 (MQ=255)
                                                                  aCGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCgg                                                                                                                    <  7:5869/104‑1 (MQ=255)
                                                                          gtcacggtcaGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCa                                                                <  2:381819/148‑1 (MQ=255)
                                                                                       cAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCtt                                                   >  4:24824/1‑148 (MQ=255)
                                                                                        aGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTa                                                                                             <  7:192233/105‑1 (MQ=255)
                                                                                        aGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTa                                                                                             >  8:192233/1‑105 (MQ=255)
                                                                                        aGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCtt                                                   >  8:286506/1‑147 (MQ=255)
                                                                                          gACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTAc                                                <  8:303669/148‑1 (MQ=255)
                                                                                           aCGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACg                                               <  7:286506/148‑1 (MQ=255)
                                                                                             ggATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTg                                             <  7:162838/148‑1 (MQ=255)
                                                                                              gATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGACAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCa                                                                <  7:192982/128‑1 (MQ=255)
                                                                                              gATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGACAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCa                                                                >  8:192982/1‑128 (MQ=255)
                                                                                                    gCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAt                                      >  4:207191/1‑148 (MQ=255)
                                                                                                    gCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAt                                      >  4:326830/1‑148 (MQ=255)
                                                                                                    gCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTAAGTGAGGACAt                                      >  2:317944/1‑148 (MQ=255)
                                                                                                        gTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTg                                  <  2:114993/148‑1 (MQ=255)
                                                                                                                 cAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTc                              >  3:158832/1‑143 (MQ=255)
                                                                                                                   ggAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCTTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCg                       <  6:165572/148‑1 (MQ=255)
                                                                                                                            gCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGaa              <  3:299338/148‑1 (MQ=255)
                                                                                                                                         ggCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCCGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGAGGATGAAGGGCACGCCCtt  >  8:237403/1‑147 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
GGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCATGGAGCCGTGGCCGTTGGACAGCACGAAGCGGTCACGGTCAGCCCAGGACGGATTCTGCGGGTTGTGTTTCAGGAAATCACGCCACAGGACTTCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTT  >  NC_000913/3081388‑3081671

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