Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,613 G→T 6.2% C253F (TGC→TTC)  ddlB → D‑alanine‑‑D‑alanine ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,526,6130GT 6.2% 255.2 / 3.1 79C253F (TGC→TTC) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base T (3/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGTGCTGCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTC  >  NC_002947/1526464‑1526761
                                                                                                                                                      |                                                                                                                                                    
ccAGGTGCTGCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTg                                                                                                                                                      <  2:1066900/150‑1 (MQ=255)
                    cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAg                                                                                                                                  >  1:1608326/1‑150 (MQ=255)
                     ccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAgg                                                                                                                                 >  2:763552/1‑150 (MQ=255)
                         ggTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                             >  1:1429542/1‑150 (MQ=255)
                         ggTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                             <  2:1429542/150‑1 (MQ=255)
                               ccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCgg                                                                                                                       >  2:1042736/1‑150 (MQ=255)
                                    cGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGa                                                                                                                  >  1:995553/1‑150 (MQ=255)
                                    cGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGa                                                                                                                  <  2:1153742/150‑1 (MQ=255)
                                    cGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGa                                                                                                                  <  1:1410796/150‑1 (MQ=255)
                                      tgtTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACg                                                                                                                <  2:1355547/150‑1 (MQ=255)
                                        ttCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTg                                                                                                              >  2:1390001/1‑150 (MQ=255)
                                             cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                         <  1:1451209/150‑1 (MQ=255)
                                             cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                         <  2:1608326/150‑1 (MQ=255)
                                             cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                         <  1:1543138/150‑1 (MQ=255)
                                              aaCTATGACGCCAAGGACATCTTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAg                                                                                                        >  2:1257693/2‑150 (MQ=255)
                                                 tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTTCGATGCCAATGGCATCgcag                                                                                                                                  >  1:379118/1‑118 (MQ=255)
                                                 tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTTCGATGCCAATGGCATCGCAg                                                                                                                                  <  2:379118/121‑1 (MQ=255)
                                                 tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGAc                                                                                                     <  1:763552/150‑1 (MQ=255)
                                                 tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCACGTGAGGGCTGGGTAGAGTCAAGCACCCGGAAAGGGTGGATCCTGCGGGGCGCGCCTTGGATGGCATGGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGAc                                                                                                     >  2:1005169/1‑150 (MQ=255)
                                                   tGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGa                                                                                                   <  2:1239629/150‑1 (MQ=255)
                                                     aCGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGc                                                                                                 <  2:609285/150‑1 (MQ=255)
                                                      cGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCa                                                                                                >  2:698076/1‑150 (MQ=255)
                                                         cAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                             >  2:1613230/1‑150 (MQ=255)
                                                            gTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCg                                                                                          <  2:1289758/150‑1 (MQ=255)
                                                            gTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCg                                                                                          >  2:1131783/1‑150 (MQ=255)
                                                             aaCATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCgg                                                                                         <  2:393445/149‑1 (MQ=255)
                                                                aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCg                                                                                          >  1:676025/1‑146 (MQ=255)
                                                                aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCg                                                                                          <  2:676025/146‑1 (MQ=255)
                                                                 tCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                     <  2:961747/150‑1 (MQ=255)
                                                                 tCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                     <  1:1042736/150‑1 (MQ=255)
                                                                   gTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                   <  1:698076/150‑1 (MQ=255)
                                                                   gTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                   >  1:216006/1‑150 (MQ=255)
                                                                       aCGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctgc                                                                               >  2:1148277/1‑150 (MQ=255)
                                                                        cGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                     <  2:222167/143‑1 (MQ=255)
                                                                        cGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                     >  1:222167/1‑143 (MQ=255)
                                                                         gATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctg                                                                                >  2:840201/1‑147 (MQ=255)
                                                                         gATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctg                                                                                <  1:840201/147‑1 (MQ=255)
                                                                           tACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgggg                                                                                                                           >  2:30480/1‑102 (MQ=255)
                                                                           tACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGa                                                                           >  2:132224/1‑150 (MQ=255)
                                                                           tACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTCATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgggg                                                                                                                           <  1:30480/102‑1 (MQ=255)
                                                                              ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTg                                                                                                                                           >  1:1497052/1‑83 (MQ=255)
                                                                              ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTg                                                                                                                                           <  2:1497052/83‑1 (MQ=255)
                                                                              ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGa                                                                                                             <  2:176679/113‑1 (MQ=255)
                                                                              ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGa                                                                                                             >  1:176679/1‑113 (MQ=255)
                                                                              ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                        <  1:1131783/150‑1 (MQ=255)
                                                                              ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                        >  2:146930/1‑150 (MQ=255)
                                                                               cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTTCGATGCCAATGGCATCgcag                                                                                                                                  >  2:513183/1‑88 (MQ=39)
                                                                               cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTTCGATGCCAATGGCATCGCAg                                                                                                                                  <  1:513183/91‑1 (MQ=39)
                                                                                aGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCa                                                                      >  1:263437/1‑150 (MQ=255)
                                                                                      gCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCg                                                                >  1:285975/1‑150 (MQ=255)
                                                                                             gTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                   >  1:495169/1‑124 (MQ=255)
                                                                                             gTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                   <  2:495169/124‑1 (MQ=255)
                                                                                             gTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAgg                                                         <  1:52859/150‑1 (MQ=255)
                                                                                               gCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCa                                                       <  1:1148277/150‑1 (MQ=255)
                                                                                                 ggagtGGATAGCGTCAAGGAACAGGAACTGATATACCTGACCGCGGGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTCATGCAGGACGAGCAGCGCCGGTTGTGGCTGCTCCAACCCCACACCGCACCAGGCacg                                                     >  2:104223/5‑148 (MQ=255)
                                                                                                   gCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGa                                                                                                                                                   <  1:149900/54‑1 (MQ=255)
                                                                                                   gCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGa                                                                                                                                                   >  2:149900/1‑54 (MQ=255)
                                                                                                        aTAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAgg                                                                                              <  2:1301920/102‑1 (MQ=255)
                                                                                                        aTAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAgg                                                                                              >  1:1301920/1‑102 (MQ=255)
                                                                                                             gTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATAACAGATCATAGc                                         <  2:263437/150‑1 (MQ=255)
                                                                                                              tCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGa                                                    <  2:1400397/138‑1 (MQ=255)
                                                                                                              tCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGa                                                    >  1:1400397/1‑138 (MQ=255)
                                                                                                               cAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCt                                       <  2:995553/150‑1 (MQ=255)
                                                                                                                     acaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGcc                                 >  2:951828/1‑150 (MQ=255)
                                                                                                                       aggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCa                               >  1:825738/1‑150 (MQ=255)
                                                                                                                         gaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCagg                             <  2:216006/150‑3 (MQ=255)
                                                                                                                         gaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATg                             >  2:305026/1‑150 (MQ=255)
                                                                                                                         gaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATg                             <  1:132224/150‑1 (MQ=255)
                                                                                                                          aacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATgg                            <  1:951828/150‑1 (MQ=255)
                                                                                                                           acTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATggc                           <  1:1535393/150‑1 (MQ=255)
                                                                                                                               aTCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATg                             <  1:1131278/144‑1 (MQ=255)
                                                                                                                               aTCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATg                             >  2:1131278/1‑144 (MQ=255)
                                                                                                                                   aCCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTg                   >  1:1610548/1‑150 (MQ=255)
                                                                                                                                            gcgcgcgcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGc           <  2:285975/150‑1 (MQ=255)
                                                                                                                                              gcgcgcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCt         <  1:1005169/150‑1 (MQ=255)
                                                                                                                                                    cTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt                                                      <  2:1475985/99‑1 (MQ=255)
                                                                                                                                                    cTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt                                                      >  1:1475985/1‑99 (MQ=255)
                                                                                                                                                     tGCGATGCCATTGGCATCGAATGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCt                                       <  1:17183/113‑1 (MQ=255)
                                                                                                                                                     tGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCt                                       >  2:17183/1‑113 (MQ=255)
                                                                                                                                                     tGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTc  >  1:351514/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCAGGTGCTGCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTC  >  NC_002947/1526464‑1526761

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