Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,068,945 T→A 6.5% I36N (ATT→AAT)  aroP‑I → aromatic amino acid transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,068,9450TA 6.5% 233.0 / 4.8 77I36N (ATT→AAT) aroP‑Iaromatic amino acid transport protein
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base A (2/3);  total (38/39)
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.18e-01

AGCAACAACACTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGAAGAGCCGGTCGCCGGCAGCTTCAGCCACTTCGCCCAC  >  NC_002947/1068802‑1069087
                                                                                                                                               |                                                                                                                                              
aGCAACAACACTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGcc                                                                                                                                            >  7:172789/1‑148 (MQ=255)
 gCAACAACACTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCATAGCACGCCCGAGCGCTTACAGCGCGCGCTGAAGAATCTCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTCTTACTGGGCATTGCCCa                                                                                                                                          >  7:120179/1‑149 (MQ=255)
   aacaacACTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGa                                                                                                                                        <  5:2330/149‑1 (MQ=255)
    acaacaCTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAg                                                                                                                                         <  6:334244/147‑1 (MQ=255)
    acaacaCTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAg                                                                                                                                         >  5:334244/1‑147 (MQ=255)
     caacaCTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCAATGCCCAggaa                                                                                                                                      <  2:127138/149‑4 (MQ=255)
      aacaCTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCa                                                                                                                                     >  1:105440/1‑149 (MQ=255)
            gCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGt                                                                                                                               >  8:298993/1‑149 (MQ=255)
             ccTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAg                                                                                                                                <  3:180069/147‑1 (MQ=255)
                  cTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTggccg                                                                                                                         >  6:227521/1‑149 (MQ=255)
                   tGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTggccgg                                                                                                                        >  7:102211/1‑149 (MQ=255)
                       cTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCACACCAACAAGTTGGCCGGcccc                                                                                                                    >  1:239168/1‑149 (MQ=255)
                        tAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCGAGACCATGAAGTTGGACGTCCACt                                                                                                                   >  8:69960/1‑149 (MQ=255)
                         aaaGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGcc                                                                                                                                            <  2:51300/123‑1 (MQ=255)
                         aaaGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGcc                                                                                                                                            >  1:51300/1‑123 (MQ=255)
                         aaaGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTc                                                                                                                  >  8:76087/1‑149 (MQ=255)
                             gACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCgtg                                                                                                              >  2:357166/1‑149 (MQ=255)
                             gACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCgtg                                                                                                              >  3:128779/1‑149 (MQ=255)
                             gACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCgtg                                                                                                              >  6:278264/1‑149 (MQ=255)
                                acTTCATGCAAGACCAGAGCACGCGCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTg                                                                                                           <  8:172789/149‑1 (MQ=255)
                                acTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTg                                                                                                           >  4:3653/1‑149 (MQ=255)
                                acTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTg                                                                                                           >  4:179151/1‑149 (MQ=255)
                                acTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTg                                                                                                           >  4:26564/1‑149 (MQ=255)
                                   tCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCAATGCCCAgg                                                                                                                                        >  7:98422/1‑116 (MQ=255)
                                   tCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCAATGCCCAgg                                                                                                                                        <  8:98422/117‑2 (MQ=255)
                                    cATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTg                                                                                                        <  2:105440/148‑1 (MQ=255)
                                     aTGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTggg                                                                                                      >  5:300322/1‑149 (MQ=255)
                                      tGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGc                                                                                                     <  5:278264/149‑1 (MQ=255)
                                        cAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTa                                                                                                   >  8:310021/1‑149 (MQ=255)
                                         aaGACCAGAGCACGCCCGAGCGGATACAGCGCGGGCCGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTAc                                                                                                  >  7:53997/1‑149 (MQ=255)
                                          aGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACg                                                                                                 <  7:69960/149‑1 (MQ=255)
                                              cAGAGCACGCCCGAGCGGTTACAGCGCTTGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCTGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTGCGTGTTGCTTGGCTACGCCAt                                                                                             <  8:203105/149‑1 (MQ=255)
                                              cAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCAt                                                                                             >  6:138695/1‑149 (MQ=255)
                                                gagCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCg                                                                                           >  4:237064/1‑149 (MQ=255)
                                                   tcagCCCGCGCTGTAACCGCTCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCg                                                                                        >  8:330786/4‑149 (MQ=255)
                                                     cGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGc                                                                                      <  7:298993/149‑1 (MQ=255)
                                                      gCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGcc                                                                                     >  6:202087/1‑149 (MQ=255)
                                                      gCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGcc                                                                                     <  3:179151/149‑1 (MQ=255)
                                                           aGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATg                                                                                <  1:357166/149‑1 (MQ=255)
                                                           aGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATg                                                                                <  2:179610/149‑1 (MQ=255)
                                                             cGGTTACAGGGCGGGCCGAAGAATCGCCATTTCCAGCTGATCGCGGTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGTCCAGACCTTCCAGTTGGCTGGCTCTTCCGTGTTGCTGGGCTACGTTATCGTCGGTCTGATGGc                                                                              <  8:53997/149‑1 (MQ=255)
                                                             cGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGc                                                                              <  8:102211/149‑1 (MQ=255)
                                                                ttACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCtt                                                                           >  8:41750/1‑149 (MQ=255)
                                                                ttACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCtt                                                                           <  5:138695/149‑1 (MQ=255)
                                                                  accgctcgGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTc                                                                          <  7:330786/142‑1 (MQ=255)
                                                                  aCAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTc                                                                          <  3:26564/148‑1 (MQ=255)
                                                                   cAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTcc                                                                         <  3:237064/148‑1 (MQ=255)
                                                                    aGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCt                                                                        >  6:112458/1‑148 (MQ=255)
                                                                      cgcgGGCTGAATAATCGCCATATCCATCTGCTCGCGCTGGGCGGGGCCATCGGTACCGGGGTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGAt                                                                     <  2:239168/149‑1 (MQ=255)
                                                                           gCTGAAGACTCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATgc                                                                <  5:227521/149‑1 (MQ=255)
                                                                           gCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATgc                                                                <  7:76087/149‑1 (MQ=255)
                                                                             tGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATgcgc                                                              >  5:227032/1‑149 (MQ=255)
                                                                              gaagaaTCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTg                                                                                                        <  8:41052/106‑1 (MQ=255)
                                                                              gaagaaTCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTg                                                                                                        >  7:41052/1‑106 (MQ=255)
                                                                               aagaaTCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATgcgc                                                              >  7:246703/1‑147 (MQ=255)
                                                                                agaaTCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAg                                                           <  7:319685/149‑1 (MQ=255)
                                                                                agaaTCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAg                                                           <  7:324335/149‑1 (MQ=255)
                                                                                    tCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTgg                                                       <  6:227032/149‑1 (MQ=255)
                                                                                      gCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGt                                                     <  1:211469/149‑1 (MQ=255)
                                                                                       ccATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGt                                                     >  8:98221/1‑148 (MQ=255)
                                                                                        cATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTg                                                    >  3:45954/1‑148 (MQ=255)
                                                                                          tatCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGaaa                                                 <  7:310021/149‑1 (MQ=255)
                                                                                          tatCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGaaa                                                 <  4:45954/149‑1 (MQ=255)
                                                                                           atCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGGTGGGTGAAAt                                                >  6:254052/1‑149 (MQ=255)
                                                                                              cAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAAtggt                                             >  5:358302/1‑149 (MQ=255)
                                                                                               aGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAAtggtg                                            <  6:300322/149‑1 (MQ=255)
                                                                                                gCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCCGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCGGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAAtggtgg                                           <  5:366445/149‑1 (MQ=255)
                                                                                                 cTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAAtggtggt                                          <  6:358302/149‑1 (MQ=255)
                                                                                                     tCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGaa                                      <  4:128779/149‑1 (MQ=255)
                                                                                                                      aTCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGAAGAGCCGGTCGCCGGCAg                     >  8:49932/1‑149 (MQ=255)
                                                                                                                          gTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGAAGAGCCGGTCGCCGGCAGCTTc                 <  8:246703/149‑1 (MQ=255)
                                                                                                                           tACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGAAGAGCCGGTCGCCGGCAGCTTCa                <  5:202087/149‑1 (MQ=255)
                                                                                                                                    gttCCTGGGCATTGCCCAGACCATCCTGTTGGCCGGCCCCTGCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGAAGAGCCGGTCGCCGGCAGCTTCAGCCACTTc        <  7:98221/148‑1 (MQ=255)
                                                                                                                                       tcTGGGCAATGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAAtggtggtgg                                        >  2:256403/2‑113 (MQ=255)
                                                                                                                                       tcTGGGCAATGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAAtggtggtgg                                        <  1:256403/112‑1 (MQ=255)
                                                                                                                                         tGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGAAGAGCCGGTCGCCGGCAGCTTCAGCCACTTCGCccac  >  6:55429/1‑149 (MQ=255)
                                                                                                                                          gggCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGAAGAGCCGGTCGCCGGCAGCTTCAGCCACTTCGCccac  >  3:32467/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                              
AGCAACAACACTGCCTCACTGAGCTAAAGGACACTTCATGCAAGACCAGAGCACGCCCGAGCGGTTACAGCGCGGGCTGAAGAATCGCCATATCCAGCTGATCGCGCTGGGCGGGGCCATCGGTACCGGGTTGTTCCTGGGCATTGCCCAGACCATCCAGTTGGCCGGCCCCTCCGTGTTGCTGGGCTACGCCATCGCCGGCCTGATGGCCTTCCTGATCATGCGCCAGCTGGGTGAAATGGTGGTGGAAGAGCCGGTCGCCGGCAGCTTCAGCCACTTCGCCCAC  >  NC_002947/1068802‑1069087

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: ‑