Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,436,937 A→T 100% D305E (GAT→GAA pdxB ← erythronate‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,436,9370AT99.0% 314.4 / ‑4.3 95D305E (GAT→GAApdxBerythronate‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (37/57);  total (38/57)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

TATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTAT  >  NC_000913/2436840‑2437031
                                                                                                 |                                                                                              
tATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGTGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTc                                                                                              <  2:2125803/100‑1 (MQ=255)
tATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTc                                                                                              <  1:764959/100‑1 (MQ=255)
 aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                             <  2:1742184/100‑1 (MQ=255)
 aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                             <  1:740268/100‑1 (MQ=255)
  tCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGa                                                                                            >  1:960157/1‑100 (MQ=255)
   cGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAg                                                                                           >  1:59062/1‑100 (MQ=255)
     aaCTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCg                                                                                         >  1:199616/1‑100 (MQ=255)
     aaCTCACCCGGTATCCCGGCGACTTTACGCAGCGGTCCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGc                                                                                          >  1:678524/1‑99 (MQ=255)
      aCTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCgg                                                                                        <  2:322302/100‑1 (MQ=255)
        tCACCCGGTATCCCGGCGTCTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                      <  1:1487297/100‑1 (MQ=255)
        tCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                      <  1:1155406/100‑1 (MQ=255)
        tCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                      <  1:88140/100‑1 (MQ=255)
          aCCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCa                                                                                    <  1:456820/100‑1 (MQ=255)
           cccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCTCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCAt                                                                                   <  1:51773/100‑1 (MQ=255)
            ccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATg                                                                                  >  1:2400732/1‑100 (MQ=255)
            ccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATg                                                                                  >  2:2491149/1‑100 (MQ=255)
            ccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATg                                                                                  <  1:1356937/100‑1 (MQ=255)
             cGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGc                                                                                 <  1:2812944/100‑1 (MQ=255)
             cGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGc                                                                                 <  2:1976672/100‑1 (MQ=255)
                tATCCCGGCGACTTTACGCATCGGTGCGTCAGCGCTGCGCACATCATACACCAAATGCTCCAGCCTTTTCAGCGTAGGTTGTTCGAGCGGGCCATGCAg                                                                               <  2:2928188/99‑1 (MQ=255)
                    ccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTaa                                                                          >  2:887679/1‑100 (MQ=255)
                     cGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAAt                                                                         <  2:399342/100‑1 (MQ=255)
                      ggCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATg                                                                        >  2:688406/1‑100 (MQ=255)
                      ggCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATg                                                                        <  2:1754680/100‑1 (MQ=255)
                      ggCGACTTTACGCAGCGGGGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATg                                                                        <  2:2066816/100‑1 (MQ=255)
                       gctaCTTTACGCATCGGTGCGTCATCGCGGCGCACATCATACACCATATGCACCTGTCTTTTCAGCGTCGGTTGTTTGAGCGGGCCATGCAGCGTAATGc                                                                       <  2:2787102/97‑1 (MQ=255)
                       gCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGc                                                                       <  1:1632705/100‑1 (MQ=255)
                          aCTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAc                                                                    <  2:1329498/100‑1 (MQ=255)
                           cTTTCCCCAGCGTAGCGTCATCGCGACTCACATCATACCCCAAATGCACCATTTTTTTCAGCGTCGGTTGTTCGATCTTGCCATGCAGCTTAATGCGAc                                                                    <  2:2577739/99‑1 (MQ=255)
                           cTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAcc                                                                   <  2:763364/100‑1 (MQ=255)
                           cTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAcc                                                                   >  1:1235052/1‑100 (MQ=255)
                           cTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAcc                                                                   <  2:454100/100‑1 (MQ=255)
                            tttACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCa                                                                  <  2:349127/100‑1 (MQ=255)
                              tACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCaaa                                                                >  1:1317859/1‑100 (MQ=255)
                               aCGCAGCGGTGCGTCATCGTGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAc                                                               <  2:304644/100‑1 (MQ=255)
                                cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCTTAATGCGACCAAAct                                                              >  1:21855/1‑100 (MQ=255)
                                cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAct                                                              >  1:2301506/1‑100 (MQ=255)
                                  cAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAActct                                                            >  2:1419414/1‑100 (MQ=255)
                                  cAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAActct                                                            >  1:1271150/1‑100 (MQ=255)
                                   aGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTg                                                           <  2:3020328/100‑1 (MQ=255)
                                    gCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTgg                                                          >  1:2514617/1‑100 (MQ=255)
                                    gCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTgg                                                          <  2:938285/100‑1 (MQ=255)
                                     cGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgc                                                         <  1:1913782/100‑1 (MQ=255)
                                      ggTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcg                                                        >  1:1558290/1‑100 (MQ=255)
                                         gCGTCATCGCGGCGCACATCATACACCAAATGCACCCGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                     <  2:2846918/100‑1 (MQ=255)
                                         gCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                     <  2:1235040/100‑1 (MQ=255)
                                         gCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                     <  2:1732813/100‑1 (MQ=255)
                                          cGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcagg                                                    >  2:2518867/1‑100 (MQ=255)
                                           gTCATCGCGGCGCACATCTTACACCAAGTTCGCCAGTCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCACTGTAATGCGACTAAACTCGGGCgaagcc                                                   <  2:1888238/100‑3 (MQ=255)
                                           gTCATCGCGGCGCACATCATACACCAAATGCACCATCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggc                                                   <  1:1809149/100‑1 (MQ=255)
                                           gTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggc                                                   >  2:457853/1‑100 (MQ=255)
                                             cATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggcag                                                 >  2:2291692/1‑100 (MQ=255)
                                              aTCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAACGGGCCATACAGCATAATGCGACCAAACTCTGGCgcaggcag                                                 >  1:1996904/1‑99 (MQ=255)
                                               tCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTa                                               <  2:963621/100‑1 (MQ=255)
                                                cgcgGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTaa                                              <  2:3043783/100‑1 (MQ=255)
                                                  cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  2:1317847/100‑1 (MQ=255)
                                                  cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  1:1724518/100‑1 (MQ=255)
                                                   ggcgCACATCATACACCAAATGCACCAGCCTTTTCAGTGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTaa                                              <  1:1887676/97‑1 (MQ=255)
                                                   ggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  2:2016169/99‑1 (MQ=255)
                                                     cgcACATCATACACCAAATGCACCATCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGTCTCAGGCAGTAATg                                            <  2:1123024/97‑1 (MQ=255)
                                                     cgcACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTAt                                         >  1:1747519/1‑100 (MQ=255)
                                                     cgcACATCATACACCAAATACACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTAt                                         >  2:1232857/1‑100 (MQ=255)
                                                       cacaTCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATcc                                       <  2:427055/100‑1 (MQ=255)
                                                        acaTCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCa                                      >  2:251690/1‑100 (MQ=255)
                                                         catcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCa                                      >  1:1448890/1‑99 (MQ=255)
                                                         catcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAg                                     >  1:1926519/1‑100 (MQ=255)
                                                            catACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgcgc                                  <  2:3019093/100‑1 (MQ=255)
                                                            catAAACCAAATGCACCAACCTTTTCAGCGTCTGTTGTTCGAGCGGGCCATGCAGCATAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgcgc                                  >  2:1482117/1‑100 (MQ=255)
                                                               acacCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAAc                               <  1:1927275/100‑1 (MQ=255)
                                                                cacCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACg                              >  2:2621099/1‑100 (MQ=255)
                                                                cacCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGAGACCAAACTCTGACGCAGGCAGTAATGTATCCAGCGCAACg                              >  1:305005/1‑100 (MQ=255)
                                                                  ccAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTg                            >  2:2810459/1‑100 (MQ=255)
                                                                   cAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGc                           <  2:1554331/100‑1 (MQ=255)
                                                                          agcAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCAt                    >  1:2045753/3‑100 (MQ=255)
                                                                          aCCAGGCTTTTCAGCGTCGGTTGTTCTGGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCATCGCAACGTGCTGTTCa                     <  2:2831166/99‑1 (MQ=255)
                                                                             aGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGcc                 >  1:2737840/1‑100 (MQ=255)
                                                                             aGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGcc                 <  2:1781653/100‑1 (MQ=255)
                                                                              gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                <  1:2985894/100‑1 (MQ=255)
                                                                              gCCTTTTCAGCGTCGCTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                >  2:356286/1‑100 (MQ=255)
                                                                                cTTTTCAGCGTCGGTTGTTCGAGCAGGCCATGCAGCGTAATGCGACCAAACTCTGGCGAAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCaa              >  2:1141311/1‑100 (MQ=255)
                                                                                  tttCAGTGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATTCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATa            <  1:1585272/100‑1 (MQ=255)
                                                                                    tCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATaaa          <  1:190600/100‑1 (MQ=255)
                                                                                    tCAGCGTAGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTCATGTATCCAGCGCAACGTGCTGTTCATGCCCAATaaa          <  2:1623503/100‑1 (MQ=255)
                                                                                     cAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAAc         >  2:2778202/1‑100 (MQ=255)
                                                                                      aGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGTGCAACGTGCTGTTCATGCCCAATAAAc         <  2:1000694/99‑1 (MQ=255)
                                                                                      aGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACt        <  1:1206260/100‑1 (MQ=255)
                                                                                      aGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACt        >  2:1469779/1‑100 (MQ=255)
                                                                                        cGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTg      <  1:2307057/100‑1 (MQ=255)
                                                                                        cGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTg      <  2:2296099/100‑1 (MQ=255)
                                                                                        cGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTg      >  1:1646554/1‑100 (MQ=255)
                                                                                         gTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGc     <  1:229476/100‑1 (MQ=255)
                                                                                         gTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGc     >  2:1928011/1‑100 (MQ=255)
                                                                                         gTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGc     <  2:159219/100‑1 (MQ=255)
                                                                                          tCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCt    <  2:1667051/100‑1 (MQ=255)
                                                                                            ggTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCtat  >  2:1528822/1‑100 (MQ=255)
                                                                                                 |                                                                                              
TATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTAT  >  NC_000913/2436840‑2437031

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