Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 896,520 T→C 6.4% P129P (CCT→CCC)  potH → putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913896,5200TC 6.4% 183.4 / 2.2 62P129P (CCT→CCC) potHputrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base C (2/2);  total (32/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 = 1.00e+00

TTAATCTCGGTAATTTTCTGCAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTGGGGGTTATCGATCAACCG  >  NC_000913/896375‑896664
                                                                                                                                                 |                                                                                                                                                
ttAATCTCGGTAATTTTCTGCAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTc                                                                                                                                                >  4:186737/1‑148 (MQ=255)
ttAATCTCGGTAATTTTCTGCAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCAGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTc                                                                                                                                                >  8:168101/1‑148 (MQ=255)
        ggTAATTTTCTGCAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTaa                                                                                                                                       >  1:119224/1‑149 (MQ=255)
        ggTAATTTTCTGCAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTaa                                                                                                                                       >  1:245411/1‑149 (MQ=255)
                   gCAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATtac                                                                                                                            >  4:327031/1‑149 (MQ=255)
                    cAACTTACCGACGATCCGCTCTATTTCGATTCTTATCTCCAGTCGCTCCAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATtac                                                                                                                            <  1:56598/148‑1 (MQ=255)
                    cAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATtac                                                                                                                            <  8:148037/148‑1 (MQ=255)
                    cAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGGGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTttatta                                                                                                                             <  5:227127/147‑1 (MQ=255)
                     aaCTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATtacta                                                                                                                          >  3:200440/1‑149 (MQ=255)
                          aCCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCTGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGt                                                                                                                     >  4:103439/1‑149 (MQ=255)
                          aCCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTgg                                                                                                                      <  8:141658/148‑1 (MQ=255)
                          aCCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTg                                                                                                                       <  2:119224/147‑1 (MQ=255)
                            cgacgaTCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGa                                                                                                                   >  3:188716/1‑149 (MQ=255)
                             gacgaTCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGAt                                                                                                                  >  7:117749/1‑149 (MQ=255)
                              acgaTCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATc                                                                                                                 <  2:64041/149‑1 (MQ=255)
                                 aTCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTg                                                                                                              <  2:169039/149‑1 (MQ=255)
                                     gCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCcgt                                                                                                          >  5:268973/1‑149 (MQ=255)
                                           tttCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTgg                                                                                                     <  1:123552/148‑1 (MQ=255)
                                           tttCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCCTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGa                                                                                                    >  6:36768/1‑149 (MQ=255)
                                           tttCGATGCTTATCTCCAGTCGCTACAGGGGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGa                                                                                                    <  1:192895/149‑1 (MQ=255)
                                                   cTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGttt                                                                                            <  4:200440/149‑1 (MQ=255)
                                                        ctcCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCATGCCGTCTTGGCCATCGTTTCTGAt                                                                                       >  8:148498/1‑149 (MQ=255)
                                                         tcCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATc                                                                                      >  8:114382/1‑149 (MQ=255)
                                                                    cAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGc                                                                           <  7:114382/149‑1 (MQ=255)
                                                                    cAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGc                                                                           >  3:56219/1‑149 (MQ=255)
                                                                      gggggCTGCGATTTCGACATTTTGCTGGTTACTGATCGGCTATCCGCTGGCCGGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTTGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCt                                                                         <  7:148498/146‑1 (MQ=255)
                                                                      ggTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTGATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCt                                                                         <  7:104708/149‑1 (MQ=255)
                                                                       gTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTg                                                                        >  2:85507/1‑149 (MQ=255)
                                                                            tGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGAt                                                                     <  7:53119/147‑1 (MQ=255)
                                                                            tGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGAt                                                                     >  8:53119/1‑147 (MQ=255)
                                                                            tGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATgg                                                                   >  2:189893/1‑149 (MQ=255)
                                                                            tGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATgg                                                                   >  3:306614/1‑149 (MQ=255)
                                                                            tGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATgg                                                                   >  6:47625/1‑149 (MQ=255)
                                                                              cGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGa                                                                 <  1:189893/149‑1 (MQ=255)
                                                                                        aTTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCTGCGTTTATGCCTGGATGGGAATATTaaaaa                                                       <  4:188716/149‑1 (MQ=255)
                                                                                         ttttGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGCCCCGTAATATTTTATTACGACTGGTGATCCTGCCGTCGTGGCCCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTaaaaaa                                                      >  8:368012/1‑149 (MQ=255)
                                                                                         ttttGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTaaaaaa                                                      >  1:89785/1‑149 (MQ=255)
                                                                                              cTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACg                                                 <  7:168101/149‑1 (MQ=255)
                                                                                               tGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACgg                                                >  5:248143/1‑149 (MQ=255)
                                                                                                 tttACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGgtg                                              >  4:169438/1‑149 (MQ=255)
                                                                                                     cTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGc                                           <  1:85507/148‑1 (MQ=255)
                                                                                                        aTCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGaa                                       >  4:95512/1‑149 (MQ=255)
                                                                                                        aTCGGCTATCCGCTGGCGTGGGCCGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTCGTGATCCTGCCGTCGTGGACCTCCTTTATGATCGGCGTTTATGCCTGGATGGGAATATTAAAAAACAACCGTGTGCTGaa                                       >  2:350611/1‑149 (MQ=255)
                                                                                                         tCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGaa                                       >  7:114334/1‑148 (MQ=255)
                                                                                                           ggCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGaataa                                    >  4:323400/1‑149 (MQ=255)
                                                                                                             cTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCCTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAAtt                                  >  1:205544/1‑149 (MQ=255)
                                                                                                                tCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTc                               <  3:327031/149‑1 (MQ=255)
                                                                                                                  cGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTctg                             >  5:82480/1‑149 (MQ=255)
                                                                                                                  cGCTGGCGTGGGCGGTGGCGCACAGTAAGCCCTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTctg                             <  2:205544/149‑1 (MQ=255)
                                                                                                                  cGCTGGCGTGGGCGGTGGCGCACAGTAAGCCCTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTct                              <  5:36768/148‑1 (MQ=255)
                                                                                                                           gggCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTg                    <  8:117749/149‑1 (MQ=255)
                                                                                                                              cGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTgggg                 <  3:323400/149‑1 (MQ=255)
                                                                                                                              cGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTgggg                 <  5:47625/149‑1 (MQ=255)
                                                                                                                               ggTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGATTAATTTTCTGCTGTGGCTggggg                <  3:95512/149‑1 (MQ=255)
                                                                                                                                    cgcACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGCCCTCGTTTCTGAACCGCGTTTATGCCGGGTGGGTAAAATTAAAAAACAACGGTGTGGTGAATAATTTTCTGCTGTGGCTGGGGGTTATc           <  6:268973/149‑1 (MQ=255)
                                                                                                                                    cgcACAGTAAGCCTTCGACAAGTAATATTTTATTACTACTTTTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGACTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTGGGGGTTATc           <  1:163207/149‑1 (MQ=255)
                                                                                                                                     gcACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATGCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTGGGGGTTATc           >  4:329615/1‑148 (MQ=255)
                                                                                                                                         aGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGttt                                                                               <  2:248801/76‑1 (MQ=255)
                                                                                                                                         aGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGttt                                                                               >  1:248801/1‑76 (MQ=255)
                                                                                                                                          gTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTGGGGGTTATCGATCaa     <  4:56219/149‑1 (MQ=255)
                                                                                                                                           tAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTGGGGGTTATCGATCAAc    >  6:139485/1‑149 (MQ=255)
                                                                                                                                             aGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTGGGGGTTATCGATCAACCg  <  5:139485/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTAATCTCGGTAATTTTCTGCAACTGACCGACGATCCGCTCTATTTCGATGCTTATCTCCAGTCGCTACAGGTGGCTGCGATTTCGACATTTTGCTGTTTACTGATCGGCTATCCGCTGGCGTGGGCGGTGGCGCACAGTAAGCCTTCGACCCGTAATATTTTATTACTACTGGTGATCCTGCCGTCGTGGACCTCGTTTCTGATCCGCGTTTATGCCTGGATGGGAATATTAAAAAACAACGGTGTGCTGAATAATTTTCTGCTGTGGCTGGGGGTTATCGATCAACCG  >  NC_000913/896375‑896664

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: ‑