Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,078,568 A→G 5.9% F98L (TTC→CTC)  cmtA ← putative mannitol‑specific PTS IIB and IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,078,5680AG 5.9% 207.4 / 2.6 68F98L (TTC→CTC) cmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base G (2/2);  total (34/34)
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.93e-01

AGCAGACAAAGGAGCATCCCCGCGATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGACAATTTTGGCG  >  NC_000913/3078423‑3078713
                                                                                                                                                 |                                                                                                                                                 
agcagACAAAGGAGCATCCCCGCGATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACAt                                                                                                                                                >  2:329243/1‑149 (MQ=255)
     aCAAAGGAGCATCCCCGCGATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGa                                                                                                                                           <  5:360746/149‑1 (MQ=255)
     aCAAAGGAGCATCCCCGCGATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGa                                                                                                                                           >  7:223807/1‑149 (MQ=255)
      cAAAGGAGCATCCCCGCGATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGAt                                                                                                                                          <  5:125926/149‑1 (MQ=255)
       aaaGGAGCATCCCCGCGATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCTGCCCCATGATCATTGAGCCAAGGAACATCGTGAt                                                                                                                                          >  2:231897/1‑148 (MQ=255)
                  cccGCGATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAAc                                                                                                                              <  8:223807/149‑1 (MQ=255)
                   ccGCGATACCTAATGAGAAGTTATTGATTACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACg                                                                                                                             <  8:173020/149‑1 (MQ=255)
                    cgcgATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGa                                                                                                                            <  8:216138/149‑1 (MQ=255)
                     gcgATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGAt                                                                                                                           <  8:308992/149‑1 (MQ=255)
                     gcgATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGAt                                                                                                                           >  4:224567/1‑149 (MQ=255)
                     gcgATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGAt                                                                                                                           <  1:195360/149‑1 (MQ=255)
                       gATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCCTGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATca                                                                                                                         <  3:82923/149‑1 (MQ=255)
                        aTACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATcac                                                                                                                        <  3:224567/149‑1 (MQ=255)
                            cTAATGAGAAGTTATTGATAACCATCTCAACACCGGCTGGTATGCGGTTTTCCAGTGCCTTATCGCCATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACAcc                                                                                                                     >  2:121621/1‑148 (MQ=255)
                              aaTGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCGTTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATc                                                                                                                          <  4:247379/141‑1 (MQ=255)
                              aaTGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCGTTGAGACAAGGAACATCGGGATCTCTGCGCCAACGATc                                                                                                                          >  3:247379/1‑141 (MQ=255)
                              aaTGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATc                                                                                                                          <  6:158871/141‑1 (MQ=255)
                              aaTGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATc                                                                                                                          >  5:158871/1‑141 (MQ=255)
                               aTGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCCACGATCACACCaa                                                                                                                   >  6:175858/1‑147 (MQ=255)
                               aTGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTATCCAGCGCCTTATCGACAAATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATa                                                                                                                 <  6:358826/149‑1 (MQ=255)
                                tGAGAAGTTATTGATTACAATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTTGGACCAACCCACCGAGCGGCCCCATCATCATTGAGCCAAGGAACATCGGGATCTCTTCGCCACCGCTCACACGAATAg                                                                                                                >  2:294590/1‑149 (MQ=255)
                                      gTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTa                                                                                                          <  1:195546/149‑1 (MQ=255)
                                        tATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTAtt                                                                                                        >  4:253786/1‑149 (MQ=255)
                                             aTAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGcc                                                                                                   >  7:275723/1‑149 (MQ=255)
                                              tAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGccc                                                                                                  <  2:143286/149‑1 (MQ=255)
                                                   aTCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGAc                                                                                             <  3:2374/149‑1 (MQ=255)
                                                                 aGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTAc                                                                               <  1:329243/149‑1 (MQ=255)
                                                                   gTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTAccg                                                                             <  4:102115/149‑1 (MQ=255)
                                                                       gCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCTGa                                                                         >  5:62596/1‑149 (MQ=255)
                                                                                   tGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTAc                                                                               <  1:144164/131‑1 (MQ=255)
                                                                                   tGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTAc                                                                               >  2:144164/1‑131 (MQ=255)
                                                                                     ccTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCt                                                           <  8:275723/149‑1 (MQ=255)
                                                                                     ccTTATCGACATATTTTATGACCAAAACCCCGCGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTCCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACAACCt                                                           >  2:315195/1‑149 (MQ=255)
                                                                                      cTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTg                                                          >  3:343141/1‑149 (MQ=255)
                                                                                             aCATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGac                                                                <  4:2839/136‑1 (MQ=255)
                                                                                             aCATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGac                                                                >  3:2839/1‑136 (MQ=255)
                                                                                              cATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAAcc                                                   >  1:275558/1‑148 (MQ=255)
                                                                                                tatTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCaa                                                 <  2:111720/148‑1 (MQ=255)
                                                                                                tatTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCaa                                                 >  1:111720/1‑148 (MQ=255)
                                                                                                         aCCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGGGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGAGCACCTGTAGAACCAATCATCACggg                                       <  5:274377/149‑1 (MQ=255)
                                                                                                         aCCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACggg                                       >  5:171629/1‑149 (MQ=255)
                                                                                                          ccAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACggg                                       >  5:304043/1‑148 (MQ=255)
                                                                                                             aCCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCaat                                   <  6:171629/149‑1 (MQ=255)
                                                                                                                 aCCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCTGACCAGATGACCACCTGTAGAACCAATCATCACGGGCaataaat                               <  6:62596/149‑1 (MQ=255)
                                                                                                                 aCCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCaataaat                               <  3:253786/149‑1 (MQ=255)
                                                                                                                 aCCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCaataaat                               <  3:281852/149‑1 (MQ=255)
                                                                                                                         gCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACTAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAatc                       <  6:304043/149‑1 (MQ=255)
                                                                                                                         gCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGACCAGCCTTTACCGCCCACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGAAat                        >  2:177263/1‑148 (MQ=255)
                                                                                                                            ccATGATCATTGAGCCAAGTAAAATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCTCGTTTACCTCCGCCCAGATGACCACCTGTACAACCAATCATAACGGGCAATAAATACGAAatcatc                    >  3:209436/1‑149 (MQ=255)
                                                                                                                            ccATGATCATTGAGCCAAGGAGCATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTAc                                                                               >  5:94294/1‑90 (MQ=255)
                                                                                                                            ccATGATCATTGAGCCAAGGAGCATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTAc                                                                               <  6:155091/90‑1 (MQ=255)
                                                                                                                            ccATGATCATTGAGCCAAGGAGCATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTAc                                                                               >  5:155091/1‑90 (MQ=255)
                                                                                                                            ccATGATCATTGAGCCAAGGAGCATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTAc                                                                               <  6:94294/90‑1 (MQ=255)
                                                                                                                            ccATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCGGACCAGATGACCACCTGTAGACCCAATCATCACGGTCAATAAATAGGTAatcatc                    >  4:87975/1‑149 (MQ=255)
                                                                                                                            ccATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAatcatc                    >  1:169517/1‑149 (MQ=255)
                                                                                                                             cATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCg                   >  2:245741/1‑149 (MQ=255)
                                                                                                                               tGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAAtcatca                                           >  2:179229/1‑123 (MQ=255)
                                                                                                                               tGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAAtcatca                                           <  1:179229/123‑1 (MQ=255)
                                                                                                                                 aTCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGcc               >  4:204516/1‑149 (MQ=255)
                                                                                                                                   cATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGa             >  6:294146/1‑149 (MQ=255)
                                                                                                                                       gAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGACAAt         >  4:1362/1‑149 (MQ=255)
                                                                                                                                        aGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGACAAtt        >  8:368778/1‑149 (MQ=255)
                                                                                                                                           cAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGACAATTTTg     <  2:275558/149‑1 (MQ=255)
                                                                                                                                             aGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGACAATTTTGGc   <  2:169517/149‑1 (MQ=255)
                                                                                                                                             aGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGACAATTTTGGc   <  5:294146/149‑1 (MQ=255)
                                                                                                                                              ggAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTaa                         <  6:90338/126‑1 (MQ=255)
                                                                                                                                              ggAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTaa                         >  5:90338/1‑126 (MQ=255)
                                                                                                                                               gAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGCCCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGCCAATTTTGGCg  >  7:353615/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AGCAGACAAAGGAGCATCCCCGCGATACCTAATGAGAAGTTATTGATAACCATCTCAAAACCGGCAGGTATGCGTTTTTCCAGTGCCTTATCGACATATTTTATGACCAACCCACCGAGCGGCCCCATGATCATTGAGCCAAGGAACATCGGGATCTCTGCGCCAACGATCACACCAATAGTACCTATTCCGCCCATGACCGCCCCGCGTTTACCGCCGACCAGATGACCACCTGTAGAACCAATCATCACGGGCAATAAATAGGTAATCATCGGGCCGACAATTTTGGCG  >  NC_000913/3078423‑3078713

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: