Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,239,8230TG44.1% 64.3 / 13.6 62G315G (GGT→GGGdadXalanine racemase, catabolic, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base T (11/22);  new base G (25/1);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 4.45e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAG  >  NC_000913/1239673‑1239971
                                                                                                                                                      |                                                                                                                                                    
cGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGt                                                                                                                                                      <  1:638937/151‑1 (MQ=255)
   cATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTAt                                                                                                                                                    <  2:465078/150‑1 (MQ=255)
       gTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGCATTggga                                                                                                                                               >  2:314653/1‑149 (MQ=255)
                gggTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGGTg                                                                                                                                      >  1:1026574/1‑151 (MQ=255)
                gggTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGGCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGGAAGGCGGTTg                                                                                                                                      >  2:930164/1‑151 (MQ=255)
                  gTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGGTGAg                                                                                                                                    >  1:1137014/1‑151 (MQ=255)
                   tACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAg                                                                                                                                    <  2:23499/150‑1 (MQ=255)
                     cGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTAAGcccggtttggtt                                                                                                                                 >  2:14924/1‑141 (MQ=255)
                     cGCCGACGGTTATCCGCGCCACGCGCCGACCGGTACGCCGGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTGGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAACtg                                                                                                                                 <  2:176291/151‑1 (MQ=255)
                           cGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCCGTTGGGCTgggggg                                                                                                                            >  2:897771/1‑150 (MQ=255)
                             gTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGGCCGTTGAGCTGTGGGGCaa                                                                                                                         >  2:655396/1‑151 (MQ=255)
                               tATCCGGGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAgg                                                                                                                       <  2:2458/151‑1 (MQ=255)
                                         ccgcgCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAAt                                                                                                             <  1:898149/150‑1 (MQ=255)
                                         aCGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACGGTCTCGATGGATATGCTAGCGGTCGATTTAGCGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGAGCAAAAt                                                                                                             <  2:425619/151‑1 (MQ=255)
                                         aCGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAAt                                                                                                             <  2:1054941/151‑1 (MQ=255)
                                            cgcCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATtga                                                                                                          <  2:435560/151‑1 (MQ=255)
                                            cgcCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTGAGCTGTGGGGCAAGGGGATCAAAATtga                                                                                                          >  1:181329/1‑151 (MQ=255)
                                                tACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGGTGAGCTGTGGGGCAAGGGGATCAAAATttatgtg                                                                                                      >  2:159194/1‑149 (MQ=255)
                                                 aCCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGatgatg                                                                                                     <  1:762724/151‑1 (MQ=255)
                                                          ccTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGTTTGGCGGGGGGGCCCGGAGATCACAATTGGTGATTTcgccccc                                                                                             >  1:738035/1‑150 (MQ=255)
                                                           cTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGTTGGGCTGTGGGGCAAGGGGATCAAAAATGATGATGTcgccgccc                                                                                            >  2:456921/1‑149 (MQ=255)
                                                           cTGTTTTAGTGGACGGCGTGCGCACCATGACGGAGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTACCGCCTTGCCCGCCGGCGGGGGTTGGGGCCCCGGCTGAGCTGTGGGGCCAGGAGATACAAAATGTAGCTCTcgccgccg                                                                                            >  2:515531/1‑150 (MQ=255)
                                                             gTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGGACGCCGGGTGAGCTGTGGGGCAAGGAGGGCAAAATTAATGATGTcgccgccc                                                                                            >  1:873517/1‑147 (MQ=255)
                                                             gTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTTGGCTGTGGGGGAAGGGGATCAAAATTTATTGTTTcgcccccg                                                                                            >  1:940339/1‑148 (MQ=255)
                                                             gTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTGAGCTGTGGGGCAAGGAGATCAAAAATGATGATGTcgccgccg                                                                                            >  2:349534/1‑148 (MQ=255)
                                                                ttAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTTGTACGCCGGTTGGGGTGTGGGGGAAGGGGATCAAAATTTATGATTTCGCCGCCGCTGCCg                                                                                      >  1:287542/1‑151 (MQ=255)
                                                                    tGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGaa                                                                                   <  2:1158430/150‑1 (MQ=255)
                                                                       aCGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACgg                                                                                <  2:1136488/150‑1 (MQ=255)
                                                                          gCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTTAGGTGTGGGGGAAGGGGGTCAAAAATGGTTATGTCGCCGCCGCGGCCGGGACgggtgg                                                                            >  2:633547/1‑147 (MQ=255)
                                                                          gCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGTTGGGCTGTGGGGCCAGGAGATCAAAATTGTTGATGTCGCCCCCGCCTCCGGGACgggggg                                                                            >  2:424838/1‑151 (MQ=255)
                                                                                 cACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAg                                                                     <  1:455953/151‑1 (MQ=255)
                                                                                 cACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGGACGGTGGGGGATGAg                                                                     >  2:536752/1‑151 (MQ=255)
                                                                                 cACCAGGACGTTGGGGACCAACTGGATGGATATGCTAGCGTTCGATTGAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATAAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCAATGa                                                                      <  2:143109/150‑1 (MQ=255)
                                                                                  aCCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGt                                                                    <  1:633861/151‑1 (MQ=255)
                                                                                  aCCAAGACGTGGGGGCCCGTCTATATGGATATGCGGGCGGTTTGTTGACCGCCGTGCCCACGGGCGGGGATTGGGACGCCGGTTGGGCTGTGGGGCAGGGACAGGAAAATGGATGTTGTCGCCGCCGCTGCCGGAACTGTGGGCTATGAGt                                                                    <  2:607798/151‑1 (MQ=255)
                                                                                   ccATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTTCGCCGGTTGAGCTGGGGGGCAAGGGGATCAAAATTTGTGATGTCGCCGCCGCCGCCCGGACGGTGGGGTATGGGtt                                                                   >  2:497459/1‑151 (MQ=255)
                                                                                      tGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGGGGGCTATGAGTTGAt                                                                >  1:1046467/1‑151 (MQ=255)
                                                                                         cGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGGACGCCGGGTGTGCTGTGGGGCAAGGGGATCCAAATTGATGATGTCGCCGCCGCTGCCGGGACGGGGGGCTATGAGTTGGtgtg                                                             >  2:78341/1‑151 (MQ=255)
                                                                                         cGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGGTGGGCCGTGGGGCCAGGAGGTCAAAAATTGTGCTGTCGCCGCCCCTGCCGGGACGGTGGGCCATGAGGTGAtgtg                                                             >  2:1053864/1‑151 (MQ=255)
                                                                                          ggTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTGAGCTGTGGGGCAAGGAGCTCAAAATTGATGATGTCGCCTCCGGTGCCGGGACGGGGGGGTTTGAGTTTAtgtgg                                                            >  2:623206/1‑150 (MQ=255)
                                                                                                    gTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGAGCAAAATTGATGATGTCGCCCCCGCTGCCGGAACGGTGGGCTATGAGGTGAGGTGGGCGCTGgccg                                                  >  2:796889/1‑149 (MQ=255)
                                                                                                    gTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGACTTGGTACGCCGGTTTAGCTTTGGGGGAAGGAGATCAAAATTTATGATGTCTCCGCCGCTCCCGGGACGGTGGGCTATGAGTTTATGTGCGCGCCGgcg                                                   >  1:706342/1‑150 (MQ=255)
                                                                                                     tctcGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGgcgc                                                  <  2:1152837/150‑1 (MQ=255)
                                                                                                     tctaGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGgcgc                                                  <  2:563301/150‑1 (MQ=255)
                                                                                                         gATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGTGGGGCAAGGGGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGGGGGCTATGAGTTGATTTGCGCGCTGGGGCTAcg                                              >  2:863074/1‑150 (MQ=255)
                                                                                                            ggATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTAcgc                                             <  2:288956/148‑1 (MQ=255)
                                                                                                            ggATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGGACGCCGGTTGAGCTGTGGGGGAAGGGGATCAAAATTGATGATGTCGCCGCCCCCGCCGGAACGGGGGGGTTTGAGTTGATGTGCGCGCTGGCGCCAcgcccg                                          >  1:289045/1‑148 (MQ=255)
                                                                                                            ggATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTGGGCTGTGGGGGAAGGGGATCAAAATTGATGATGTTGCCGCCGCTGCCGGAACGGGGGGGTTTGAGTTGATGTGCGCGCTGGCGCTTCGCGt                                           >  1:1105576/1‑150 (MQ=255)
                                                                                                            ggATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGTTGAGCTGTGGGGGAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGGGGGCTATGAGTTGATGTGGGCGCTGGCGCTAcgcgg                                           >  1:566838/1‑149 (MQ=255)
                                                                                                              atatGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTTGTACGCCGGGTGGGCTGTGGGGCCAGGGGATCAAAATTGATGATGTCGCCGGCGCTGCCGGGACGGGGGGGGATGGGTTGATGTGCGCGCTGGCGCTACGCGttc                                         >  2:680984/1‑148 (MQ=255)
                                                                                                                    tAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGCTTGCTACGCCGGGTTAAACGTTGGGCTAGGAGTCCCAAATTGATGATGTTGCCACCCCCGCCCGAACGGTGGGCTTCTAATTTATATGCGCGATGGCGCTAccccgctccgtcgg                                  >  2:455453/1‑138 (MQ=255)
                                                                                                                             aTTTAACGCCTTGCCCGCAGGCGGGTATTTGTACGCCGGGTGAGCTGTGGGGCAAGGAGGTCAAAATTGATGAGGTCGCCGCCGCTGCCGGGACGGGGGGCTATGAGGTGATGTGGGCGCTGGCGCTACGCGTCCCGGGTGTGACGGGGTa                         >  2:690123/1‑151 (MQ=255)
                                                                                                                             aTTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGGTGAGCTGTGGGGCAAGGAGATCCAAATTGATGATGTCGCCGCCGCTGCCGGGACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGGGTa                         >  2:650307/1‑151 (MQ=255)
                                                                                                                             aTTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGGTGATGTCGCCGCCGCTGCCGGAAAGGTGGGGTATGAGGTGATGTGCGGGCGGGCGCCACGCGTCCCGGGTGGGAAGGGGTa                         >  1:758523/1‑151 (MQ=255)
                                                                                                                             aTTTAACGCCTTGCCCGCAGGCGGGCATTGGTACGCCGGTTGAGGTGTGGGGGAAGGAGAGGAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGCGGGGTATGAGTTGATCTGCGCGCCGGCGGTACGCGTCCCCGGTGTGACGGGGTa                         >  2:1158296/1‑151 (MQ=255)
                                                                                                                                 aaCGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACtt                     <  2:1170157/151‑1 (MQ=255)
                                                                                                                                 aaCGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACtt                     <  2:181278/151‑1 (MQ=255)
                                                                                                                                 aaCGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACtt                     <  2:873143/151‑1 (MQ=255)
                                                                                                                                   cGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACttgt                   <  1:65068/151‑1 (MQ=255)
                                                                                                                                     ccTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACttgttg                 <  2:118655/151‑1 (MQ=255)
                                                                                                                                                    ggTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAg  >  2:682099/1‑151 (MQ=255)
                                                                                                                                                    ggTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAg  >  1:971701/1‑151 (MQ=255)
                                                                                                                                                    ggTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAg  >  1:224859/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAG  >  NC_000913/1239673‑1239971

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