Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,785,057 G→A 5.1% intergenic (+25/+81) PP_1587 → / ← dapC NhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,785,0570GA 5.1% 375.0 / 6.1 118intergenic (+25/+81)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (48/64);  new base A (4/2);  total (52/66)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 123 total aligned reads displayed.

TGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAG  >  NC_002947/1784911‑1785206
                                                                                                                                                  |                                                                                                                                                     
tGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGgcgcgc                                                                                                                                                    >  2:3638032/1‑150 (MQ=255)
tGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgc                                                                                                                                                    <  1:5056374‑M1/150‑4 (MQ=255)
      cgcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTc                                                                                                                                              <  1:1148869/150‑1 (MQ=255)
       gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGCTGAGGGGGTGCCGCAAGTACTGGGTGAGCTGGAGATGAGTGAGGCGAACCTGGGTCCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaa                                                                                                                                             <  2:3921725‑M1/150‑11 (MQ=255)
       gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGcgaa                                                                                                                                             <  1:143486/150‑4 (MQ=255)
       gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGcgaa                                                                                                                                             <  2:4930782/150‑4 (MQ=255)
       gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                             <  1:206495/150‑1 (MQ=255)
        cAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTccc                                                                                                                                            <  1:4119437/150‑1 (MQ=255)
        cAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTccc                                                                                                                                            >  1:3487778/1‑150 (MQ=255)
        cAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTccc                                                                                                                                            <  2:783696/150‑1 (MQ=255)
        cAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTccc                                                                                                                                            >  1:3487852/1‑150 (MQ=255)
           cGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCcaca                                                                                                                                         <  2:3524612/150‑1 (MQ=255)
             gTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAgg                                                                                                                                       <  1:2956640/150‑1 (MQ=255)
              ttGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAACCCCTTCGCGGGCAcgcccgcgaagggccttg                                                                                                                                     <  2:2730293‑M1/150‑19 (MQ=255)
              ttGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAACCCCTTCGCGGGCAcgcccgcgaagggccttg                                                                                                                                     <  2:2730231‑M1/150‑19 (MQ=255)
                    cTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcg                                                                                                                                >  2:3254530/1‑150 (MQ=255)
                     tGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCACTCCCACAGGTActgcgc                                                                                                                               >  2:1637133/1‑150 (MQ=255)
                     tGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCACTCCCACAGGTActgcgc                                                                                                                               >  2:1637440/1‑150 (MQ=255)
                      gTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGTAaagcgcg                                                                                                                              >  1:1491008/1‑149 (MQ=255)
                         cAACCTGCCTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTTCTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgcggcg                                                                                                                           <  2:5347286/150‑1 (MQ=255)
                         cAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgctgcg                                                                                                                           <  2:3263095/150‑1 (MQ=255)
                          aaCCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggtttt                                                                                                                          <  2:1114950‑M1/150‑30 (MQ=255)
                           aCCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTgcgctgcgct                                                                                                                         <  1:1433687/150‑1 (MQ=255)
                            ccTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAcagcgctgatttc                                                                                                                        <  2:4550392/150‑6 (MQ=255)
                             cTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCttt                                                                                                                       >  1:4201766/1‑149 (MQ=255)
                             cTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTg                                                                                                                       >  2:4215472/1‑150 (MQ=255)
                             cTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACTGCGCTGCGCTTg                                                                                                                       >  1:1332985/1‑150 (MQ=255)
                              tGCATCGTCAGCTCCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGa                                                                                                                      <  2:3541973/150‑1 (MQ=255)
                              tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGa                                                                                                                      <  1:5614180/150‑1 (MQ=255)
                              tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGa                                                                                                                      <  2:3707885/150‑1 (MQ=255)
                              tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGa                                                                                                                      >  2:4602716/1‑150 (MQ=255)
                              tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGa                                                                                                                      <  1:1266325/150‑1 (MQ=255)
                                 aTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAAt                                                                                                                   <  1:1747718/150‑1 (MQ=255)
                                   cGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCa                                                                                                                 <  2:3694778/150‑1 (MQ=255)
                                    gTCAGCACCTGGGGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACCCGTACTGCGCTGCGCGGGATATCAg                                                                                                                >  1:1106445/1‑150 (MQ=255)
                                     tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgtccgcgaaggggttgatcagcggttttgcaggaaggcg                                                                                                               <  2:1544691‑M1/150‑41 (MQ=255)
                                     tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGGACAGCGCGGCGCATGaaaaaacc                                                                                                               >  1:4027689/1‑145 (MQ=255)
                                     tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgc                                                                                                               <  1:3672582/150‑1 (MQ=255)
                                        gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCt                                                                                                            <  1:979835/150‑1 (MQ=255)
                                         cACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTg                                                                                                           <  1:3801424/150‑1 (MQ=255)
                                           ccTGGTGGGGGATGTGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                         <  1:4461628/150‑1 (MQ=255)
                                           ccTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                         <  1:3483941/150‑1 (MQ=255)
                                           ccTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                         <  1:3483959/150‑1 (MQ=255)
                                            cTGGTGGGGGATGTGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTAAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAc                                                                                                        <  1:2796099/150‑1 (MQ=255)
                                             tggtggGGGATGTGGTGGTGCGGCAAGTACTGTGTGAGCTGGATATGAGTGAGTCGAAGCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcagtaaggcgcgaatgcg                                                                                                       <  1:4524373‑M1/150‑49 (MQ=255)
                                              ggtggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCt                                                                                                      <  1:2365274/150‑1 (MQ=255)
                                                tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                    <  2:3099823/150‑1 (MQ=255)
                                                tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                    <  2:3852985/150‑1 (MQ=255)
                                                tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                    <  1:2997955/150‑1 (MQ=255)
                                                tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                    <  1:149064/150‑1 (MQ=255)
                                                 gggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgtg                                                                                                   <  1:4645279/150‑1 (MQ=255)
                                                  ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                  <  1:4347749/150‑1 (MQ=255)
                                                  ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                  <  1:4416004/150‑1 (MQ=255)
                                                  ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                  <  1:4347784/150‑1 (MQ=255)
                                                  ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                  <  2:1073279/150‑1 (MQ=25)
                                                    ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                                <  1:3638032/150‑1 (MQ=255)
                                                    ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                                <  1:3770089/150‑1 (MQ=255)
                                                    ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                                <  1:1039171/150‑1 (MQ=255)
                                                    ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGa                                                                                                <  2:47394/150‑1 (MQ=21)
                                                    ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGa                                                                                                <  1:3351957/150‑1 (MQ=21)
                                                    ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGTACAGCGCTGCGCTTGAAAAAAGCGCAGTAACCGTgggg                                                                                                >  1:5373592/1‑149 (MQ=255)
                                                        gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAATGCGCTGCGCTTGAAAACAGCGCTGTACCTGTGGGAGCgg                                                                                            >  1:5583064/1‑150 (MQ=255)
                                                          ggtggtGCGGCAAGTACTGGGTGAGCTGGATATGTGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGc                                                                                          <  1:5129381/150‑1 (MQ=21)
                                                                 cGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCGTGAAATAAGCGCTGTACCTGTGGGAGCGGGCGTGcccc                                                                                   >  1:3188627/1‑149 (MQ=255)
                                                                   gCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcgc                                                                                                             >  2:911815/1‑122 (MQ=255)
                                                                   gCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcgc                                                                                                             <  1:911815/122‑1 (MQ=255)
                                                                        tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                            >  2:3641592/1‑150 (MQ=255)
                                                                          cTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGCTTAAGCGCTTCGCTTTGAAACAGCGCTTTACCAGTGGGATCGGGCGTGCCCGCGAAGGGGt                                                                          >  2:731875/1‑150 (MQ=255)
                                                                                  gCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGc                                                                  >  2:4664975/1‑150 (MQ=255)
                                                                                   cTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGTACTGCGCTGCGCGTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCg                                                                 >  1:3257400/1‑150 (MQ=255)
                                                                                   cTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTTCTGCGGTGGGCTTGATATGAGCGCGGTAGCTGTGGGGGCGGGGGTGGCCGCGTAGGGGTTGGTGAGCg                                                                 >  1:2332834/1‑150 (MQ=255)
                                                                                   cTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGCAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCg                                                                 >  2:5526469/1‑150 (MQ=255)
                                                                                   cTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCGTGAAATAAGCGCGGTTCCGGTGGGGGCGGGCGGGCCCGCGAAGGGGGTGTTCAGCg                                                                 >  2:2607961/1‑150 (MQ=37)
                                                                                      gATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCATTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATAAGCGGtt                                                              >  2:2626138/1‑150 (MQ=255)
                                                                                          tgagtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCCAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgccggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccat                                                          <  2:5071220‑M1/150‑94 (MQ=255)
                                                                                          tgagtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccat                                                          <  2:1817671‑M1/150‑94 (MQ=255)
                                                                                           gagtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatg                                                         >  1:289605‑M1/1‑56 (MQ=255)
                                                                                           gagtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACGGCGCTGCGCTTGAAATCAGCGCTGTACCTGGGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCACCGGTGTTGcc                                                         >  2:323602/1‑149 (MQ=255)
                                                                                              tgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTCCTGCGCTGCGCTTGGAATCAGCGCGGTCCCGGGGGGCGCGGGCGGGCCCGCGACGGGGGTGATCAGCGGGTTTGCAGGa                                                      >  2:5494188/1‑150 (MQ=255)
                                                                                                  gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACGGCGCTGCGCGTGAAATCAGCCCTGTACCTGTGGGAGCGGGCGGGCCCGCGGAGGGGGTGATCAGCGGGTTTGCCGGAAGgc                                                  >  2:675460/1‑150 (MQ=38)
                                                                                                       ccTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccg                                             <  1:2688765‑M1/150‑107 (MQ=255)
                                                                                                            gTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCGTTAAATCAGCGCTGTAACTGTGGGCGCGGGCGTGCCCGCGAAGGGGTTGTTCAGCGGGTTTGCAGGAAGGCGCGAATgcgc                                        >  2:5263736/1‑150 (MQ=38)
                                                                                                             ttCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGACGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGACTGCGCt                                       >  2:3678602/1‑150 (MQ=255)
                                                                                                                  tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCGGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTccgcc                                  >  2:1772601/1‑150 (MQ=255)
                                                                                                                  tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACAGCGCTGCGCGTTAAATCAGCGCGGTACCAGTGGGGGCGGGCGTGCCCCCGAAGGGGTTGTTCAGCCGTTTTGCAGGAAGGCGCCCAAGCGCAccgcc                                  >  2:1126766/1‑150 (MQ=255)
                                                                                                                      gTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCcgccgcgt                              >  2:1243190/1‑148 (MQ=255)
                                                                                                                      gTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACTGCGCTGCGCTTGAAATCAGCGCGGTACCTGTGGGCGCGGGCGTGCCCGCGAAGGGGGTGGTCAGCGGTTTTGCCGGAAGGCGCGCATGCGCGCCGCCGCtg                              >  2:3541883/1‑149 (MQ=37)
                                                                                                                            ccTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatcc                        >  1:1894142‑M1/1‑23 (MQ=255)
                                                                                                                            ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGTAAGGCGCGAATGCGCTCCGCCGCTTCGATGc                        >  2:4076938/1‑150 (MQ=255)
                                                                                                                                    ctgggggcgcgggcgtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacg                <  2:1333091‑M1/143‑136 (MQ=255)
                                                                                                                                    cggggggcgcgggcgtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacg                <  2:1332985‑M1/143‑136 (MQ=255)
                                                                                                                                       ccTTCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcga             >  2:4975552‑M1/1‑12 (MQ=255)
                                                                                                                                        cTTCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgactgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgac            >  2:4995664‑M1/1‑11 (MQ=255)
                                                                                                                                        cTTCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgac            >  1:4872393‑M1/1‑11 (MQ=255)
                                                                                                                                                gCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacaggtacga    >  2:2792296‑M1/1‑3 (MQ=255)
                                                                                                                                                gCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacaggtacga    <  2:5373592‑M1/150‑148 (MQ=255)
                                                                                                                                                gCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacaggtacga    <  2:3257400‑M1/150‑148 (MQ=255)
                                                                                                                                                gCGtgcccgcgaagggggtgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacaggtacga    >  1:5357953‑M1/1‑3 (MQ=255)
                                                                                                                                                  gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTAGGCCAGCGGTGCAACCAg  >  1:4062892/1‑150 (MQ=255)
                                                                                                                                                  gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCAGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTTCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg  >  1:2156064/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAG  >  NC_002947/1784911‑1785206

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