Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,035,745 A→G 7.1% intergenic (+34/‑461) PP_4433 → / → dadA‑I amino acid MFS transporter/D‑amino acid dehydrogenase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,035,7450AG 7.1% 699.7 / 27.8 225intergenic (+34/‑461)PP_4433/dadA‑Iamino acid MFS transporter/D‑amino acid dehydrogenase small subunit
Reads supporting (aligned to +/- strand):  ref base A (109/100);  new base G (7/9);  total (116/109)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 226 total aligned reads displayed.

GCACCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGC  >  NC_002947/5035598‑5035894
                                                                                                                                                   |                                                                                                                                                     
 cACCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATc                                                                                                                                                    <  1:460095/150‑1 (MQ=255)
  aCCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                   <  1:3947073/150‑1 (MQ=255)
  aCCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                   <  2:6606970/150‑1 (MQ=255)
  aCCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                   <  1:3097048/150‑1 (MQ=255)
  aCCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATACGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                   >  1:3678447/1‑150 (MQ=255)
  aCCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATCGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                   <  2:3898099/150‑1 (MQ=255)
     tGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCa                                                                                                                                                <  2:870573/150‑1 (MQ=255)
      gCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCAc                                                                                                                                               <  1:6289063/150‑1 (MQ=255)
        ctactaCATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCAccc                                                                                                                                             >  2:1416160/1‑150 (MQ=255)
         tactacATCCTGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCa                                                                                                                                            <  1:5432921/150‑1 (MQ=255)
          actacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAg                                                                                                                                           >  1:4006999/1‑150 (MQ=255)
           ctacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAgg                                                                                                                                          <  1:5022002/150‑1 (MQ=255)
           ctacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAgg                                                                                                                                          <  1:6230205/150‑1 (MQ=255)
           ctacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAgg                                                                                                                                          <  1:6753033/150‑1 (MQ=255)
            tacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAGGc                                                                                                                                         <  2:2259232/150‑1 (MQ=255)
             acATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAGGCt                                                                                                                                        <  2:5741710/150‑1 (MQ=255)
                tcatGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGtt                                                                                                                                     <  1:203008/150‑1 (MQ=255)
                     gTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGaa                                                                                                                                <  2:5339846/150‑1 (MQ=255)
                          aTCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTcc                                                                                                                           >  2:542188/1‑150 (MQ=255)
                            cAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCAt                                                                                                                         <  2:5619796/150‑1 (MQ=255)
                                  ttGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctcc                                                                                                                   >  2:1170084/1‑150 (MQ=255)
                                  ttGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctcc                                                                                                                   >  2:1793558/1‑150 (MQ=255)
                                  ttGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctcc                                                                                                                   >  2:6221464/1‑150 (MQ=255)
                                  ttGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctcc                                                                                                                   >  1:1606886/1‑150 (MQ=255)
                                    gggATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCgg                                                                                                                 <  1:1211402/150‑1 (MQ=255)
                                      gATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGc                                                                                                               >  2:1965908/1‑150 (MQ=255)
                                      gATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGc                                                                                                               >  1:5151370/1‑150 (MQ=255)
                                               gttgttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGAcc                                                                                                      >  2:6062171/1‑150 (MQ=255)
                                                 tgttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCgg                                                                                                    <  2:2189204/150‑1 (MQ=255)
                                                    tCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTa                                                                                                 >  2:1957179/1‑150 (MQ=255)
                                                      tttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCg                                                                                                                  <  2:1350502/131‑1 (MQ=255)
                                                      tttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCg                                                                                                                  >  1:1350502/1‑131 (MQ=255)
                                                      tttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACg                                                                                               >  1:6184628/1‑150 (MQ=255)
                                                      tttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACg                                                                                               >  1:601365/1‑150 (MQ=255)
                                                       ttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGt                                                                                              <  2:3341340/150‑1 (MQ=255)
                                                        tgtgATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTc                                                                                             >  2:3853438/1‑150 (MQ=255)
                                                          tgATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGa                                                                                           <  2:6565621/150‑1 (MQ=255)
                                                             tCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGatat                                                                                        >  2:5692477/1‑150 (MQ=255)
                                                             tCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGatat                                                                                        >  2:3604762/1‑150 (MQ=255)
                                                               gagaCTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATaa                                                                                      >  1:6341953/1‑150 (MQ=255)
                                                                agaCTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAag                                                                                     <  2:5235977/150‑1 (MQ=255)
                                                                 gaCTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagc                                                                                    <  2:4916266/150‑1 (MQ=255)
                                                                   cTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcag                                                                                  <  2:4230940/150‑1 (MQ=255)
                                                                    tAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagc                                                                                 >  1:4772082/1‑150 (MQ=255)
                                                                    tAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagc                                                                                 >  2:5333301/1‑150 (MQ=255)
                                                                    tAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagc                                                                                 >  2:1612740/1‑150 (MQ=255)
                                                                     aaaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACTGGCTACGTCGATATAagcagca                                                                                <  2:636462/150‑1 (MQ=255)
                                                                     aaaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagca                                                                                <  1:4462746/150‑1 (MQ=255)
                                                                     aaaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagca                                                                                <  1:4060738/150‑1 (MQ=255)
                                                                      aaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagcag                                                                               >  1:5347022/1‑150 (MQ=255)
                                                                        ggTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAgcg                                                                             >  2:2588974/1‑150 (MQ=255)
                                                                             cTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGatat                                                                                        >  1:3337889/1‑134 (MQ=255)
                                                                             cTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCa                                                                        >  2:6784217/1‑150 (MQ=255)
                                                                              tCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCAc                                                                       >  2:6650764/1‑150 (MQ=255)
                                                                                 aTTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACcct                                                                    >  1:1009646/1‑150 (MQ=255)
                                                                                  ttCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACcctt                                                                   >  1:5845289/1‑150 (MQ=255)
                                                                                     gcgGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctc                                                                                                                    <  2:4495058/98‑1 (MQ=255)
                                                                                     gcgGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctc                                                                                                                    >  1:4495058/1‑98 (MQ=255)
                                                                                     gcgGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACccttcct                                                                <  2:396866/150‑1 (MQ=255)
                                                                                        gggATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCtttt                                                             <  1:4325327/150‑1 (MQ=255)
                                                                                          gATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTa                                                           <  2:4143974/150‑1 (MQ=255)
                                                                                          gATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTa                                                           <  2:2774686/150‑1 (MQ=255)
                                                                                           aTGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTAt                                                          >  1:1714490/1‑150 (MQ=255)
                                                                                             ggAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTc                                                        <  1:5621848/150‑1 (MQ=255)
                                                                                              gAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCt                                                       >  2:477036/1‑150 (MQ=255)
                                                                                                aTACCTGGGGCAAAATAACCTAATAGGCAGATTTTATAGAGCCCGGGTGCTAATCAGCACCAAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAATCCGACCGGCTACGTCGATATAAGCAGCATCGATCCACCCTTCCTTTTTATTCTTa                                                     <  1:2503598/150‑1 (MQ=255)
                                                                                                           aaaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCtt                                          >  2:758445/1‑150 (MQ=255)
                                                                                                           aaaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCtt                                          >  1:2798373/1‑150 (MQ=255)
                                                                                                           aaaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCtt                                          >  1:2254297/1‑150 (MQ=255)
                                                                                                            aaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTg                                         <  1:3221170/150‑1 (MQ=255)
                                                                                                             aaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATGAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTgg                                        >  1:6422607/1‑150 (MQ=255)
                                                                                                             aaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCAc                                               >  1:4171388/1‑143 (MQ=255)
                                                                                                             aaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCAc                                               <  2:4171388/143‑1 (MQ=255)
                                                                                                                  ccTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCt                                   <  1:3551548/150‑1 (MQ=255)
                                                                                                                  ccTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCt                                   <  1:229457/150‑1 (MQ=255)
                                                                                                                   cTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTg                                  <  2:3792390/150‑1 (MQ=255)
                                                                                                                    tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTgg                                 >  1:5707319/1‑150 (MQ=255)
                                                                                                                    tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTgg                                 >  2:4576092/1‑150 (MQ=255)
                                                                                                                    tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTgg                                 >  2:5024477/1‑150 (MQ=255)
                                                                                                                     aaTAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGc                                >  1:607811/1‑150 (MQ=255)
                                                                                                                      aTAGGCAGTTGGTATAGATCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGcc                               <  1:6735191/150‑1 (MQ=255)
                                                                                                                          gCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGa                           <  1:2181004/150‑1 (MQ=255)
                                                                                                                            aGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATa                         >  1:937055/1‑150 (MQ=255)
                                                                                                                            aGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGAGCTCCACCCTTCCTTTTTATTCTTATATCACCTCTTGGCATCTGGCCCTGATa                         >  2:6304974/1‑150 (MQ=255)
                                                                                                                               tGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACa                      <  2:3580550/150‑1 (MQ=255)
                                                                                                                                ggTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAg                     <  2:2607082/150‑1 (MQ=255)
                                                                                                                                  tataGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAgtg                   >  2:3306805/1‑150 (MQ=255)
                                                                                                                                    taGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTa                                                     >  1:4551759/1‑114 (MQ=255)
                                                                                                                                    taGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTc                 <  1:4405186/150‑1 (MQ=255)
                                                                                                                                    taGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTc                 >  1:964205/1‑150 (MQ=255)
                                                                                                                                       agCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTc              <  1:1876370/150‑1 (MQ=255)
                                                                                                                                         ccTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGc            >  1:1850778/1‑150 (MQ=255)
                                                                                                                                          cTGGGTGCTGATTAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCt           >  1:1180867/1‑150 (MQ=255)
                                                                                                                                          cTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCt           <  1:859935/150‑1 (MQ=255)
                                                                                                                                           tGGGTGCTGATTAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTa          >  2:2829795/1‑150 (MQ=255)
                                                                                                                                           tGGGTGCTGATTAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTa          >  1:5399405/1‑150 (MQ=255)
                                                                                                                                              gTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTAtct       <  1:2436995/150‑1 (MQ=255)
                                                                                                                                                 cTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCtt    <  1:6641282/150‑1 (MQ=255)
                                                                                                                                                  tAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTg   >  1:5552602/1‑150 (MQ=255)
                                                                                                                                                   aaTCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGc  <  2:918819/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCACCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGC  >  NC_002947/5035598‑5035894

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