Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,750 T→A 8.1% E198D (GAA→GAT)  dapC ← N‑succinyl‑L,L‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,780,7500TA 8.1% 168.6 / 5.1 62E198D (GAA→GAT) dapCN‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base A (3/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AGCGCTTGGACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTT  >  NC_002947/1780602‑1780893
                                                                                                                                                    |                                                                                                                                               
aGCGCTTGGACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCgtgt                                                                                                                                                 <  4:224052/149‑1 (MQ=255)
aGCGCTTGGACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCgtgt                                                                                                                                                 >  1:265778/1‑149 (MQ=255)
aGCGCTTGGACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCgtgt                                                                                                                                                 >  1:334226/1‑149 (MQ=255)
 gcgcTTGGACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGtt                                                                                                                                                >  1:347445/1‑149 (MQ=255)
  cgcTTGGACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTc                                                                                                                                               >  6:346840/1‑149 (MQ=255)
    cTTGGACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGt                                                                                                                                             >  7:109832/1‑149 (MQ=255)
        gACAGGCTGTGGAACACCGCGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGc                                                                                                                                         <  4:354476/149‑1 (MQ=255)
        gACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATCCCCACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCTGCGGCGATCACGAAGTCGTCTACGTCCGc                                                                                                                                         >  8:341245/1‑149 (MQ=255)
         aCAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGc                                                                                                                                         >  1:114333/1‑148 (MQ=255)
              cTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCgtgt                                                                                                                                                 <  2:284393/135‑1 (MQ=255)
              cTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCgtgt                                                                                                                                                 >  1:284393/1‑135 (MQ=255)
               tgtgGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGc                                                                                                                                   >  7:9941/1‑148 (MQ=255)
                   gAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCa                                                                                                                              <  1:273636/149‑1 (MQ=255)
                     acacCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAgg                                                                                                                                     >  6:64384/1‑140 (MQ=255)
                     acacCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAgg                                                                                                                                     <  5:64384/140‑1 (MQ=255)
                      caccacGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACATCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGCTCGTCAGCCAGGGCAATCAGct                                                                                                                           <  7:233713/149‑1 (MQ=255)
                        ccacGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGcttc                                                                                                                         >  6:49797/1‑149 (MQ=255)
                         cacGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGcttct                                                                                                                        <  5:346840/149‑1 (MQ=255)
                          acGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGcttctt                                                                                                                       <  5:23565/149‑1 (MQ=255)
                          acGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGcttctt                                                                                                                       <  4:43342/149‑1 (MQ=255)
                           cGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCttcttc                                                                                                                      <  4:227160/149‑1 (MQ=255)
                            gCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCttcttc                                                                                                                      <  1:315408/148‑1 (MQ=255)
                             cAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAg                                                                                                                    <  1:291427/149‑1 (MQ=255)
                                      aaaTCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCAt                                                                                                           <  1:157123/149‑1 (MQ=255)
                                       aaTCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATc                                                                                                          >  1:254875/1‑149 (MQ=255)
                                        aTCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCg                                                                                                         <  2:347445/149‑1 (MQ=255)
                                               cGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCa                                                                                                   <  2:265778/148‑1 (MQ=255)
                                                gACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAgg                                                                                                                                     <  7:95289/113‑1 (MQ=255)
                                                gACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAgg                                                                                                                                     >  8:95289/1‑113 (MQ=255)
                                                 aCCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATc                                                                                                          >  1:15141/1‑139 (MQ=255)
                                                 aCCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATc                                                                                                          <  2:15141/139‑1 (MQ=255)
                                                  ccAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAgcg                                                                                                >  2:171621/1‑148 (MQ=255)
                                                   cAAGCCCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCa                                                                                              >  4:7068/1‑149 (MQ=255)
                                                          ggCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAg                                                                                       <  3:43300/149‑1 (MQ=255)
                                                          ggCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAg                                                                                       >  7:261320/1‑149 (MQ=255)
                                                             gcATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGt                                                                                    <  6:320648/149‑1 (MQ=255)
                                                                     aCAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATc                                                                                                          <  4:79903/119‑1 (MQ=255)
                                                                         cagGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGc                                                                        <  2:114333/149‑1 (MQ=255)
                                                                            gCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGcaca                                                                     <  8:109832/149‑1 (MQ=255)
                                                                             ccTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCacag                                                                    <  2:334226/149‑1 (MQ=255)
                                                                              cTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCacag                                                                    >  4:179877/1‑148 (MQ=255)
                                                                               tGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGACAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCacagg                                                                   >  6:355873/1‑148 (MQ=255)
                                                                                  cGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGACAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCacaggaaca                                                               <  5:355873/149‑1 (MQ=255)
                                                                                       gcgcATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCacag                                                                    >  5:337873/1‑139 (MQ=255)
                                                                                       gcgcATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCacag                                                                    <  6:337873/139‑1 (MQ=255)
                                                                                          cATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTgg                                                       <  3:7068/149‑1 (MQ=255)
                                                                                          cATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTgg                                                       <  2:254875/149‑1 (MQ=255)
                                                                                              cTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAAc                                                   >  8:33005/1‑149 (MQ=255)
                                                                                                 gtcgAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCt                                                >  3:140771/1‑149 (MQ=255)
                                                                                                   cgAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTc                                              >  2:93619/1‑149 (MQ=255)
                                                                                                    gAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTc                                              <  2:205283/148‑1 (MQ=255)
                                                                                                     aaGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCca                                            <  8:9941/149‑1 (MQ=255)
                                                                                                     aaGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCATGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCca                                            >  4:14127/1‑149 (MQ=255)
                                                                                                       gTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCcaca                                          >  1:112105/1‑149 (MQ=255)
                                                                                                         aCAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACAAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACAcc                                        >  5:253670/1‑149 (MQ=255)
                                                                                                             cTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCTTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCgg                                    <  1:222199/149‑1 (MQ=255)
                                                                                                                gCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGcc                                  <  5:49797/148‑1 (MQ=255)
                                                                                                                        aCTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCg                         <  4:140771/149‑1 (MQ=255)
                                                                                                                          tcgtcgGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTc                       >  8:133604/1‑149 (MQ=255)
                                                                                                                                 tgacgAACACGACTTCGTGATCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTcc                                             <  1:63817/117‑1 (MQ=255)
                                                                                                                                 tgacgAACACGACTTCGTGATCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTcc                                             >  2:63817/4‑120 (MQ=255)
                                                                                                                                 tgacgAACACGACTTCGTGATCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTcc                >  4:24947/4‑149 (MQ=255)
                                                                                                                                   acgAACACGACTTCGTGATCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCgg              >  4:95511/2‑149 (MQ=255)
                                                                                                                                    cgATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCggg             <  8:261320/149‑1 (MQ=255)
                                                                                                                                          cgacttCGTGATCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGaaa       <  3:24947/144‑1 (MQ=255)
                                                                                                                                           gAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGACTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTACGGCTTGaaaa      >  5:26930/1‑148 (MQ=255)
                                                                                                                                               tCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCgtt  <  6:253670/149‑1 (MQ=255)
                                                                                                                                               tCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCgtt  <  2:112105/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AGCGCTTGGACAGGCTGTGGAACACCACGCAGCGCGCAAAATCGCTACGACCAAGCTCGGCGCATGCCGACAGCAGGCCTGGCGGTGGCGCATCCTCGTCGAAGTACAGCTCGCTGTAGCACTCGTCGGCGGCGATCACGAAGTCGTGTTCGTCAGCCAGGGCAATCAGCTTCTTCAGGGTGTTCATCGGCACCAGCGCACCGGTAGGGTTGCCTGGCGAGCACAGGAACAGGATCTGGCAACGCTTCCACACCTCGGCCGGTACGGCGTCGAAGTCCGGGTTGAAACCGTT  >  NC_002947/1780602‑1780893

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