Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,030,799 C→T 9.0% intergenic (+37/‑458) 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,030,7990CT 9.0% 214.6 / 10.5 78intergenic (+37/‑458)PP_4433/dadA‑Iamino acid MFS transporter/D‑amino acid dehydrogenase small subunit
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base T (4/3);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

ACCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGCC  >  NC_002947/5030651‑5030946
                                                                                                                                                    |                                                                                                                                                    
aCCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                     >  1:450290/1‑150 (MQ=255)
   tGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAGTCAGCa                                                                                                                                                  >  2:1682149/1‑150 (MQ=255)
     ccTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGAAGATGGTATAGAGCCTGGAAGCTAATCATCAcc                                                                                                                                                >  2:1342889/1‑150 (MQ=255)
       tactacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCa                                                                                                                                              >  1:833190/1‑150 (MQ=255)
         ctacATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAgg                                                                                                                                            <  2:1547380/150‑1 (MQ=255)
         ctacATCATGGTTACATCCGCACTTAAGCGTATTAAGTTTGCCAAGTTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAgg                                                                                                                                            <  2:1231808/150‑1 (MQ=255)
              tcatGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                     >  1:2093977/1‑136 (MQ=255)
              tcatGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCa                                                                                                                                                     <  2:2093977/136‑1 (MQ=255)
              tcatGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGtt                                                                                                                                       <  2:2250925/150‑1 (MQ=255)
               catGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTg                                                                                                                                      >  2:1864031/1‑150 (MQ=255)
                atGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTga                                                                                                                                     >  1:114948/1‑149 (MQ=255)
                 tGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGc                                                                                                                                           <  1:829997/143‑1 (MQ=255)
                 tGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGc                                                                                                                                           >  2:829997/1‑143 (MQ=255)
                    ttGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAGGCt                                                                                                                                          >  1:188837/1‑141 (MQ=255)
                    ttGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCAGGCt                                                                                                                                          <  2:188837/141‑1 (MQ=255)
                    ttGCATCAGCACTTGGGATAATTGCGTTGTTCTGTGTGATCGAGACTAAAGGTTGATCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCGAATAGGCAGATGGTATCGAGCCTGGGAACTAATCAGCACCCAGGCTGTTGCGAAg                                                                                                                                 >  2:915180/1‑150 (MQ=255)
                     tGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGc                                                                                                                                <  2:541504/150‑1 (MQ=255)
                           aGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCAtt                                                                                                                          <  1:711148/150‑1 (MQ=255)
                             cACTTGGGATTATTTCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcc                                                                                                                        <  2:1026234/150‑1 (MQ=255)
                               cTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTAGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctc                                                                                                                      <  1:1342889/150‑1 (MQ=255)
                                ttGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctcc                                                                                                                     >  1:2009684/1‑150 (MQ=255)
                                ttGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGCGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcctcc                                                                                                                     >  2:662629/1‑150 (MQ=255)
                                         ttGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGcc                                                                                                            <  2:450290/150‑1 (MQ=255)
                                              ttgttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGtt                                                                                                                                       >  1:1495810/1‑118 (MQ=255)
                                              ttgttCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGtt                                                                                                                                       <  2:1495810/118‑1 (MQ=255)
                                                    tttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCa                                                                                                                                              >  2:768488/1‑105 (MQ=255)
                                                    tttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCa                                                                                                                                              <  1:768488/105‑1 (MQ=255)
                                                    tttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCgg                                                                                                                   <  1:210761/132‑1 (MQ=255)
                                                    tttGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCgg                                                                                                                   >  2:210761/1‑132 (MQ=255)
                                                         gATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAGGCCGACCGGCTACGTCGat                                                                                            <  2:278431/150‑1 (MQ=255)
                                                         gATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGat                                                                                            <  2:1879026/150‑1 (MQ=255)
                                                             gagaCTAAGGGTTGATTGATTTGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATaa                                                                                        <  2:386624/150‑1 (MQ=255)
                                                                  tAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGat                                                                                            >  1:183816/1‑141 (MQ=255)
                                                                  tAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGat                                                                                            <  2:183816/141‑1 (MQ=255)
                                                                  tAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagc                                                                                  >  2:1876157/1‑150 (MQ=255)
                                                                    aaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCg                                                                                                       >  1:1208603/1‑128 (MQ=255)
                                                                    aaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCg                                                                                                       <  2:1208603/128‑1 (MQ=255)
                                                                    aaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGata                                                                                           <  2:478586/140‑1 (MQ=255)
                                                                    aaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGata                                                                                           >  1:478586/1‑140 (MQ=255)
                                                                    aaGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGACTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATCTAagcagcag                                                                                >  2:896687/1‑150 (MQ=255)
                                                                        ttGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAgcgat                                                                            >  1:2089434/1‑148 (MQ=255)
                                                                         tGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTc                                                                           <  2:747008/150‑1 (MQ=255)
                                                                           cTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCa                                                                         <  2:114948/150‑1 (MQ=255)
                                                                            tCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCAc                                                                        <  2:2009684/150‑1 (MQ=255)
                                                                            tCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCAc                                                                        >  2:191108/1‑150 (MQ=255)
                                                                                  cgcgGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCTCTCCACccttcc                                                                  <  1:1864031/150‑1 (MQ=255)
                                                                                        gATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGc                                                                                                             >  1:1506248/1‑102 (MQ=255)
                                                                                        gATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGc                                                                                                             <  2:1506248/102‑1 (MQ=255)
                                                                                        gATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTa                                                            >  1:1924100/1‑150 (MQ=255)
                                                                                                  cTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCa                                                                                                                                              >  1:2122849/1‑59 (MQ=37)
                                                                                                  cTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTGATTAGCACCCa                                                                                                                                              <  2:2122849/59‑1 (MQ=37)
                                                                                                  cTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTAtctc                                                  >  2:2386225/1‑150 (MQ=255)
                                                                                                         aaaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACctc                                             >  1:2085913/1‑148 (MQ=255)
                                                                                                         aaaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACctc                                             <  2:2085913/148‑1 (MQ=255)
                                                                                                         aaaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCtt                                           <  2:1248067/150‑1 (MQ=255)
                                                                                                          aaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTc                                                         >  1:2391817/1‑135 (MQ=255)
                                                                                                          aaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTc                                                         <  2:2391817/135‑1 (MQ=255)
                                                                                                          aaaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTg                                          <  1:2386225/150‑1 (MQ=255)
                                                                                                           aaTAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTgg                                         >  2:1276658/1‑150 (MQ=255)
                                                                                                                ccTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAg                                                                                                                                             <  2:1580195/46‑1 (MQ=255)
                                                                                                                ccTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAg                                                                                                                                             >  1:1580195/1‑46 (MQ=255)
                                                                                                                  tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAg                                                                                                                                             <  2:1771328/44‑1 (MQ=255)
                                                                                                                  tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAg                                                                                                                                             >  1:1771328/1‑44 (MQ=255)
                                                                                                                  tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcc                                                                                                                        <  1:1920986/65‑1 (MQ=255)
                                                                                                                  tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATtcc                                                                                                                        >  2:1920986/1‑65 (MQ=255)
                                                                                                                  tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagca                                                                                 <  1:1658795/103‑1 (MQ=255)
                                                                                                                  tAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAagcagca                                                                                 >  2:1658795/1‑103 (MQ=255)
                                                                                                                        gCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGa                            >  1:1509988/1‑150 (MQ=255)
                                                                                                                          aGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGa                            <  1:1560324/148‑1 (MQ=255)
                                                                                                                          aGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGa                            >  2:1560324/1‑148 (MQ=255)
                                                                                                                          aGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGa                            >  1:1274116/1‑148 (MQ=255)
                                                                                                                          aGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCCGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGa                            <  2:1274116/148‑1 (MQ=255)
                                                                                                                                 ataGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAgtgt                   <  1:191108/150‑1 (MQ=255)
                                                                                                                                         tGGGTGCTGATTAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCCAGTGTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTa           >  1:1956870/1‑150 (MQ=255)
                                                                                                                                         tGGGTGCTGATTAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTa           >  2:231070/1‑150 (MQ=255)
                                                                                                                                                  aTCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGcc  >  1:1550387/1‑150 (MQ=255)
                                                                                                                                                   tCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGcc                                >  1:819799/1‑119 (MQ=255)
                                                                                                                                                   tCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGcc                                <  2:819799/119‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCTGCCTACTACATCATGGTTGCATCAGCACTTGGGATTATTGCGTTGTTCTTTGTGATCGAGACTAAAGGTTGCTCGATTCGCGGGGATGGAATACCTGGGGCAAAATAACCTAATAGGCAGATGGTATAGAGCCTGGGTGCTAATCAGCACCCAGGCTGTTGCGAAGCTCCATTCCTCCGGGCAAGCCGACCGGCTACGTCGATATAAGCAGCAGCGCTCCACCCTTCCTTTTTATTCTTATCTCACCTCTTGGCATCTGGCCCTGATAACAGTGTCTTCGCTATCTCTTGCC  >  NC_002947/5030651‑5030946

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