Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,961,628 T→A 5.1% I175N (ATT→AAT)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,961,6280TA 5.1% 310.3 / 4.3 97I175N (ATT→AAT) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (53/39);  new base A (2/3);  total (55/42)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

ACCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGC  >  NC_002947/5961483‑5961773
                                                                                                                                                  |                                                                                                                                                 
aCCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgc                                                                                                                                                >  1:226517/1‑149 (MQ=255)
aCCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgc                                                                                                                                                <  5:402044/149‑1 (MQ=255)
  cgacgaGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGGCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCGTTAAAGAAAACAAGTTCGGCATCGCCATtgccg                                                                                                                                              <  6:440212/149‑1 (MQ=255)
   gacgaGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccga                                                                                                                                             >  7:383646/1‑149 (MQ=255)
        gCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccgatgcc                                                                                                                                         >  5:13838/1‑148 (MQ=255)
         cTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGGCAAtggcgatgccga                                                                                                                                       <  2:213776/149‑1 (MQ=255)
         cTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGCACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACGAGTTCGGCATCGCCATtgccgatgcc                                                                                                                                         >  5:308055/1‑147 (MQ=255)
            gAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAg                                                                                                                                      <  8:331716/147‑1 (MQ=255)
            gAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAg                                                                                                                                      <  6:13838/147‑1 (MQ=255)
             aGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGcc                                                                                                                                   <  8:118537/149‑1 (MQ=255)
                 ccTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATc                                                                                                                                >  8:62451/1‑148 (MQ=255)
                 ccTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATc                                                                                                                                >  1:389059/1‑148 (MQ=255)
                        gtagtCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtg                                                                                                                                                 <  1:234959/124‑1 (MQ=255)
                        gtagtCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtg                                                                                                                                                 >  2:234959/1‑124 (MQ=255)
                           gtCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACgcgc                                                                                                                     >  1:134248/1‑149 (MQ=255)
                            tCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCCCCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgcc                                                                                                                    >  3:338997/1‑149 (MQ=255)
                            tCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgcc                                                                                                                    >  6:147898/1‑149 (MQ=255)
                              gccgccGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgccg                                                                                                                   >  3:202941/1‑148 (MQ=255)
                               ccgccgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGca                                                                                                                 <  3:412493/149‑1 (MQ=255)
                                 gccgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGCCGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGca                                                                                                                 >  1:142996/1‑147 (MQ=255)
                                  ccgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAg                                                                                                              >  1:23283/1‑149 (MQ=255)
                                  ccgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAg                                                                                                              >  1:21337/1‑149 (MQ=255)
                                   cgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGc                                                                                                             <  5:363778/149‑1 (MQ=255)
                                      aaTGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAc                                                                                                          <  5:321068/149‑1 (MQ=255)
                                               ggTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTgg                                                                                                 <  7:62451/149‑1 (MQ=255)
                                                 ttGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                               <  5:330761/149‑1 (MQ=255)
                                                    ccccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTg                                                                                            <  3:46237/149‑1 (MQ=255)
                                                     cccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTgg                                                                                           >  4:306608/1‑149 (MQ=255)
                                                       cGGTGTCGCTGCGGGTCAACCCGGACGTCGATTCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACGAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTc                                                                                         <  6:308055/149‑1 (MQ=255)
                                                       cGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTc                                                                                         <  5:70365/149‑1 (MQ=255)
                                                       cGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTc                                                                                         <  2:185924/149‑1 (MQ=255)
                                                        ggTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                        <  8:166204/149‑1 (MQ=255)
                                                        ggTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                        <  2:134248/149‑1 (MQ=255)
                                                              gCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTc                                                                                   >  6:120764/1‑148 (MQ=255)
                                                                  cGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACt                                                                               <  5:224200/148‑1 (MQ=255)
                                                                   gggTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGc                                                                             <  8:383646/149‑1 (MQ=255)
                                                                    ggTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGcc                                                                            <  2:389059/149‑1 (MQ=255)
                                                                      tCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCg                                                                                  >  7:248367/1‑141 (MQ=255)
                                                                      tCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCg                                                                                  <  8:248367/141‑1 (MQ=255)
                                                                        aaCCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACAt                                                                        <  5:147898/149‑1 (MQ=255)
                                                                         aCCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATc                                                                       >  5:179080/1‑149 (MQ=255)
                                                                           ccGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCgg                                                                     >  6:237661/1‑149 (MQ=255)
                                                                             ggACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCt                                                                   <  2:226517/149‑1 (MQ=255)
                                                                                  tCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAg                                                              <  2:21337/149‑1 (MQ=255)
                                                                                      tGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGa                                                          <  3:352942/149‑1 (MQ=255)
                                                                                                cacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTc                                                  >  3:194182/1‑147 (MQ=255)
                                                                                                cacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGa                                                <  2:142996/149‑1 (MQ=255)
                                                                                                  cccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAg                                               >  5:358096/1‑148 (MQ=255)
                                                                                                  cccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAg                                               <  6:358096/148‑1 (MQ=255)
                                                                                                     gTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGCACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCt                                           >  5:13280/1‑149 (MQ=255)
                                                                                                        cATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATAGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTc                                         <  6:13280/148‑1 (MQ=255)
                                                                                                         aTCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCt                                       >  7:379609/1‑149 (MQ=255)
                                                                                                          tctcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTc                                      >  7:334703/1‑149 (MQ=255)
                                                                                                           ctcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTc                                      >  4:265906/1‑148 (MQ=255)
                                                                                                             ccACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGCCCCCCGTCGAGCCCTTCCTCgcc                                   >  4:62866/1‑147 (MQ=255)
                                                                                                             ccACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGAc                                   >  5:227175/1‑149 (MQ=255)
                                                                                                              cACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACGGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACg                                  >  1:139020/1‑149 (MQ=255)
                                                                                                              cACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGAc                                   <  3:265906/148‑1 (MQ=255)
                                                                                                              cACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACg                                  >  5:432756/1‑149 (MQ=255)
                                                                                                               aCCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGc                                 >  3:223866/1‑149 (MQ=255)
                                                                                                                ccGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCGGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGcc                                >  1:296279/1‑149 (MQ=255)
                                                                                                                ccGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGcc                                <  8:334703/149‑1 (MQ=255)
                                                                                                                ccGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGcc                                >  6:381518/1‑149 (MQ=255)
                                                                                                                    ccTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                        <  2:352987/89‑1 (MQ=255)
                                                                                                                    ccTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                        >  1:352987/1‑89 (MQ=255)
                                                                                                                    ccTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTc                             >  4:427646/1‑148 (MQ=255)
                                                                                                                    ccTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTc                             >  4:13454/1‑148 (MQ=255)
                                                                                                                     cTTAAAGAAAACAAGTTCGGCATCGCGCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTc                             >  1:5010/1‑148 (MQ=255)
                                                                                                                      ttAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGAc                          >  5:355969/1‑149 (MQ=255)
                                                                                                                      ttAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGAc                          >  7:160551/1‑149 (MQ=255)
                                                                                                                        aaaGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCg                        >  2:112962/1‑149 (MQ=255)
                                                                                                                          aGAAAACAAGTTCGTCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCg                                     <  4:134086/134‑1 (MQ=255)
                                                                                                                          aGAAAACAAGTTCGTCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCg                                     >  3:134086/1‑134 (MQ=255)
                                                                                                                           gAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCct                     >  5:15570/1‑149 (MQ=255)
                                                                                                                              aaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctgct                  >  3:188977/1‑149 (MQ=255)
                                                                                                                               aCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGcc                                                                            <  8:161196/90‑1 (MQ=255)
                                                                                                                               aCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGcc                                                                            >  7:161196/1‑90 (MQ=255)
                                                                                                                               aCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCACACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctgctg                 >  4:5623/1‑149 (MQ=255)
                                                                                                                                 aaGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggt               >  3:252493/1‑149 (MQ=255)
                                                                                                                                 aaGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGAGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGCGCCCTTCCTCGACGCCCTCGACCGCCtgctggt               >  6:244758/1‑149 (MQ=255)
                                                                                                                                   gTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc             <  2:23283/149‑1 (MQ=255)
                                                                                                                                   gTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGGCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc             <  4:338997/149‑1 (MQ=255)
                                                                                                                                    ttCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgct            >  8:274056/1‑149 (MQ=255)
                                                                                                                                    ttCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgct            >  2:418163/1‑149 (MQ=255)
                                                                                                                                    ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAAcc                                                                                                    >  2:123361/2‑61 (MQ=25)
                                                                                                                                    ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAAcc                                                                                                    <  1:123361/60‑1 (MQ=25)
                                                                                                                                     tCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgct            <  3:427646/148‑1 (MQ=255)
                                                                                                                                      cGGCATCGGCAATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctg           <  3:122843/148‑1 (MQ=255)
                                                                                                                                      cGGCATCGGCAATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctg           >  4:122843/1‑148 (MQ=255)
                                                                                                                                       ggCATCGCCATTGCCGATGCCGAGGCCATCTCCGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTgctggtgctggt         <  3:62866/149‑1 (MQ=255)
                                                                                                                                       ggCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctgg          <  3:306608/148‑1 (MQ=255)
                                                                                                                                       ggCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTgctggtgctggt         >  6:77139/1‑149 (MQ=255)
                                                                                                                                        gCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTc        <  1:418163/149‑1 (MQ=255)
                                                                                                                                         cATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCg       >  7:5419/1‑149 (MQ=255)
                                                                                                                                          aTCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGa      >  8:253230/1‑149 (MQ=255)
                                                                                                                                           tCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAACGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGAc     >  8:265581/1‑149 (MQ=255)
                                                                                                                                              ccATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGc  >  8:253601/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACCGC  >  NC_002947/5961483‑5961773

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: