Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 87,887 C→T 8.5% intergenic (+39/‑141) ilvH → / → cra acetolactate synthase/acetohydroxybutanoate synthase, regulatory subunit/DNA‑binding transcriptional dual regulator Cra

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091387,8870CT 8.5% 244.3 / 8.9 82intergenic (+39/‑141)ilvH/craacetolactate synthase/acetohydroxybutanoate synthase, regulatory subunit/DNA‑binding transcriptional dual regulator Cra
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base T (2/5);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

AAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGTGAAAC  >  NC_000913/87744‑88034
                                                                                                                                               |                                                                                                                                                   
aaGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTggg                                                                                                                                               >  2:1506200/1‑150 (MQ=255)
 aGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGc                                                                                                                                              <  1:737942/150‑1 (MQ=255)
  gCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGCt                                                                                                                                             <  2:883452/150‑1 (MQ=255)
  gCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCt                                                                                                                                             >  1:1103216/1‑150 (MQ=255)
  gCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCt                                                                                                                                             <  2:299977/150‑1 (MQ=255)
  gCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCt                                                                                                                                             >  1:1665074/1‑150 (MQ=255)
      gATGCATTTTTAGCATCGATTCGCTATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGCtgtaa                                                                                                                                         <  2:1953836/150‑5 (MQ=255)
        tGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGCtgtaaat                                                                                                                                       <  1:1880193/150‑7 (MQ=255)
         gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGCtgtaaatt                                                                                                                                      <  2:680893/150‑8 (MQ=255)
          cATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGc                                                                                                                                     >  1:630800/1‑150 (MQ=255)
                tAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAATTc                                                                                                                               <  2:561820/150‑1 (MQ=255)
                   cATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                            <  1:1893945/150‑1 (MQ=255)
                    aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTg                                                                                                                           <  1:1506200/150‑1 (MQ=255)
                     tCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTgg                                                                                                                          >  1:440353/1‑150 (MQ=255)
                       gATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGa                                                                                                                        >  2:1875628/1‑150 (MQ=255)
                                tgtgGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAAt                                                                                                               <  2:833323/150‑1 (MQ=255)
                                tgtgGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAAt                                                                                                               <  1:1000887/150‑1 (MQ=255)
                                   ggCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaa                                                                                                            <  2:205174/150‑1 (MQ=255)
                                    gCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaaa                                                                                                           >  1:142685/1‑150 (MQ=255)
                                             gtgGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGTCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGcc                                                                                                  >  2:1603571/1‑150 (MQ=255)
                                              tgGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCAGGAATAAAAGCAGTTGCCg                                                                                                 >  2:1068326/1‑150 (MQ=255)
                                               ggAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGc                                                                                                >  1:1666153/1‑150 (MQ=255)
                                                gAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTATTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCa                                                                                               >  2:894250/1‑150 (MQ=255)
                                                 aGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAg                                                                                              <  2:1945986/150‑1 (MQ=255)
                                                    ttGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGTCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTa                                                                                           <  1:1603571/150‑1 (MQ=255)
                                                      gctcgctcTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAAt                                                                                         >  1:110867/1‑150 (MQ=255)
                                                       ctcgctcTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAAtt                                                                                        >  1:610745/1‑150 (MQ=255)
                                                           ctctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATATTCt                                                                                    >  1:1255923/1‑150 (MQ=255)
                                                             ctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATTATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTgc                                                                                  <  2:1665074/150‑1 (MQ=255)
                                                                  gtgGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTa                                                                             >  2:787687/1‑150 (MQ=255)
                                                                   tgGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAg                                                                            >  1:956606/1‑150 (MQ=255)
                                                                   tgGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAg                                                                            <  1:776197/150‑1 (MQ=255)
                                                                    ggTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGa                                                                           >  1:1437363/1‑150 (MQ=255)
                                                                    ggTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGa                                                                           >  2:1236132/1‑150 (MQ=255)
                                                                     gTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                          >  1:1333799/1‑150 (MQ=255)
                                                                     gTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTATTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                          <  1:894250/150‑1 (MQ=255)
                                                                       cGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGt                                                                        >  1:998227/1‑150 (MQ=255)
                                                                        ggACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGtt                                                                       <  2:29620/150‑1 (MQ=255)
                                                                        ggACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGtt                                                                       <  2:301979/150‑1 (MQ=255)
                                                                            actcccgccgcGAAAAAAAAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGACACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTACTGCGCTTAGATGTTAATg                                                                   <  2:1260905/145‑1 (MQ=255)
                                                                             ttCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGa                                                                  <  2:440353/150‑1 (MQ=255)
                                                                                 cgcgGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAAttt                                                              >  1:1945475/1‑150 (MQ=255)
                                                                                     gcgATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAAcc                                                          <  2:998227/150‑1 (MQ=255)
                                                                                       gATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCa                                                        >  1:490776/1‑150 (MQ=255)
                                                                                       caTAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTATTTGAGAAATCAGTGGGAACCTGGAAAAAAAGCAGTTGCCGCAGTTAATTTTATGCGCTTAGATGTTAAAGAATATAACCCa                                                        >  1:1026930/2‑150 (MQ=255)
                                                                                         tAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATa                                                      <  1:471528/150‑1 (MQ=255)
                                                                                           aaaTAATGCGTTGAGTATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATAcc                                                    <  2:177982/150‑1 (MQ=255)
                                                                                           aaaTAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATAcc                                                    <  2:1103216/150‑1 (MQ=255)
                                                                                             ataatGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAg                                                  <  1:1682456/150‑1 (MQ=255)
                                                                                              taatGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGt                                                 >  2:125539/1‑150 (MQ=255)
                                                                                               aatGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTa                                                <  1:125539/150‑1 (MQ=255)
                                                                                                 tGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACa                                              >  2:401343/1‑150 (MQ=255)
                                                                                                   cGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAAt                                            >  1:315101/1‑150 (MQ=255)
                                                                                                   cGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAAt                                            >  2:886711/1‑150 (MQ=255)
                                                                                                   cGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAAt                                            >  2:1101626/1‑150 (MQ=255)
                                                                                                    gTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGTAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATg                                                                   <  1:73257/126‑1 (MQ=255)
                                                                                                    gTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATg                                                                   >  2:73257/1‑126 (MQ=255)
                                                                                                           aTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATgg                                    >  1:1993222/1‑150 (MQ=255)
                                                                                                           aTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATgg                                    >  2:1869499/1‑150 (MQ=255)
                                                                                                            tGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGt                                   <  1:886711/150‑1 (MQ=255)
                                                                                                             gATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGtt                                  >  2:1414496/1‑150 (MQ=255)
                                                                                                                ctcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttttt                               <  1:401343/150‑1 (MQ=255)
                                                                                                                ctcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttttt                               <  2:732742/150‑1 (MQ=255)
                                                                                                                ctcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttt                                 >  2:1938538/1‑148 (MQ=255)
                                                                                                                ctcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttt                                 <  1:1938538/148‑1 (MQ=255)
                                                                                                                 tcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTa                              <  2:110867/150‑1 (MQ=255)
                                                                                                                 tcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTa                              <  2:630800/150‑1 (MQ=255)
                                                                                                                     tGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACAtt                          <  1:968987/150‑1 (MQ=255)
                                                                                                                        gcAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTa                       <  1:479456/150‑1 (MQ=255)
                                                                                                                        gcAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTa                       >  2:1717925/1‑150 (MQ=255)
                                                                                                                        gcAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTa                       >  2:1159189/1‑150 (MQ=255)
                                                                                                                          caaaaaaaGCCCAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACg                     >  2:1124812/8‑150 (MQ=255)
                                                                                                                          aaTTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACg                     <  1:568098/150‑1 (MQ=255)
                                                                                                                             ttACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCaa                  >  2:748067/1‑150 (MQ=255)
                                                                                                                                cAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAggg               <  2:142685/150‑1 (MQ=255)
                                                                                                                                cAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGCTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAggg               <  2:1860612/150‑1 (MQ=255)
                                                                                                                                  gCCCAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGc             >  1:38603/1‑150 (MQ=255)
                                                                                                                                         aTGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATtgtg      <  1:1414496/150‑1 (MQ=255)
                                                                                                                                          tGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGTGa     <  2:315101/150‑1 (MQ=255)
                                                                                                                                          tGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGTGa     <  2:1666153/150‑1 (MQ=255)
                                                                                                                                             cACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGTGAAAc  >  2:695211/1‑150 (MQ=255)
                                                                                                                                             cACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGTGAAAc  >  1:1010880/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
AAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGTGAAAC  >  NC_000913/87744‑88034

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