Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,078,184 C→T 9.2% G225D (GGC→GAC)  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,1840CT 9.2% 137.6 / 5.4 56G225D (GGC→GAC) cmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/3);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

CATGGCAATAATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGG  >  NC_000913/3078036‑3078332
                                                                                                                                                    |                                                                                                                                                    
cATGGCAATAATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTATCGAACAAGTTAAACGCCAGCAGCAGCCCCAGGc                                                                                                                                                      <  1:794453/149‑1 (MQ=255)
 aTGGCAATAATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAATGTAAACGCCAGCAGCAGCCCCAGGcc                                                                                                                                                     >  2:610610/1‑149 (MQ=255)
    gCAATACGGGGCAGCGGCGTCACCAGCACATACGGGAAATAGAGTTCGGGGATCCCACCGAGGAAGTGAAGAATCATCGCCCCGGGCGCAGAACGTATACACACCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCgg                                                                                                                                                  <  1:70420/149‑1 (MQ=255)
          atGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTgg                                                                                                                                            <  1:368167/149‑1 (MQ=255)
          atGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTCCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTgg                                                                                                                                            <  2:323213/149‑1 (MQ=255)
            ggTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTTGTGGCTCCCACCGAGGAAGTGAATAATCCTACCCCCGTGCGCAGCAGGTCGTCTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTgggg                                                                                                                                          <  1:189958/149‑2 (MQ=255)
             gTCAGCGGCTTCATCAGCATATACGGGAAACACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggc                                                                                                                                         <  2:663084/149‑3 (MQ=255)
              tCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggct                                                                                                                                        <  1:345007/149‑4 (MQ=255)
              tCAGCGGCTCCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACTGAGGAAGCAAATAATCATTCCCCCGGGCACACTACGCGTAGCGCTCCCTTTACCAACGAAGGTAGACGCCAGCCGCAGCCCCAGGCCCGGCCCGGGGCTa                                                                                                                                        >  2:326072/1‑149 (MQ=255)
                aGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTaga                                                                                                                                      <  2:788957/149‑1 (MQ=255)
                      ttCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccac                                                                                                                                <  1:742586/149‑1 (MQ=255)
                       tcatcaGCACATACGGAAAATGCAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacc                                                                                                                               <  1:491408/149‑1 (MQ=255)
                       tcagccgCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacc                                                                                                                               <  1:319867/143‑1 (MQ=255)
                        catcaGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCACCAGGCCCGGACCTGGGTTAGAGGccacca                                                                                                                              <  1:418572/149‑1 (MQ=255)
                           caGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaa                                                                                                                           <  1:739801/149‑1 (MQ=255)
                             gCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaag                                                                                                                         <  1:691570/149‑1 (MQ=255)
                                caTACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGAAAACGCCAGCAGCAGCCCCAGGCCCGGGCCTGGGTTAGAGGCCCCCATAaagaag                                                                                                                      >  1:344135/1‑149 (MQ=255)
                                      ggAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCAACGCCCCGGGACCAGAACGGTTACTCATCCCCTCACCAAACAAGGGAAACGCCAGCAGCAACCCCAGGCCCGCACCTGGGTTAAAGGCCACCAAAAACAACATGGAt                                                                                                                >  2:349862/1‑149 (MQ=255)
                                       gAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGAtt                                                                                                               >  1:128493/1‑149 (MQ=255)
                                         aaTACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTa                                                                                                             <  1:534036/149‑1 (MQ=255)
                                          aTACAGTTCGGGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGTCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTAc                                                                                                            <  1:389154/149‑1 (MQ=255)
                                              aGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGtt                                                                                                        <  1:316108/149‑1 (MQ=255)
                                               gTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTa                                                                                                       >  2:677194/1‑149 (MQ=255)
                                               gTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAAAGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCACCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTa                                                                                                       >  2:177445/1‑149 (MQ=255)
                                                ttCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                                                                      <  1:230797/149‑1 (MQ=255)
                                                   gTGGATCCCACCGAGGAAGTGAATAATCATAGCCCCGAGCGCATCCAGTTTACTCACACGTGTACCAAACAAGGTAAACGCAAGCATCACTCTCACGCTCGATCTTGGCTTAGAGGCCATCTTTAAGAGGTTGCGTTTACCGTTCAcca                                                                                                   >  2:790387/1‑148 (MQ=255)
                                                    tGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGa                                                                                                  <  1:70598/149‑1 (MQ=255)
                                                     ggATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGaa                                                                                                 <  1:7680/149‑1 (MQ=255)
                                                      gATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGACAGCAGCAGCCCCAGGCCCGGACGTGGGTGAGAGGCCACCATACAGAAGACGGATTTACCGATAACCGAAg                                                                                                >  2:645243/1‑149 (MQ=255)
                                                            aCCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCCGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGtt                                                                                          <  1:614853/149‑1 (MQ=255)
                                                             ccGAGGAAGTGCATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGGGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTg                                                                                         <  2:610978/149‑1 (MQ=255)
                                                                    aGTGACCAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTccc                                                                                  >  2:531675/1‑149 (MQ=255)
                                                                      tGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTGAGACTCATGGCACACTGCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGGCTGTTGCATTCCCAg                                                                                >  1:790229/1‑149 (MQ=255)
                                                                      tGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCCGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAg                                                                                >  1:464863/1‑149 (MQ=255)
                                                                       gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGTGGGATTTCCCTTTAACCGAAGCCTGTTGCATTCCAAGc                                                                               >  1:150120/1‑149 (MQ=255)
                                                                       gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGc                                                                               >  1:197756/1‑149 (MQ=255)
                                                                       gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGc                                                                               >  2:394028/1‑149 (MQ=255)
                                                                               atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataata                                                                       <  1:121177/149‑1 (MQ=255)
                                                                                tcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCCTTCCCAGCGGATAATAg                                                                      >  2:79828/1‑149 (MQ=255)
                                                                                               gaCGTTTCCTCAGCCCTTTACCAAACAAGGTAAACGCCAGCAGCACCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTcatccctc                                                       <  2:450301/148‑5 (MQ=255)
                                                                                                 cGTTTACTCATCCCTTTACCAAACAACGTAAACAAAAATTGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGc                                                     >  1:792953/1‑149 (MQ=255)
                                                                                                         cctCCCTTTACCAAACAAGGTCAACGCCAGCAGCAGCCCCAGGCCCCGCCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGa                                             <  2:153157/147‑1 (MQ=255)
                                                                                                         cATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGa                                             <  1:73120/149‑1 (MQ=255)
                                                                                                          aTCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAg                                            >  1:464359/1‑149 (MQ=255)
                                                                                                              cTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGGTCGATCGCCTTATTGAGGaaa                                        >  2:683015/1‑149 (MQ=255)
                                                                                                                  aCCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATGCAGGAAAAGCa                                    >  2:538607/1‑149 (MQ=255)
                                                                                                                    cAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGGTAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCAcg                                  >  1:713652/1‑148 (MQ=255)
                                                                                                                    cAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACt                                  >  1:799745/1‑149 (MQ=255)
                                                                                                                     aaaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACtt                                 >  2:735821/1‑149 (MQ=255)
                                                                                                                      aaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACttt                                >  2:42141/1‑149 (MQ=255)
                                                                                                                       aCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAGGCACTTTc                               <  1:88853/149‑1 (MQ=255)
                                                                                                                        cAAGGTAAACGCTATCAGCATCCCCAGGCCCGGACGTGGGTTAGATGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCg                              <  1:219634/149‑1 (MQ=255)
                                                                                                                                    cagcagcagCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATtga                  >  1:72451/1‑149 (MQ=255)
                                                                                                                                        ctctaacccAGGTCCGGGCCTGGGGTAGAGGCCAGCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATGGCATGCTAGAGGAAAAGCACTTTCGCCGGCTCATGTAtgac              >  1:403272/7‑149 (MQ=255)
                                                                                                                                           agcCCCAGGACCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgactga           <  1:499815/149‑1 (MQ=255)
                                                                                                                                             acccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGaca         >  1:739735/2‑149 (MQ=255)
                                                                                                                                               ccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGacacc       >  1:273569/1‑147 (MQ=255)
                                                                                                                                                cAGGACCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGacaaca      <  2:581729/149‑1 (MQ=255)
                                                                                                                                                    cccGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATgg  <  1:685745/149‑1 (MQ=255)
                                                                                                                                                    cccGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTATCCCAGCGTTATTGGTGACGGacaaccaggg  >  1:742957/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGGCAATAATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGG  >  NC_000913/3078036‑3078332

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑