Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 122,070 A→C intergenic (‑519/‑22) aroP ← / → pdhR aromatic amino acid transporter/pyruvate dehydrogenase complex repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913122,0700AC95.5% 584.5 / 12.3 179intergenic (‑519/‑22)aroP/pdhRaromatic amino acid transporter/pyruvate dehydrogenase complex repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (6/2);  new base C (88/82);  total (94/85)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 179 total aligned reads displayed.

TTCATGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGAA  >  NC_000913/121974‑122166
                                                                                                |                                                                                                
ttCATGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGt                                                                                               >  1:404828/1‑100 (MQ=255)
ttCATGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGt                                                                                               >  1:773741/1‑100 (MQ=255)
 tCATGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTg                                                                                              >  1:2734126/1‑100 (MQ=255)
 tCATGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTg                                                                                              <  1:3843195/100‑1 (MQ=255)
  cATGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTgg                                                                                             >  1:3106067/1‑100 (MQ=255)
   aTGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTTGCGTCACAGACATGAAATTGCTAAGACGAATTGACTTCGGCCAGTGGc                                                                                            <  2:3984295/100‑1 (MQ=255)
   aTGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGc                                                                                            >  1:3445075/1‑100 (MQ=255)
   aTGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGc                                                                                            >  1:2574939/1‑100 (MQ=255)
     gATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATTAAATTGGTAATACCGATTGACTTCGGCCAGTGGCtt                                                                                          >  2:3784451/1‑100 (MQ=255)
     gATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCtt                                                                                          <  1:1935157/100‑1 (MQ=255)
      aTTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTa                                                                                         >  1:1768966/1‑100 (MQ=255)
        ttCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTaan                                                                                       <  2:1515015/100‑2 (MQ=255)
        ttCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAg                                                                                       >  2:1035398/1‑100 (MQ=255)
         tcaatcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGa                                                                                      >  1:419786/1‑100 (MQ=255)
           aatcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACa                                                                                    >  2:1476877/1‑100 (MQ=255)
           aatcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCAAGTGGCTTAAGACa                                                                                    >  1:839596/1‑100 (MQ=255)
            atcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAg                                                                                   >  2:549840/1‑100 (MQ=255)
            atcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCAAGTGGCTTAAGACAg                                                                                   >  1:491557/1‑100 (MQ=255)
            atcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATGGGTAAGACCAATTGACTTCGGCCATTGGCTTAAGAAAg                                                                                   <  2:2051567/100‑1 (MQ=255)
             tcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAgg                                                                                  >  2:2080427/1‑100 (MQ=255)
             tcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAgg                                                                                  <  1:3117329/100‑1 (MQ=255)
             tcaaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGACCAGTGGCTTAAGACAgg                                                                                  >  2:3137871/1‑100 (MQ=255)
              caaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGa                                                                                 >  2:2874593/1‑100 (MQ=255)
              caaAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGa                                                                                 >  1:2985176/1‑100 (MQ=255)
                   ccTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGCCAGGCACTCa                                                                            >  2:2398707/1‑100 (MQ=255)
                    cTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCAt                                                                           >  1:3626747/1‑100 (MQ=255)
                     tGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATg                                                                          >  2:2465484/1‑100 (MQ=255)
                     tGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATg                                                                          >  1:1023728/1‑100 (MQ=255)
                     tGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATg                                                                          <  1:1521074/100‑1 (MQ=255)
                      gTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATgg                                                                         <  1:1019064/100‑1 (MQ=255)
                        aTGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGcc                                                                       >  2:2318158/1‑100 (MQ=255)
                        aTGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGcc                                                                       <  2:1002229/100‑1 (MQ=255)
                         tGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCt                                                                      <  1:2337241/100‑1 (MQ=255)
                            aCATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACa                                                                   <  1:3733546/100‑1 (MQ=255)
                             cATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAg                                                                  <  1:2857896/100‑1 (MQ=255)
                             cATAAGGTGAATACTTTGTTACTTTAGCGTCAAAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAg                                                                  >  2:3597131/1‑100 (MQ=255)
                               tAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCa                                                                >  1:2193164/1‑100 (MQ=255)
                               tAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCa                                                                <  1:3340145/100‑1 (MQ=255)
                                aaGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCaa                                                               <  2:592569/100‑1 (MQ=255)
                                aaGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCaa                                                               <  2:3857406/100‑1 (MQ=255)
                                 aGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCaaa                                                              >  2:3427812/1‑100 (MQ=255)
                                        tACTTTGTTACTTTAGCGTCACAGACACGAAATTTGTAAGTCCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATTGCCTACAGCAAAATCCGcc                                                       <  2:3822814/100‑1 (MQ=255)
                                          cTTTGTTACTTTAGCGTCACAGCCATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGccaa                                                     >  1:3415308/1‑100 (MQ=255)
                                          cTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGccaa                                                     >  2:2008038/1‑100 (MQ=255)
                                            ttGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGccaacc                                                   >  1:2927168/1‑100 (MQ=255)
                                              gTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGccaaccaa                                                 >  2:3069246/1‑100 (MQ=255)
                                                 aCTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAAc                                              >  1:2912765/1‑100 (MQ=255)
                                                 aCTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAAc                                              <  2:1547142/100‑1 (MQ=255)
                                                 aCTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAAc                                              <  2:1934145/100‑1 (MQ=255)
                                                  cTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAAct                                             >  2:2655381/1‑100 (MQ=255)
                                                      aGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTcc                                         >  1:3030577/1‑100 (MQ=255)
                                                       gCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCg                                        <  1:1401618/100‑1 (MQ=255)
                                                        cGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGa                                       <  2:2503836/100‑1 (MQ=255)
                                                         gTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAt                                      >  1:1817125/1‑100 (MQ=255)
                                                         gTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAt                                      >  2:2801857/1‑100 (MQ=255)
                                                         gTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAt                                      >  1:2849844/1‑100 (MQ=255)
                                                         gTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAt                                      >  1:2792873/1‑100 (MQ=255)
                                                           cacaGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCACAATCCGCCAACCAAAACTCCCCGAtgt                                    >  2:696831/1‑100 (MQ=255)
                                                           cacaGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAtgt                                    >  1:2458589/1‑100 (MQ=255)
                                                            acaGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGAtgtg                                   >  2:2535646/1‑100 (MQ=255)
                                                             caGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGa                                  >  1:104422/1‑100 (MQ=255)
                                                              aGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAACACTCTCCGATGTGAt                                 >  1:1626695/1‑100 (MQ=255)
                                                              aGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGAt                                 <  1:1291929/100‑1 (MQ=255)
                                                              aGACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGAt                                 <  2:3724388/100‑1 (MQ=255)
                                                               gACATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGAtt                                >  2:3610408/1‑100 (MQ=255)
                                                                 cATGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGa                              <  1:3681031/100‑1 (MQ=255)
                                                                  aTGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGag                             <  1:1541245/100‑1 (MQ=255)
                                                                  aTGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGag                             <  2:3566798/100‑1 (MQ=255)
                                                                  aTGAAATTGGTAAGACCAATTGACATCGGCGAGTGGCTTAAGACATGAACGCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTCACTGag                             <  2:3768388/100‑1 (MQ=255)
                                                                   tGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGagc                            <  2:1838332/100‑1 (MQ=255)
                                                                   tGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGagc                            <  1:1502836/100‑1 (MQ=255)
                                                                   tGAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGagc                            <  1:2659681/100‑1 (MQ=255)
                                                                    gAAATTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGagca                           <  1:2345664/100‑1 (MQ=255)
                                                                     aaaTTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGagcag                          <  1:1313334/100‑1 (MQ=255)
                                                                     aaaTTGGTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGagcag                          >  1:3839608/1‑100 (MQ=255)
                                                                      aaTTGGTAAGACCGATTGACTTCGGCCGGTGGCTTAAGACAGGAACTCATGGCCTACAGGAAAATCCGCCAACCAAAACTCTCCGATGAGATTGagcagc                         >  1:3593581/1‑100 (MQ=255)
                                                                          ggTAAGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACt                     >  1:1252302/1‑100 (MQ=255)
                                                                           gTAAGACCAATTGACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTg                    >  2:2751281/1‑100 (MQ=255)
                                                                            tAAGACCAATTGACTTCGGACAGTGGCTTAAGACAGGAACTCATGGCCTACAGTAAAATCCCCCAACGAAAACTCTCCGATGTGATTGAGCAGCAACTgg                   <  1:3617729/100‑1 (MQ=255)
                                                                             aaGACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGa                  >  2:3250579/1‑100 (MQ=255)
                                                                               naCCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGt                >  1:3230951/2‑100 (MQ=255)
                                                                               gACCAATTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGt                >  1:1191909/1‑100 (MQ=255)
                                                                                   aaTTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGttttt            <  2:1973320/100‑1 (MQ=255)
                                                                                   aaTTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGttttt            <  1:2503984/100‑1 (MQ=255)
                                                                                   aaTTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGttttt            <  2:2985546/100‑1 (MQ=255)
                                                                                   aaTTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGttttt            <  2:1095001/100‑1 (MQ=255)
                                                                                    aTTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTg           <  1:336810/100‑1 (MQ=255)
                                                                                    aTTGACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTg           <  1:2946975/100‑1 (MQ=255)
                                                                                       gACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATc        <  2:3173363/100‑1 (MQ=255)
                                                                                       gACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATc        >  2:1146633/1‑100 (MQ=255)
                                                                                       gACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGGTGTTGATc        >  1:188125/1‑100 (MQ=255)
                                                                                       gACTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGAGGTCTTTCATc        >  2:1146473/1‑100 (MQ=255)
                                                                                        aCTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATcc       <  1:3519333/100‑1 (MQ=255)
                                                                                        aCTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATcc       <  2:1279978/100‑1 (MQ=255)
                                                                                         cTTCGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGCTTTTGATCCt      >  2:3796341/1‑100 (MQ=255)
                                                                                            cGGCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGa   >  2:814991/1‑100 (MQ=255)
                                                                                             ggCCAGTGGCTTAAGACAGGACCTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGCTGTGATTGAGCAGCACCTGGAGTTTTTGATCCTCGaa  >  2:1784716/1‑100 (MQ=255)
                                                                                             ggCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGaa  <  1:3128900/100‑1 (MQ=255)
                                                                                             ggCCAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGaa  >  2:916668/1‑100 (MQ=255)
                                                                                             ggCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGaa  >  1:2532387/1‑100 (MQ=255)
                                                                                                |                                                                                                
TTCATGATTTCAATCAAAACCTGTATGGACATAAGGTGAATACTTTGTTACTTTAGCGTCACAGACATGAAATTGGTAAGACCAATTGACTTCGGCAAGTGGCTTAAGACAGGAACTCATGGCCTACAGCAAAATCCGCCAACCAAAACTCTCCGATGTGATTGAGCAGCAACTGGAGTTTTTGATCCTCGAA  >  NC_000913/121974‑122166

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