Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,591,970 A→T 5.0% T122S (ACG→TCG)  dapE → N‑succinyl‑diaminopimelate deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,591,9700AT 5.0% 245.0 / 2.9 80T122S (ACG→TCG) dapEN‑succinyl‑diaminopimelate deacylase
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base T (2/2);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGAAGTGGTAGGT  >  NC_000913/2591824‑2592116
                                                                                                                                                  |                                                                                                                                                  
gCGACGCCGATCGTTGGATCAATCCCCCGTTTGACCCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACg                                                                                                                                                  >  5:171998/1‑149 (MQ=255)
 cGACGCCGATCGTTGGATCAATCCCCCGTTTTAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGTCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACgg                                                                                                                                                 <  6:45388/149‑1 (MQ=255)
  gACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACggg                                                                                                                                                <  2:285947/149‑1 (MQ=255)
  gACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACggg                                                                                                                                                >  2:106038/1‑149 (MQ=255)
      ccGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGa                                                                                                                                            <  2:2540/149‑1 (MQ=255)
      ccGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCGCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGa                                                                                                                                            <  8:101360/149‑1 (MQ=255)
        gATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACt                                                                                                                                          <  8:47100/149‑1 (MQ=255)
           cGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTgg                                                                                                                                        <  1:106038/148‑1 (MQ=255)
             ttGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACCACATCCAAACCATACGGGGCGACTGGCAt                                                                                                                                     >  6:178669/1‑149 (MQ=255)
             taggATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGACGCACAACATCCCAACCATACGGGGCGACTGGCAt                                                                                                                                     >  8:200178/3‑149 (MQ=255)
               ggATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGTCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGACGCACAACATCCCAACCATACGGGGCGACTGGCAttt                                                                                                                                   <  7:200178/149‑1 (MQ=255)
                 aTCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAAACATACGGGGCGACTGGCATTTCt                                                                                                                                 >  8:327143/1‑149 (MQ=255)
                    aaTCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCATCCATACGAGGCGACAGGCATTTCTGAt                                                                                                                              >  4:84870/1‑149 (MQ=255)
                    aaTCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGAt                                                                                                                              <  1:356885/149‑1 (MQ=255)
                         cccGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTg                                                                                                                                >  2:146675/1‑142 (MQ=255)
                         cccGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTg                                                                                                                                <  1:146675/142‑1 (MQ=255)
                             tttGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCtga                                                                                                                     <  3:341789/149‑1 (MQ=255)
                                gAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCtgatg                                                                                                                   >  2:183972/1‑148 (MQ=255)
                                   cccACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTc                                                                                                                                  <  4:170625/130‑1 (MQ=255)
                                   cccACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTc                                                                                                                                  >  3:170625/1‑130 (MQ=255)
                                    ccaccaTTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATgaagaa                                                                                                              >  8:89044/1‑149 (MQ=255)
                                    ccaccaTTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATgaagaa                                                                                                              <  2:93914/149‑1 (MQ=255)
                                     caccaTTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGaagaag                                                                                                             >  2:59973/1‑149 (MQ=255)
                                      gccaTTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGc                                                                                                            >  1:48773/2‑149 (MQ=255)
                                      accaTTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGc                                                                                                            >  5:173900/1‑149 (MQ=255)
                                        caTTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCa                                                                                                          <  8:272122/149‑1 (MQ=255)
                                              tGACGGCATGTTATTCTGGCGCTGTGCGGCAGATAAGAAAGCCTAGCTGGCGGCGATGAACGTCGCGGCAGAACGTTTTGTCGCACAAAATCACAACCAAACGGGCCCACTGGCATTTCTCATCACCTCTGATGAAGAACCCAGTGcc                                                                                                     >  8:47179/1‑148 (MQ=255)
                                              tGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGccc                                                                                                    <  6:214552/149‑1 (MQ=255)
                                                  ggCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACaa                                                                                                >  7:50770/1‑149 (MQ=255)
                                                  ggCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACaa                                                                                                >  8:380514/1‑149 (MQ=255)
                                                  ggCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTACTGATCACCTCTGATGAAGAAGCCAGTGCCCACaa                                                                                                >  6:135257/1‑149 (MQ=255)
                                                  ggCATGTTATTCGGGCGCGGGGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACaa                                                                                                <  8:84406/149‑1 (MQ=255)
                                                      tGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCCCAACATCCCAACCATACGGGGCGACTGGCATTTCTGCTCACCTCTGATGAAGAAGCCAGTGCCCACAacgg                                                                                             >  8:95316/1‑148 (MQ=255)
                                                      tGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACACCATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAacggt                                                                                            >  6:87928/1‑149 (MQ=255)
                                                       gTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAacggta                                                                                           >  6:253797/1‑149 (MQ=255)
                                                        ttATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAacggta                                                                                           <  5:135257/148‑1 (MQ=255)
                                                          aTTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAacggtacg                                                                                         <  1:349704/148‑1 (MQ=255)
                                                               ggCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTaaaa                                                                                   <  3:279199/149‑1 (MQ=255)
                                                               ggCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTaaaa                                                                                   <  7:89044/149‑1 (MQ=255)
                                                                  gcgGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAgtc                                                                                >  6:103870/1‑149 (MQ=255)
                                                                            aGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGt                                                                      <  1:183972/149‑1 (MQ=255)
                                                                             gATATGAAAGGCACGCTGGCGGCGATCGTGGAGGCGCCAGAACGTTTTGTCGCACAACATCCCCACCATACGGCGCGACTGGCATATCTCATCCCCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTTAAACTCGTCGAAGCga                                                                      >  1:353941/1‑147 (MQ=255)
                                                                              atatGAAAGGCTCGCGGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTa                                                                    <  6:132238/149‑1 (MQ=255)
                                                                                atGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAgtc                                                                                <  2:52183/135‑1 (MQ=255)
                                                                                atGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAgtc                                                                                >  1:52183/1‑135 (MQ=255)
                                                                                  gAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCtga                                                                                                                     <  3:279185/96‑1 (MQ=255)
                                                                                  gAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCtga                                                                                                                     >  4:279185/1‑96 (MQ=255)
                                                                                    aaGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACATCATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCa                                                              <  4:167002/149‑1 (MQ=255)
                                                                                         tCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTaa                                                         >  4:13095/1‑149 (MQ=255)
                                                                                          cGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAAt                                                        <  5:103870/149‑1 (MQ=255)
                                                                                                   cgATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTATGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGtctc                                               <  7:327143/149‑1 (MQ=255)
                                                                                                          ggtggCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGACGAAGCCAGTGCCCACAACGGTACGGTAAATGTCGTCGAAGCGTTAATGGCACGTATTGAGCGTCTCGATTACt                                        <  6:77519/149‑1 (MQ=255)
                                                                                                            tggCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGc                                      <  6:173900/149‑1 (MQ=255)
                                                                                                               cggcAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTCATGAAGAAGCCAGTGCCCACATCGGTACGGTAAAAGTCGTCGAAGCGTTAATTGCACGTAATGAGCGTCTCGATTACTTCCTg                                   <  8:153318/149‑1 (MQ=255)
                                                                                                                ggcAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTgg                                  <  2:48773/149‑1 (MQ=255)
                                                                                                                ggcAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTgg                                  <  5:253797/149‑1 (MQ=255)
                                                                                                                     aaCGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGc                             >  6:335902/1‑149 (MQ=255)
                                                                                                                      aCGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGt                                                                      <  5:354638/107‑1 (MQ=255)
                                                                                                                      aCGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGt                                                                      >  6:354638/1‑107 (MQ=255)
                                                                                                                      aCGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGc                             >  3:22972/1‑148 (MQ=255)
                                                                                                                         ttttGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAAc                         >  1:115225/1‑149 (MQ=255)
                                                                                                                             gTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTc                     >  5:321816/1‑149 (MQ=255)
                                                                                                                                   ctatctCCCAACCTTTCGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGGCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTa                <  7:278985/143‑1 (MQ=255)
                                                                                                                                   caacaTCCCAACCATACGGGGCGACTGGCATTTCTGAACACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTAt               >  2:257637/1‑149 (MQ=255)
                                                                                                                                    aacaTCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATc              >  5:79684/1‑149 (MQ=255)
                                                                                                                                      cttCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGa            <  7:380514/147‑1 (MQ=255)
                                                                                                                                       aTCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGaa           <  8:50770/149‑1 (MQ=255)
                                                                                                                                       aTCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGaa           <  5:87928/149‑1 (MQ=255)
                                                                                                                                         cccAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGAAGt         >  3:320113/1‑149 (MQ=255)
                                                                                                                                          ccAACCATTCGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGCGCGACTCGATTACTGCCTGGTTGGCGAAGCGGCGAGTATAGatcg         >  8:273794/1‑145 (MQ=255)
                                                                                                                                          ccAACCATTCGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGACTGGTTGGCGAACCGTCGAgca                >  8:278985/1‑139 (MQ=255)
                                                                                                                                          ccAACCATTCGGGGAGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTa                <  7:273794/141‑1 (MQ=255)
                                                                                                                                          ccAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGAAGTg        >  7:277762/1‑149 (MQ=255)
                                                                                                                                           cAACGATACGGGGCGACTGGCATTTCTGAGCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGAAGTgg       <  3:84870/149‑1 (MQ=255)
                                                                                                                                             aCCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGAAGTgg       >  2:277526/1‑147 (MQ=255)
                                                                                                                                              ccATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAACGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGGTGGCGAACCGTCGAGTATCGAAGTGGTAg    <  1:59973/149‑1 (MQ=255)
                                                                                                                                               cATCCGGGTCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGAAGTGGTAgg   <  6:171998/149‑1 (MQ=255)
                                                                                                                                                aTACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGtctc                                               <  7:76093/104‑1 (MQ=255)
                                                                                                                                                aTACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGtctc                                               >  8:76093/1‑104 (MQ=255)
                                                                                                                                                aTACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGCACCGTCGAGTATCGAAGTGGTAGGt  >  4:270780/1‑149 (MQ=255)
                                                                                                                                                  aCGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGt                                 <  5:146638/116‑1 (MQ=255)
                                                                                                                                                  aCGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGt                                 >  6:146638/1‑116 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGACGCCGATCGTTGGATCAATCCCCCGTTTGAACCCACCATTCGTGACGGCATGTTATTCGGGCGCGGTGCGGCAGATATGAAAGGCTCGCTGGCGGCGATGGTGGTGGCGGCAGAACGTTTTGTCGCACAACATCCCAACCATACGGGGCGACTGGCATTTCTGATCACCTCTGATGAAGAAGCCAGTGCCCACAACGGTACGGTAAAAGTCGTCGAAGCGTTAATGGCACGTAATGAGCGTCTCGATTACTGCCTGGTTGGCGAACCGTCGAGTATCGAAGTGGTAGGT  >  NC_000913/2591824‑2592116

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