Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,839,621 G→A 6.4% R33H (CGC→CAC)  ascF → cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,839,6210GA 6.4% 193.7 / 2.4 62R33H (CGC→CAC) ascFcellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base A (2/2);  total (28/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 = 1.00e+00

GACTGATAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTG  >  NC_000913/2839477‑2839766
                                                                                                                                                |                                                                                                                                                 
gACTGATAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCttt                                                                                                                                               <  1:26460/149‑1 (MQ=255)
gACTGATAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCtt                                                                                                                                                <  1:306410/148‑1 (MQ=255)
 aCTGATAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTg                                                                                                                                              >  7:116212/1‑149 (MQ=255)
  cTGATAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTCGCGCTTTGt                                                                                                                                             <  3:280384/149‑1 (MQ=255)
   tGATAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGtt                                                                                                                                            <  4:29296/149‑1 (MQ=255)
     aTAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCTCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTTATAACATCTCGGCGGTCACGCACTTTATGACTCGGTTGCGCTTTGtt                                                                                                                                            <  4:88639/147‑1 (MQ=255)
      tAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTAt                                                                                                                                          >  6:138157/1‑148 (MQ=255)
       aacaacTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGGATGACGCGGTTGCGCTTTGTTATCa                                                                                                                                        >  8:332525/1‑148 (MQ=255)
         caacTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGCCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCttt                                                                                                                                               <  8:254942/140‑1 (MQ=255)
         caacTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGCCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCttt                                                                                                                                               >  7:254942/1‑140 (MQ=255)
          aacTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCa                                                                                                                                        >  1:225716/1‑146 (MQ=255)
             tACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAgat                                                                                                                                   >  2:352890/1‑148 (MQ=255)
               cATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAgatga                                                                                                                                 <  3:209427/148‑1 (MQ=255)
                   tACCCTACTGATAACAGTATAAAATCCGATGGCCAAAAATTATGCGTCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACTCACTGTATGACTCGGTTGCGCTTTGTTATCAAAGATGATGCaa                                                                                                                            >  4:215778/1‑148 (MQ=255)
                     cccTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACtt                                                                                                                          <  3:215778/149‑1 (MQ=255)
                       cTACTGATATCAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTAt                                                                                                                        <  1:352890/149‑1 (MQ=255)
                       cTACTGATAACAGGATAAAATCGGATGGACAAAAATTATGCGGCGATGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGATGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTAt                                                                                                                        <  5:222781/149‑1 (MQ=255)
                        tACTGATAACAGGATAAAATCCGATGGCCATAAATTATGCGGCGCTGGCACGCTCGGAGATAGCGGCACTGGGCGGCGTTGATAACATCACGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATAAAAGCTGATGCACTTATc                                                                                                                       >  6:329282/1‑149 (MQ=255)
                          cTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGa                                                                                                                     <  7:302729/149‑1 (MQ=255)
                             aTAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACa                                                                                                                   <  5:138157/148‑1 (MQ=255)
                                aCAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGccc                                                                                                               >  4:166757/1‑149 (MQ=255)
                                aCAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGccc                                                                                                               >  8:42196/1‑149 (MQ=255)
                                  aGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGa                                                                                                             >  6:32226/1‑149 (MQ=255)
                                   ggATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGAAAGAGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGTGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGAc                                                                                                            >  3:339858/1‑149 (MQ=255)
                                      tAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGatgatg                                                                                                                               <  6:211904/127‑1 (MQ=255)
                                      tAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGatgatg                                                                                                                               >  5:211904/1‑127 (MQ=255)
                                       aaaaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAgat                                                                                                                                   >  6:310214/1‑122 (MQ=255)
                                       aaaaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAgat                                                                                                                                   <  5:310214/122‑1 (MQ=255)
                                       aaaaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACa                                                                                                                   <  3:112477/138‑1 (MQ=255)
                                       aaaaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACa                                                                                                                   >  4:112477/1‑138 (MQ=255)
                                         aaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTgcgc                                                                                                                                                  <  6:311966/105‑1 (MQ=255)
                                         aaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTgcgc                                                                                                                                                  >  5:311966/1‑105 (MQ=255)
                                         aaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACa                                                                                                                   >  2:287511/1‑136 (MQ=255)
                                         aaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACa                                                                                                                   <  1:287511/136‑1 (MQ=255)
                                         aaTCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTaa                                                                                                       >  1:287233/1‑148 (MQ=255)
                                             cGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAAc                                                                                                   >  2:214588/1‑148 (MQ=255)
                                              gATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCa                                                                                                 >  4:122100/1‑149 (MQ=255)
                                               aTGGCCAAATATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCAt                                                                                                <  7:42196/149‑1 (MQ=255)
                                               aTGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCAt                                                                                                <  3:166757/149‑1 (MQ=255)
                                                tGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCttt                                                                                                                                               <  6:32001/101‑1 (MQ=255)
                                                tGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCttt                                                                                                                                               >  5:32001/1‑101 (MQ=255)
                                                tGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATc                                                                                               <  3:122100/149‑1 (MQ=255)
                                                  gcataaaaTTATGCGGCGCTGGCACGCTCGGTGAAAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGTGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATccc                                                                                             <  4:339858/145‑1 (MQ=255)
                                                   ccAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCa                                                                                                                             >  8:341369/1‑116 (MQ=255)
                                                   ccAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCa                                                                                                                             <  7:341369/116‑1 (MQ=255)
                                                     aaaaaTTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAg                                                                                                                                     <  1:94726/106‑1 (MQ=255)
                                                     aaaaaTTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAg                                                                                                                                     >  2:94726/1‑106 (MQ=255)
                                                     aaaaaTTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCACTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCgg                                                                                          >  8:360229/1‑149 (MQ=255)
                                                     aaaaaTTATGCGGCGCTGGCACGCTCGGTGATAGCGGCAAGGGGCGGCGTTGATAACATCTCGGCGGTCCCGCACTGTATGACGCGGTTGCGCTTTGTTAACAAATATGCTGCACTTATCGCAAGCCCGACGTTAAAACCCATCCCCgg                                                                                          >  5:67621/1‑149 (MQ=255)
                                                      aaaaTTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCaa                                                                                                                                       <  7:75705/103‑1 (MQ=255)
                                                      aaaaTTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCaa                                                                                                                                       >  8:75705/1‑103 (MQ=255)
                                                                ggcgCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACACACTGTATGACGCGGTTGCACTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCg                                                                               >  1:54387/1‑149 (MQ=255)
                                                                 gcgcTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGt                                                                              <  5:28098/149‑1 (MQ=255)
                                                                               ggTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGa                                                                  <  2:225716/147‑1 (MQ=255)
                                                                                      gCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAgt                                                          <  5:32226/148‑1 (MQ=255)
                                                                                       cGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCACTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAgtgt                                                        <  7:360229/149‑1 (MQ=255)
                                                                                         tcACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTTCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAgtgt                                                        >  6:103755/2‑147 (MQ=255)
                                                                                          cACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAg                                                     <  2:287233/149‑1 (MQ=255)
                                                                                           aCTGGGCGGCGTTGATAACATCTCGGCGGTCACACACTGTATGACGCGGTTGCACTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAgg                                                    <  2:54387/149‑1 (MQ=255)
                                                                                                           aaCATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTgg                                                                            <  6:101046/109‑1 (MQ=255)
                                                                                                                 tcGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGttt                               >  1:169414/1‑148 (MQ=255)
                                                                                                                   ggcggTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTcac                            <  1:79718/149‑1 (MQ=255)
                                                                                                                             cACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCa                  <  7:332525/149‑1 (MQ=255)
                                                                                                                                    tgccgcgGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCAGGAAgt            <  8:116212/145‑1 (MQ=255)
                                                                                                                                              gcgcTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATGCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTAGCGCCAATACCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCctg  <  5:103755/148‑1 (MQ=255)
                                                                                                                                                gcTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTgg                                                                            <  6:260644/72‑1 (MQ=255)
                                                                                                                                                gcTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCAGGAAGTCGTCAGc     <  4:198697/143‑1 (MQ=255)
                                                                                                                                                gcTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCAGGAAGTCGTCAGc     >  3:198697/1‑143 (MQ=255)
                                                                                                                                                gcTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCACGGCGTGCTCGGCGTgg                                                                            >  5:260644/1‑72 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GACTGATAACAACTACATCTACCCTACTGATAACAGGATAAAATCCGATGGCCAAAAATTATGCGGCGCTGGCACGCTCGGTGATAGCGGCACTGGGCGGCGTTGATAACATCTCGGCGGTCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTG  >  NC_000913/2839477‑2839766

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: