Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,415,793 G→A 6.4% A11V (GCA→GTA) ‡ gabP‑V ← gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,415,7930GA 6.4% 219.5 / 4.8 78A11V (GCA→GTA) ‡gabP‑Vgamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base G (31/42);  new base A (2/3);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

GTACCACCAGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTGTTGTTGATATTTGGCACG  >  NC_002947/5415647‑5415942
                                                                                                                                                   |                                                                                                                                                     
gTACCACCAGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCaa                                                                                                                                                    >  2:2123698/1‑150 (MQ=255)
gTACCACCAGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCaa                                                                                                                                                    <  2:25587/150‑1 (MQ=255)
 tACCACCAGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAg                                                                                                                                                   <  2:680457/150‑1 (MQ=255)
  accaccAGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAgg                                                                                                                                                  <  1:472799/150‑1 (MQ=255)
    caccaGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGACACCTTGTGCAAGGtt                                                                                                                                                >  2:828097/1‑150 (MQ=255)
       caGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGcc                                                                                                                                             <  2:527086/150‑1 (MQ=255)
        aGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGCACAAGGtttcaa                                                                                                                                            <  1:1719847/150‑5 (MQ=255)
          ggTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATc                                                                                                                                          <  1:941209/150‑1 (MQ=255)
           gTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCa                                                                                                                                         <  2:1751091/150‑1 (MQ=255)
           gTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCa                                                                                                                                         <  2:198499/150‑1 (MQ=255)
             aCCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGa                                                                                                                                       <  2:2146109/150‑1 (MQ=255)
                      gAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGgt                                                                                                                              >  1:358558/1‑150 (MQ=255)
                       aaaTAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGgtg                                                                                                                             <  1:2305408/150‑1 (MQ=255)
                          tAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATTGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTc                                                                                                                          >  2:1278617/1‑150 (MQ=255)
                            ggAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCAt                                                                                                                        >  2:80624/1‑150 (MQ=255)
                             gAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGCTGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCACGGTTCCCATCAGATCCCCCCGTGTTCAtt                                                                                                                       >  2:915763/1‑150 (MQ=255)
                                aTTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAAACGTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGAt                                                                                                                    >  1:1591347/1‑150 (MQ=255)
                                 ttATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAGCCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGAtg                                                                                                                   <  2:2179862/150‑1 (MQ=255)
                                 ttATGGTTGCTGGGCCGGATGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACTTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGAtg                                                                                                                   <  2:532354/150‑1 (MQ=255)
                                      gTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGAAGCCTAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATgttgtt                                                                                                              <  1:511959/150‑1 (MQ=255)
                                             ggCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATc                                                                                                                                          <  1:1019706/115‑1 (MQ=255)
                                             ggCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATc                                                                                                                                          >  2:1019706/1‑115 (MQ=255)
                                             ggCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCtttt                                                                                                       >  1:1129627/1‑150 (MQ=255)
                                               ccGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCAtt                                                                                                                       <  2:527035/132‑1 (MQ=255)
                                               ccGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCAtt                                                                                                                       >  1:527035/1‑132 (MQ=255)
                                               ccGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTg                                                                                                     >  1:1531735/1‑150 (MQ=255)
                                                cGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGc                                                                                                    <  1:985575/150‑1 (MQ=255)
                                                 ggCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCt                                                                                                   <  2:1805741/150‑1 (MQ=255)
                                                   cTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGa                                                                                                 <  1:80624/150‑1 (MQ=255)
                                                   cTGCCGCGATGGCATGGCCAGAGCCGACCAAAAGCCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGa                                                                                                 <  1:2302859/150‑1 (MQ=255)
                                                       cgcgATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCttttt                                                                                                      >  1:454310/1‑141 (MQ=255)
                                                       cgcgATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCttttt                                                                                                      <  2:454310/141‑1 (MQ=255)
                                                        gcgATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCttttt                                                                                                      <  1:1471425/140‑1 (MQ=255)
                                                        gcgATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCttttt                                                                                                      >  2:1471425/1‑140 (MQ=255)
                                                        gcgATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACaa                                                                                            <  2:732628/150‑1 (MQ=255)
                                                        gcgATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCATTGTTGGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACaa                                                                                            <  1:1205136/150‑1 (MQ=255)
                                                            tggcatggcCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCATGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACaaggaa                                                                                        <  1:420450/150‑1 (MQ=255)
                                                            gggcatggcCAGAGCCGACGAAAAGTACGGCGCCGCTGCTGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACaaggaa                                                                                        <  2:1591347/149‑1 (MQ=255)
                                                                atggcCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGAc                                                                                    >  2:1679046/1‑150 (MQ=255)
                                                                atggcCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGAc                                                                                    <  1:1609051/150‑1 (MQ=255)
                                                                atggcCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGCACAAGGtttca                                                                                                                                             <  1:1022303/93‑4 (MQ=39)
                                                                atggcCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGCACAAGGtttca                                                                                                                                             >  2:1022303/1‑90 (MQ=255)
                                                                     cAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAAGCTTGTGCAAGGTTGCCATCAGATCGCACGGTGTTCATTGATGTTGTTCTCTTTTCAATGCTACACGGAACGACGGGCa                                                                               >  2:733694/1‑150 (MQ=255)
                                                                      agagCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAGACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAg                                                                              <  1:828097/150‑1 (MQ=255)
                                                                      agagCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTTCTGACACAAGGAAGGACGGGCAg                                                                              <  2:1697585/150‑1 (MQ=255)
                                                                       gagCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAgg                                                                             <  1:473477/150‑1 (MQ=255)
                                                                            gacgaAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAg                                                                        <  2:295163/150‑1 (MQ=255)
                                                                             acgaAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAgg                                                                       >  1:2128609/1‑150 (MQ=255)
                                                                               gaAAAGTCCGGCGCCGATGATGCCGGCGATTGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAGGTACGGGCAGGCGTAGGCg                                                                     >  1:120794/1‑150 (MQ=255)
                                                                               gaAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACATTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGa                                                                                                 <  1:560333/122‑1 (MQ=255)
                                                                               gaAAAGTACGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACATTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGa                                                                                                 >  2:560333/1‑122 (MQ=255)
                                                                                    gTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGaaaaa                                                                >  2:2023184/1‑150 (MQ=255)
                                                                                         gcgcCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCTCACGGTGTTCCTTGATGTTGGTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATAtct                                                           <  1:2123698/150‑1 (MQ=255)
                                                                                           gcCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTcc                                                         <  1:856857/150‑1 (MQ=255)
                                                                                               atgatgCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTg                                                     <  2:1475221/150‑1 (MQ=255)
                                                                                                  atgCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGaaa                                                  <  2:1531735/150‑1 (MQ=255)
                                                                                                        ggCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAAt                                             <  1:1888959/150‑1 (MQ=255)
                                                                                                        ggCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAAt                                             >  2:364196/1‑150 (MQ=255)
                                                                                                              ggACAGCATCGTTACGTGACGTGGCGTGAAACCTGGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCg                                       >  2:1969933/1‑150 (MQ=255)
                                                                                                               gACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGAtg                                                                                                                   >  2:2409911/1‑73 (MQ=255)
                                                                                                               gACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGAtg                                                                                                                   <  1:2409911/73‑1 (MQ=255)
                                                                                                               gACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGa                                      <  2:120794/150‑1 (MQ=255)
                                                                                                                  aGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGacac                                   >  1:667076/1‑150 (MQ=255)
                                                                                                                    cATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAat                                                               <  1:390256/120‑1 (MQ=255)
                                                                                                                    cATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAat                                                               >  2:390256/1‑120 (MQ=255)
                                                                                                                     aTCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACgtt                                >  2:698649/1‑150 (MQ=255)
                                                                                                                         ttACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTc                            <  2:1129627/150‑1 (MQ=255)
                                                                                                                                    tgatggcaaCCTTGCACAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCtgttg                 >  2:425358/8‑150 (MQ=255)
                                                                                                                                      cTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTgttgtt               <  1:2119893/150‑1 (MQ=255)
                                                                                                                                      atggcaaCCTTGCACAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTgttgtt               <  1:425358/145‑1 (MQ=255)
                                                                                                                                       ttGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTGttgttg              >  1:221540/1‑150 (MQ=255)
                                                                                                                                       ttGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTGttgttg              >  2:1701213/1‑150 (MQ=255)
                                                                                                                                          aaaCCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTGTTGTTGata           >  2:327771/1‑150 (MQ=255)
                                                                                                                                           aaCCTTGTGCAAGGTTGCCATAAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGaca                                    >  1:2286153/1‑124 (MQ=255)
                                                                                                                                           aaCCTTGTGCAAGGTTGCCATAAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGaca                                    <  2:2286153/124‑1 (MQ=255)
                                                                                                                                            aCCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTgttgtt               <  2:2132694/144‑1 (MQ=255)
                                                                                                                                            aCCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTgttgtt               >  1:2132694/1‑144 (MQ=255)
                                                                                                                                             ccTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTGTTGTTGATAttt        >  2:1859786/1‑150 (MQ=255)
                                                                                                                                                tgtgCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTACTAACAGAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCCCGCTGTTGTTGATATTTGGc     <  1:1207198/150‑1 (MQ=255)
                                                                                                                                                   gCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTTTTTTTTTTTTTAATGATGGGCGGGCTTAAGCGAAAAATATATCCTTTGAAATAAATAATTCTACGCGTTGTTCTCGCGCTGTTTTTGATATTTGGcgcc  >  2:162981/1‑147 (MQ=255)
                                                                                                                                                   gCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCttt                                                      <  2:1146861/98‑1 (MQ=255)
                                                                                                                                                   gCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCttt                                                      >  1:1146861/1‑98 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTACCACCAGGGTACCGGCCACGAAATAGGAAATTATGGTTGCTGGGCCGGCTGCCGCGATGGCATGGCCAGAGCCGACGAAAAGTCCGGCGCCGATGATGCCGGCGATGGACAGCATCGTTACGTGACGTGGCTTGAAACCTTGTGCAAGGTTGCCATCAGATCCCACGGTGTTCATTGATGTTGTTCTCTTTTTGCTGACACAAGGAAGGACGGGCAGGCGTAGGCGAAAAATATCTCCTTTGAAATCAATAAGTCGACACGTTGTTCGCGCGCTGTTGTTGATATTTGGCACG  >  NC_002947/5415647‑5415942

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