Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,654,962 G→A 5.0% T217T (ACG→ACA cysB → transcriptional regulator CysB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,654,9620GA 5.0% 247.7 / 2.0 79T217T (ACG→ACAcysBtranscriptional regulator CysB
Reads supporting (aligned to +/- strand):  ref base G (36/39);  new base A (2/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

TTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGTGAGCTGTTTGAA  >  NC_002947/2654820‑2655109
                                                                                                                                              |                                                                                                                                                    
agcgggtgCCCCTGCGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAACCTTTCAATAACCGTGGCCTCACGCCaa                                                                                                                                                  >  8:369115/5‑147 (MQ=255)
   tggtgCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGt                                                                                                                                              <  5:78074/148‑1 (MQ=255)
   tggtgCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTg                                                                                                                                             <  7:288070/149‑1 (MQ=255)
     gtgCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCATCACGCAAAAAGTGGt                                                                                                                                           >  8:306735/1‑149 (MQ=255)
          gCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCa                                                                                                                                      >  4:317297/1‑149 (MQ=255)
          gCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCa                                                                                                                                      >  4:169279/1‑149 (MQ=255)
                    ccGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGa                                                                                                                           >  3:63705/1‑149 (MQ=255)
                       cTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGc                                                                                                                        <  4:396487/149‑1 (MQ=255)
                       cTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGc                                                                                                                        >  3:396487/1‑149 (MQ=255)
                          gCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGa                                                                                                                     >  4:47215/1‑149 (MQ=255)
                           ccAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGAc                                                                                                                    <  7:134430/149‑1 (MQ=255)
                           ccAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGAc                                                                                                                    <  4:345990/149‑1 (MQ=255)
                            cAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAACTCGTCTTCACCGCGGCTGACGCCGACg                                                                                                                   >  6:74791/1‑149 (MQ=255)
                                   ccGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCaa                                                                                                            <  1:200138/149‑1 (MQ=255)
                                        aCTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACg                                                                                                        >  2:329125/1‑148 (MQ=255)
                                        aCTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGt                                                                                                       <  5:399097/149‑1 (MQ=255)
                                         cTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTa                                                                                                      <  7:281536/149‑1 (MQ=255)
                                            aCCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCCAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTAtgt                                                                                                   <  1:294406/149‑1 (MQ=255)
                                               cTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTgcg                                                                                                >  3:207590/1‑149 (MQ=255)
                                                tGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTgcgc                                                                                               >  7:34361/1‑149 (MQ=255)
                                                tGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTgcgc                                                                                               >  3:54798/1‑149 (MQ=255)
                                                  gAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCt                                                                                             <  3:317297/149‑1 (MQ=255)
                                                    acccGTCGCCGAGTTACCGATACTCAACTACGTGTTCGGGTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTgg                                                                                           <  6:160518/147‑1 (MQ=255)
                                                    aGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTgg                                                                                           >  7:357139/1‑149 (MQ=255)
                                                       cGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGt                                                                                        <  6:152875/149‑1 (MQ=255)
                                                             cGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACACCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCg                                                                                  >  6:130110/1‑149 (MQ=255)
                                                             cGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACACCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCg                                                                                  >  6:124303/1‑149 (MQ=255)
                                                              gAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGt                                                                                 <  1:397914/149‑1 (MQ=255)
                                                                gTACCCGATCGTCACCTTCGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCATAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTgg                                                                               <  5:74791/149‑1 (MQ=255)
                                                                gTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTgg                                                                               <  6:221722/149‑1 (MQ=255)
                                                                gTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTgg                                                                               <  3:107896/149‑1 (MQ=255)
                                                                    ccGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTAt                                                                           >  6:334097/1‑149 (MQ=255)
                                                                    ccGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCATTCAATCACCGTGGCCTCACGCCAAAAGTCGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACCTATCTGCGCCTGGGGTTGGGCGTGTGTAt                                                                           >  6:286582/1‑149 (MQ=255)
                                                                     cGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATc                                                                          <  3:47215/149‑1 (MQ=255)
                                                                     cGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATc                                                                          >  4:94745/1‑149 (MQ=255)
                                                                          gTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGc                                                                     <  8:34361/149‑1 (MQ=255)
                                                                          gTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGc                                                                     <  3:169279/149‑1 (MQ=255)
                                                                           tCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGcc                                                                    <  1:170418/149‑1 (MQ=255)
                                                                                 agttgtTCGGTTTCACCGGGCGCTCAAAGCTGGATCAACCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTggg                                                                                    <  6:304291/124‑1 (MQ=255)
                                                                                 aCGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTggg                                                                                    >  5:304291/1‑127 (MQ=255)
                                                                                     gtTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGTGCCTGGGGTTGGGCGTgg                                                                               >  7:343097/1‑128 (MQ=255)
                                                                                     gtTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGTGCCTGGGGTTGGGCGTgg                                                                               <  8:343097/128‑1 (MQ=255)
                                                                                       tCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGc                                                                 <  3:170295/140‑1 (MQ=255)
                                                                                       tCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGc                                                                 >  4:170295/1‑140 (MQ=255)
                                                                                          gTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGAc                                                     <  5:334097/149‑1 (MQ=255)
                                                                                           tttCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTCACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCg                                                                   >  2:343856/1‑134 (MQ=255)
                                                                                             tCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGc                                                   >  6:75217/1‑148 (MQ=255)
                                                                                              cACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGcc                                                  >  2:213594/1‑148 (MQ=255)
                                                                                               aCCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGcc                                                                    <  7:139037/129‑1 (MQ=255)
                                                                                               aCCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGcc                                                                    >  8:139037/1‑129 (MQ=255)
                                                                                               aCCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGcc                                                  >  6:239608/1‑147 (MQ=255)
                                                                                               aCCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCaa                                                <  3:94745/149‑1 (MQ=255)
                                                                                                cctttctcTCAATGCTGTATGTAGGTTTCAATCAGCGTGGCCTCAGGGCAAATGTGGTCTTCACCGGGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCaa                                                <  8:155761/141‑1 (MQ=255)
                                                                                                ccGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAg                                               >  2:222143/1‑149 (MQ=255)
                                                                                                ccGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAg                                               >  7:45661/1‑149 (MQ=255)
                                                                                                       cTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACa                                        <  3:193037/149‑1 (MQ=255)
                                                                                                        tCAAAGCTGGATGAAGCTTTCAATCACCTTGGCCTCACGCCAAAAGTGGTCTTCACCGCCGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTCGAGATCGCCGTGGGTATCGTAGCCGGCATGCCTCTGGACGCCaa                                                >  7:155761/1‑140 (MQ=255)
                                                                                                             gCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGCCGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACt                                  >  5:292767/1‑149 (MQ=255)
                                                                                                                 gATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTtgg                               <  8:45661/148‑1 (MQ=255)
                                                                                                                 gATGAAGCTTTCAATCACCGTGGCCTCACACCAAAAGTGGTCTTCACCGCGGCTGACTCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTtg                                <  5:130110/147‑1 (MQ=255)
                                                                                                                 gATGAAGCTTTCAATCACCGTGGCCTCACACCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTtg                                <  5:124303/147‑1 (MQ=255)
                                                                                                                     aaGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGgc                           <  3:347182/148‑1 (MQ=255)
                                                                                                                       gCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCt                        >  1:268996/1‑149 (MQ=255)
                                                                                                                         tttCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGATGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTg                       >  3:228802/1‑148 (MQ=255)
                                                                                                                         tttCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTtggtg                             <  1:208128/142‑1 (MQ=255)
                                                                                                                         tttCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTtggtg                             >  2:208128/1‑142 (MQ=255)
                                                                                                                         tttCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTg                       >  5:127177/1‑148 (MQ=255)
                                                                                                                             aaTCACCTTGGCCTGACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCTAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGc                  <  5:286582/149‑1 (MQ=255)
                                                                                                                             aaTCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCaa                                                <  5:164255/119‑1 (MQ=255)
                                                                                                                             aaTCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCaa                                                >  6:164255/1‑119 (MQ=255)
                                                                                                                             aaTCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATg                   <  4:63705/148‑1 (MQ=255)
                                                                                                                                  ccGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGt              >  4:38519/1‑148 (MQ=255)
                                                                                                                                   cGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCGTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGt              >  1:60039/1‑148 (MQ=255)
                                                                                                                                     tGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGTGAGc          <  8:357139/149‑1 (MQ=255)
                                                                                                                                      ggCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGTGAGCt         >  8:276532/1‑149 (MQ=255)
                                                                                                                                         cTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGTGAGCTGtt      <  7:276532/149‑1 (MQ=255)
                                                                                                                                            aCGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGTGAGCTGttt     <  1:329125/147‑1 (MQ=255)
                                                                                                                                             cGCCAAAAGTGGTCTTCACCGCGTCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCTGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGTGAGCTGTTTGaa  <  1:403610/149‑1 (MQ=255)
                                                                                                                                             cGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCTCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGTCGCTGGATGCCAGTGAGCTGTTTGaa  <  7:369115/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TTGTGGTGCCGCAGGGGCACCCGCTGGCCAAGCTGCCGAAACTGACCCTGGAAGCCGTCGCCGAGTACCCGATCGTCACCTACGTGTTCGGTTTCACCGGGCGCTCAAAGCTGGATGAAGCTTTCAATCACCGTGGCCTCACGCCAAAAGTGGTCTTCACCGCGGCTGACGCCGACGTGATCAAGACGTATGTGCGCCTGGGGTTGGGCGTGGGTATCGTTGCCGGCATGGCTGTGGACGCCAAGCTCGACAGTGACTTGGTGGCGCTGGATGCCAGTGAGCTGTTTGAA  >  NC_002947/2654820‑2655109

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: