Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 589,213 G→C 6.3% intergenic (‑568/‑664) PP_0498 ← / → PP_0500 hypothetical protein/dTDP‑4‑rhamnose reductase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947589,2130GC 6.3% 190.9 / 2.6 63intergenic (‑568/‑664)PP_0498/PP_0500hypothetical protein/dTDP‑4‑rhamnose reductase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (25/34);  new base C (2/2);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

GTAGGCGGCGATCAGGATCAGCAGTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAGGCCACCGCCCCCGGCAATAGCATCAAAAAAGCCGG  >  NC_002947/589065‑589359
                                                                                                                                                    |                                                                                                                                                  
gTAGGCGGCGATCAGGATCAGCAGTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCg                                                                                                                                                    <  8:68395/149‑1 (MQ=255)
gTAGGCGGCGATCAGGATCAGCAGTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCg                                                                                                                                                    <  2:30912/149‑1 (MQ=255)
 tAGGCGGCGATCAGGATCAGCAGTGGAGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCgg                                                                                                                                                   <  8:298745/149‑1 (MQ=255)
         gatcaggatcagCAGTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGgc                                                                                                                                           <  5:274268/149‑1 (MQ=255)
             aggatcagCAGTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTCGCCCGCAAAGgcatga                                                                                                                                       <  5:3007/149‑5 (MQ=255)
             aggatcagCAGTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTCGCCCGCAAAGgcatga                                                                                                                                       >  6:3007/1‑145 (MQ=255)
                   agcagTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGc                                                                                                                                  <  7:65827/148‑1 (MQ=255)
                    gcagTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAg                                                                                                                                             <  1:314377/136‑1 (MQ=255)
                    gcagTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAg                                                                                                                                             >  2:314377/1‑136 (MQ=255)
                    gcagTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGcc                                                                                                                                 >  7:33311/1‑148 (MQ=255)
                       gtggtgGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGCGCGAAGGCGCACGTAGCCGaaa                                                                                                                             <  7:203451/149‑1 (MQ=255)
                          gtggGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCg                                                                                                                          >  2:247484/1‑149 (MQ=255)
                               acgcacgcCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTCGCCCGCAAAGgca                                                                                                                                          >  5:39430/1‑127 (MQ=255)
                               acgcacgcCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTCGCCCGCAAAGgca                                                                                                                                          <  6:39430/128‑2 (MQ=255)
                                 gcacgcCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAAc                                                                                                                            <  3:302642/140‑1 (MQ=255)
                                 gcacgcCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAAc                                                                                                                            >  4:302642/1‑140 (MQ=255)
                                   acgcCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAg                                                                                                                 <  8:31228/149‑1 (MQ=255)
                                     gcCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGACATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCt                                                                                                               >  4:193255/1‑149 (MQ=255)
                                          cAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCt                                                                                                           <  7:192852/148‑1 (MQ=255)
                                           aaaaCGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTg                                                                                                          >  4:5115/1‑148 (MQ=255)
                                                gTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACtt                                                                                                    >  7:61790/1‑149 (MQ=255)
                                                gTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACtt                                                                                                    <  8:61790/149‑1 (MQ=255)
                                                   tttGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGt                                                                                                  <  1:247484/148‑1 (MQ=255)
                                                       gggACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTg                                                                                             >  3:10651/1‑149 (MQ=255)
                                                        ggACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTgg                                                                                            <  1:22932/149‑1 (MQ=255)
                                                          aCAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCg                                                                                          <  4:2057/149‑1 (MQ=255)
                                                           cAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGa                                                                                         <  3:5115/149‑1 (MQ=255)
                                                            aaaGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGa                                                                                         >  3:254087/1‑148 (MQ=255)
                                                            aaaGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGa                                                                                         >  4:169699/1‑148 (MQ=255)
                                                              aGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCTCGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAg                                                                                      >  5:93271/1‑149 (MQ=255)
                                                                 tATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTg                                                                                   >  6:166522/1‑149 (MQ=255)
                                                                   tCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGc                                                                                  <  5:166522/148‑1 (MQ=255)
                                                                     gACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGTACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGcc                                                                                 >  8:72175/1‑147 (MQ=255)
                                                                     gACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGcc                                                                                 >  3:165934/1‑147 (MQ=255)
                                                                     gACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCaa                                                                               <  7:21000/149‑1 (MQ=255)
                                                                        aaGGCCCCCAGCGCCCAACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAggg                                                                            <  7:72175/149‑1 (MQ=255)
                                                                          ggCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGtt                                                                          <  3:169699/149‑1 (MQ=255)
                                                                              cccAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTac                                                                      <  8:216192/149‑1 (MQ=255)
                                                                                cAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTacac                                                                    >  8:1486/1‑149 (MQ=255)
                                                                                    gcCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGTAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGc                                                                <  3:193255/149‑1 (MQ=255)
                                                                                     ccccACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCa                                                               <  4:165934/149‑1 (MQ=255)
                                                                                         aCCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATa                                                            >  2:106751/1‑148 (MQ=255)
                                                                                          ccGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATaaa                                                          <  7:1486/149‑1 (MQ=255)
                                                                                           cGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACAcc                                                                   >  8:233942/1‑139 (MQ=255)
                                                                                           cGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACAcc                                                                   <  7:233942/139‑1 (MQ=255)
                                                                                           cGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATaaa                                                          >  5:70425/1‑148 (MQ=255)
                                                                                           cGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAaac                                                         <  6:70425/149‑1 (MQ=255)
                                                                                            gATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAaaca                                                        <  4:10651/149‑1 (MQ=255)
                                                                                                     cATGACCGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGa                                               <  7:134809/149‑1 (MQ=255)
                                                                                                       tGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTg                                                                                             <  1:271766/101‑1 (MQ=255)
                                                                                                       tGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTg                                                                                             >  2:271766/1‑101 (MQ=255)
                                                                                                           aGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATaaa                                                          <  7:87118/132‑1 (MQ=255)
                                                                                                              cATCGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTATCCGAAGCTGCCTGGAACTTGTTCGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATaa                                      <  4:4816/149‑1 (MQ=255)
                                                                                                                tGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAgg                                    >  1:252767/1‑149 (MQ=255)
                                                                                                                   gACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCt                                                                                    >  1:288804/1‑98 (MQ=255)
                                                                                                                   gACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCt                                                                                    <  2:288804/98‑1 (MQ=255)
                                                                                                                   gACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAGGCCa                                 >  3:96706/1‑149 (MQ=255)
                                                                                                                    aCTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAGGCCAc                                >  2:274975/1‑149 (MQ=255)
                                                                                                                                  aGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTATCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAACGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAGGCCACCGCCCCCGGCAATa                  <  4:237204/149‑1 (MQ=255)
                                                                                                                                  aGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAggg                                                                            <  1:41174/91‑1 (MQ=255)
                                                                                                                                  aGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAggg                                                                            >  2:41174/1‑91 (MQ=255)
                                                                                                                                  aGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAGGCCACCGCCCCCGGCAATa                  <  3:221527/149‑1 (MQ=255)
                                                                                                                                   gTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGc                                                                >  4:165735/1‑102 (MQ=255)
                                                                                                                                   gTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGc                                                                <  3:165735/102‑1 (MQ=255)
                                                                                                                                                 tGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAGGCCACCGCCCCCGGCAATAGCATCAAAAAAGc     <  1:274975/147‑1 (MQ=255)
                                                                                                                                                   cGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAGGCCACCGCCCCCGGCAATAGCATCAAAAAAGCCgg  <  8:33311/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTAGGCGGCGATCAGGATCAGCAGTGGTGGGACGCACGCCTGCAAAACGTTTTTGGGGACAAAGTTATCGACAAGGCCCCCAGCGCCCCACCGATAAACGCCATGACAGCCATGGGACTGCCTCGTCTCCAGTCGATCATGCCTTTGCGGGCGAAGGCGCACGTAGCCGAAACCGTAGCCGAAGCTGCCTGGAACTTGTTAGTGGCGATAGCTGCCAAGGGTTCTACACCGGCAATAAACAACACCGGGAGCGTGATAAGGCCACCGCCCCCGGCAATAGCATCAAAAAAGCCGG  >  NC_002947/589065‑589359

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: