Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,622 T→A 5.1% I256N (ATT→AAT)  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,6220TA 5.1% 251.5 / 4.3 79I256N (ATT→AAT) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

CCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTG  >  NC_002947/1526484‑1526770
                                                                                                                                           |                                                                                                                                                    
cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAg                                                                                                                                           >  1:1608326/1‑150 (MQ=255)
 ccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAgg                                                                                                                                          >  2:763552/1‑150 (MQ=255)
     ggTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                                      <  2:1429542/150‑1 (MQ=255)
     ggTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                                      >  1:1429542/1‑150 (MQ=255)
           ccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCgg                                                                                                                                >  2:1042736/1‑150 (MQ=255)
                cGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGa                                                                                                                           <  2:1153742/150‑1 (MQ=255)
                cGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGa                                                                                                                           <  1:1410796/150‑1 (MQ=255)
                cGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGa                                                                                                                           >  1:995553/1‑150 (MQ=255)
                  tgtTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACg                                                                                                                         <  2:1355547/150‑1 (MQ=255)
                    ttCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTg                                                                                                                       >  2:1390001/1‑150 (MQ=255)
                         cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                  <  2:1608326/150‑1 (MQ=255)
                         cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                  <  1:1451209/150‑1 (MQ=255)
                         cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                  <  1:1543138/150‑1 (MQ=255)
                          aaCTATGACGCCAAGGACATCTTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAg                                                                                                                 >  2:1257693/2‑150 (MQ=255)
                             tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTTCGATGCCAATGGCATCgcag                                                                                                                                           >  1:379118/1‑118 (MQ=255)
                             tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTTCGATGCCAATGGCATCGCAg                                                                                                                                           <  2:379118/121‑1 (MQ=255)
                             tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGAc                                                                                                              <  1:763552/150‑1 (MQ=255)
                             tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCACGTGAGGGCTGGGTAGAGTCAAGCACCCGGAAAGGGTGGATCCTG‑CGGGGCGCGCCTTGGATGGCATGGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGAc                                                                                                              >  2:1005169/1‑150 (MQ=255)
                               tGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGa                                                                                                            <  2:1239629/150‑1 (MQ=255)
                                 aCGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGc                                                                                                          <  2:609285/150‑1 (MQ=255)
                                  cGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCa                                                                                                         >  2:698076/1‑150 (MQ=255)
                                     cAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                                      >  2:1613230/1‑150 (MQ=255)
                                        gTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCg                                                                                                   <  2:1289758/150‑1 (MQ=255)
                                        gTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCg                                                                                                   >  2:1131783/1‑150 (MQ=255)
                                         aaCATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCgg                                                                                                  <  2:393445/149‑1 (MQ=255)
                                            aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCg                                                                                                   >  1:676025/1‑146 (MQ=255)
                                            aTCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCg                                                                                                   <  2:676025/146‑1 (MQ=255)
                                             tCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                              <  1:1042736/150‑1 (MQ=255)
                                             tCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                              <  2:961747/150‑1 (MQ=255)
                                               gTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                            >  1:216006/1‑150 (MQ=255)
                                               gTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                            <  1:698076/150‑1 (MQ=255)
                                                   aCGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctgc                                                                                        >  2:1148277/1‑150 (MQ=255)
                                                    cGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                              <  2:222167/143‑1 (MQ=255)
                                                    cGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                              >  1:222167/1‑143 (MQ=255)
                                                     gATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctg                                                                                         <  1:840201/147‑1 (MQ=255)
                                                     gATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctg                                                                                         >  2:840201/1‑147 (MQ=255)
                                                       tACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgggg                                                                                                                                    >  2:30480/1‑102 (MQ=255)
                                                       tACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGa                                                                                    >  2:132224/1‑150 (MQ=255)
                                                       tACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTCATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgggg                                                                                                                                    <  1:30480/102‑1 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTg                                                                                                                                                    >  1:1497052/1‑83 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTg                                                                                                                                                    <  2:1497052/83‑1 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGa                                                                                                                      <  2:176679/113‑1 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGa                                                                                                                      >  1:176679/1‑113 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                 >  2:146930/1‑150 (MQ=255)
                                                          ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                 <  1:1131783/150‑1 (MQ=255)
                                                           cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTTCGATGCCAATGGCATCgcag                                                                                                                                           >  2:513183/1‑88 (MQ=39)
                                                           cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTTCGATGCCAATGGCATCGCAg                                                                                                                                           <  1:513183/91‑1 (MQ=39)
                                                            aGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCa                                                                               >  1:263437/1‑150 (MQ=255)
                                                                  gCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCg                                                                         >  1:285975/1‑150 (MQ=255)
                                                                         gTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                            >  1:495169/1‑124 (MQ=255)
                                                                         gTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                            <  2:495169/124‑1 (MQ=255)
                                                                         gTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAgg                                                                  <  1:52859/150‑1 (MQ=255)
                                                                           gCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCa                                                                <  1:1148277/150‑1 (MQ=255)
                                                                             ggagtGGATAGCGTCAAGGAACAGGAACTGATATA‑CCTGACCGCGGGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTCATGCAGGACGAGCAGCGCCGGTTGTGGCTGCTCCAACCCCACACCGCACCAGGCacg                                                              >  2:104223/5‑148 (MQ=255)
                                                                                    aTAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAgg                                                                                                       <  2:1301920/102‑1 (MQ=255)
                                                                                    aTAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAgg                                                                                                       >  1:1301920/1‑102 (MQ=255)
                                                                                         gTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATAACAGATCATAGc                                                  <  2:263437/150‑1 (MQ=255)
                                                                                          tCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGa                                                             >  1:1400397/1‑138 (MQ=255)
                                                                                          tCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGa                                                             <  2:1400397/138‑1 (MQ=255)
                                                                                           cAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCt                                                <  2:995553/150‑1 (MQ=255)
                                                                                                 acaggaacTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGcc                                          >  2:951828/1‑150 (MQ=255)
                                                                                                   aggaacTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCa                                        >  1:825738/1‑150 (MQ=255)
                                                                                                     gaacTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCagg                                      <  2:216006/150‑3 (MQ=255)
                                                                                                     gaacTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATg                                      >  2:305026/1‑150 (MQ=255)
                                                                                                     gaacTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATg                                      <  1:132224/150‑1 (MQ=255)
                                                                                                      aacTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATgg                                     <  1:951828/150‑1 (MQ=255)
                                                                                                       acTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATggc                                    <  1:1535393/150‑1 (MQ=255)
                                                                                                           aTCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATg                                      >  2:1131278/1‑144 (MQ=255)
                                                                                                           aTCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATg                                      <  1:1131278/144‑1 (MQ=255)
                                                                                                               a‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTg                            >  1:1610548/1‑150 (MQ=255)
                                                                                                                        gcgcgcgcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGc                    <  2:285975/150‑1 (MQ=255)
                                                                                                                          gcgcgcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCt                  <  1:1005169/150‑1 (MQ=255)
                                                                                                                                cTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt                                                               >  1:1475985/1‑99 (MQ=255)
                                                                                                                                cTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt                                                               <  2:1475985/99‑1 (MQ=255)
                                                                                                                                 tGCGATGCCATTGGCATCGAATGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCt                                                <  1:17183/113‑1 (MQ=255)
                                                                                                                                 tGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCt                                                >  2:17183/1‑113 (MQ=255)
                                                                                                                                 tGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTc           >  1:351514/1‑150 (MQ=255)
                                                                                                                                   cGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCa         <  1:1390001/150‑1 (MQ=255)
                                                                                                                                          aTTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTg  <  1:305026/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGA‑CCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTG  >  NC_002947/1526484‑1526770

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: ‑