Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,899,253 T→G 7.9% I72I (ATA→ATC PP_4307 ← DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,899,2530TG 7.9% 191.8 / 5.2 63I72I (ATA→ATCPP_4307DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base G (3/2);  total (34/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 = 9.76e-01

GGTGATCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGC  >  NC_002947/4899106‑4899401
                                                                                                                                                    |                                                                                                                                                     
ggTGATCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTa                                                                                                                                                      >  8:167827/1‑149 (MQ=255)
     tCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGc                                                                                                                                                  <  4:128059/148‑1 (MQ=255)
       ttttCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAg                                                                                                                                               >  4:102196/1‑149 (MQ=255)
       ttttCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAg                                                                                                                                               >  7:25652/1‑149 (MQ=255)
             cccAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATc                                                                                                                                                    >  5:196934/1‑138 (MQ=255)
             cccAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATc                                                                                                                                                    <  6:196934/138‑1 (MQ=255)
                     cAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGc                                                                                                                                 >  8:14131/1‑149 (MQ=255)
                       gcgcCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCgg                                                                                                                               >  1:62043/1‑149 (MQ=255)
                         gccgccAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCt                                                                                                                             >  6:384829/1‑149 (MQ=255)
                           cgccAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGa                                                                                                                           <  6:128676/149‑1 (MQ=255)
                            gccAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGAt                                                                                                                          >  7:144385/1‑149 (MQ=255)
                                  cccACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCAcccc                                                                                                                    >  5:11813/1‑149 (MQ=255)
                                   ccaccaGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCAccccc                                                                                                                   <  5:212747/149‑1 (MQ=255)
                                     accaGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCgg                                                                                                                 <  7:330569/149‑1 (MQ=255)
                                       caGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGgc                                                                                                               >  3:243234/1‑149 (MQ=255)
                                         gCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGgcgc                                                                                                             <  7:167827/149‑1 (MQ=255)
                                          cATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCa                                                                                                            <  4:243234/149‑1 (MQ=255)
                                                  gCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGc                                                                                                                                 >  7:109442/1‑120 (MQ=255)
                                                  gCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGc                                                                                                                                 <  8:109442/120‑1 (MQ=255)
                                                   cTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCa                                                                                                   >  3:99211/1‑148 (MQ=255)
                                                       aCGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAAcgc                                                                                              >  2:347264/1‑149 (MQ=255)
                                                        cGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAAcgcg                                                                                             <  8:25652/149‑1 (MQ=255)
                                                        cGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAAcgcg                                                                                             >  1:196513/1‑149 (MQ=255)
                                                            gTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACgcggcg                                                                                          >  8:53712/1‑148 (MQ=255)
                                                              acacCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCa                                                                                        >  8:160129/1‑148 (MQ=255)
                                                               cacCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCaaa                                                                                      <  5:384829/149‑1 (MQ=255)
                                                                   gccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGc                                                                                  <  7:14131/149‑1 (MQ=255)
                                                                   gccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGc                                                                                  <  6:67304/149‑1 (MQ=255)
                                                                   gccgTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGc                                                                                  <  4:99211/149‑1 (MQ=255)
                                                                          cGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAg                                                                           <  1:347264/149‑1 (MQ=255)
                                                                            aTGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGt                                                                         >  1:89906/1‑149 (MQ=255)
                                                                              ggCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATc                                                                                    <  1:128433/136‑1 (MQ=255)
                                                                              ggCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATc                                                                                    >  2:128433/1‑136 (MQ=255)
                                                                               gCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCAc                                                                      >  5:204536/1‑149 (MQ=255)
                                                                                cAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCgg                                                                                                                 >  1:253568/1‑106 (MQ=255)
                                                                                cAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCgg                                                                                                                 <  2:253568/106‑1 (MQ=255)
                                                                                     gACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCAc                                                                 <  8:144385/148‑1 (MQ=255)
                                                                                         cccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                           >  4:231425/1‑149 (MQ=255)
                                                                                          ccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                           <  7:301082/149‑1 (MQ=255)
                                                                                          ccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                           <  2:196513/149‑1 (MQ=255)
                                                                                          ccGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATa                                                           >  4:69198/1‑149 (MQ=255)
                                                                                                gTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgc                                                     >  2:170611/1‑149 (MQ=255)
                                                                                                gTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgc                                                     <  3:102196/149‑1 (MQ=255)
                                                                                                 tATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgcg                                                    <  7:160129/149‑1 (MQ=255)
                                                                                                  aTGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAgcgc                                                   >  3:98328/1‑149 (MQ=255)
                                                                                                     aaGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCa                                                 <  7:235514/148‑1 (MQ=255)
                                                                                                       gCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGaa                                              <  2:62043/149‑1 (MQ=255)
                                                                                                               cagcagCCAGTTGCCGACAATGGCCACCCCGCATGGGTTTCGCCTCCAAAATTCCCGTCCGGGTGTACACCCCCGGCCCACAAAGCCCAACCCCGGCGAAAAACGCCGGAAAGGAACGGGCCCCGATTACGCGCGGCCAGacaccccc                                        >  2:350273/1‑141 (MQ=255)
                                                                                                               cagcagCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACAt                                      >  6:375039/1‑149 (MQ=255)
                                                                                                                      cAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCAc                               <  4:98328/149‑1 (MQ=255)
                                                                                                                        gTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCGAGGTCCCGGCACCGAATACGCAGCGCCAGAACACCACAGCCAGCACCg                             >  6:369992/1‑149 (MQ=255)
                                                                                                                           gCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCt                          >  5:111502/1‑149 (MQ=255)
                                                                                                                                 aaTGGCCACCCCGCTGGCTATCGCCAGCAAAAATGACGTGCGGCTCATCACACCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGCTCACGGCACCGAAAACCCAGCGCCAGAACAACACATCCAGAACCCGCTGGccc                     >  4:332398/1‑147 (MQ=255)
                                                                                                                                   tGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGac                  <  5:369992/149‑1 (MQ=255)
                                                                                                                                     gCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGacac                <  3:69198/149‑1 (MQ=255)
                                                                                                                                       atttttgCTGGCGATAGCCAGCGAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGcc                      <  6:288947/135‑1 (MQ=255)
                                                                                                                                       atttttgCTGGCGATAGCCAGCGAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGcc                      >  5:288947/7‑141 (MQ=255)
                                                                                                                                         ccccGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCa           <  3:231425/149‑1 (MQ=255)
                                                                                                                                            tgCTGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGc                                                        <  3:74250/101‑1 (MQ=255)
                                                                                                                                            tgCTGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGc                                                        >  4:74250/2‑102 (MQ=255)
                                                                                                                                            cGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCAc          >  5:187287/1‑148 (MQ=255)
                                                                                                                                             gCTGGCGATAGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACg         >  3:331646/1‑148 (MQ=255)
                                                                                                                                                tGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAAc      >  8:296673/1‑149 (MQ=255)
                                                                                                                                                    tATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGc  >  5:230726/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTGATCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGC  >  NC_002947/4899106‑4899401

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: