Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 87,881 A→G 8.6% intergenic (+33/‑147) 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,8810AG 8.6% 153.6 / 5.6 65intergenic (+33/‑147)ilvH/craacetolactate synthase/acetohydroxybutanoate synthase, regulatory subunit/DNA‑binding transcriptional dual regulator Cra
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base G (3/3);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-02

GGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGG  >  NC_000913/87741‑88021
                                                                                                                                            |                                                                                                                                              
ggTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGCGGTTGATCTCAATGCGCAATTTACAGCCCAACTTg                                                                                                                                              >  1:112522/1‑143 (MQ=255)
 gTAAGCTTGATGCTTTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAACTCATAAGATGTGAATGATCTCAATGCGCAATTTACAGCCCAACATGt                                                                                                                                             <  1:479338/143‑1 (MQ=255)
   taGCGTGCTGCATTTTTAGCAGCGATTGGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCa                                                                                                                                           <  2:148123/142‑1 (MQ=255)
       ttGAGGCGTGTTTGGCATCGATTCGCGAGGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGtt                                                                                                                                       <  2:434508/143‑1 (MQ=255)
          aTGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAGTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTggg                                                                                                                                    >  2:909176/1‑143 (MQ=255)
          aTGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTggg                                                                                                                                    <  1:954188/143‑1 (MQ=255)
          aTGCATTGTTAGCATCGACTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCGCTGGTGTGGTCGGACTTTCGCGCGGCGATACAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTggg                                                                                                                                    <  2:560379/143‑1 (MQ=255)
           tGCCTTTTTAGCATCGATGCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGc                                                                                                                                   <  1:475784/143‑1 (MQ=255)
                  ttAGCCTCGATTCGCGATATGGCGAAAATTGTGGAGGTGGCTCGCTCTTGTGTGGTCGGACTTTCGCGCGACGATAAAATAATGCGTTGAGAATTATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttt                                                                                                                           <  1:486482/143‑1 (MQ=255)
                  tgggggtCGAGTCGGGGTGTGGCGAATATTGGGGAGGTTGCGCGCTCTGGTGGGGTCGGACTTTCGCGGGGTGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttt                                                                                                                           <  2:201231/137‑1 (MQ=255)
                       aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATCATGCTTTGAGAATGATCTGCATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTGTTTTGCGaa                                                                                                                      >  1:704979/1‑143 (MQ=255)
                        tCGATTCGCGATGTGGCGAAAGTTGTGGAGGGTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGaaa                                                                                                                     <  2:1071933/143‑1 (MQ=255)
                        tCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGTGCAATTTACAGGCGGACGTGTGATGTGGGGCTTTTTTTGCGcaa                                                                                                                     >  1:153064/1‑143 (MQ=255)
                         cGATTCGTGATGTGGCGAAAATTGTGGAGGTTGCTCGCTATGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAAt                                                                                                                    <  1:810066/143‑1 (MQ=255)
                             tCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGGGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                >  2:912941/1‑143 (MQ=255)
                              cggtaTGTTGCTACAATTGTGGAGGTTGCTCGCTCTTGTGTCGTTGGACTTTCCCGCATCGATAATATAATGCGTTGAGAATGATCTCAATTCGCAATTTACATCCCATCATGTCACGTTGGGCTTTTTTTGCGAAATCAGTg                                                                                                               <  2:27829/139‑1 (MQ=255)
                              cgcgATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGAGTTTGGAGCACCAATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTg                                                                                                               <  1:519471/143‑1 (MQ=255)
                               tggATGTGGCGAAAATTGTGGAGGTTGCTCGCTCGGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTgg                                                                                                              <  2:43507/141‑1 (MQ=255)
                                  attgggCTAAAATTGTGGATGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCTGCGATAAAATGATGTGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGaa                                                                                                           <  2:917882/139‑1 (MQ=255)
                                       gCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCACTTTACACCCCAACATCTCACCTTGGGCTTTTTTTGCGAAACCAGTGGGGCCCTgg                                                                                                      >  2:203780/1‑143 (MQ=255)
                                         gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCCGTGGGAACATGGaa                                                                                                    >  2:276577/1‑143 (MQ=255)
                                              tgggtgAGGTTGCGCGCTGTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaaa                                                                                              <  2:768006/138‑1 (MQ=255)
                                              gtgtgGAGGTTGCTCGCTGTGGTGTGGTCGGACGGTCGCGCGGCGATAAAATAAGGCGGTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaaa                                                                                              <  1:226538/142‑1 (MQ=255)
                                                      gTTGCCAGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGc                                                                                      <  1:379944/143‑1 (MQ=255)
                                                       gtgccggctcGGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAGTTTACAGCGTAGCATGTCACGTTGGGGTTTTTTTGCGTAATCAGTGGGAACCTGGAATAAAAGCAGTTGcc                                                                                     <  1:1050680/137‑1 (MQ=255)
                                                        ttctcgttcTGGTGGGGTCGGACTTTGGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACCGGCCAACATGTCACGTTGGGCTATTTTTGCTAAAGCAGTAGGAACCTGGAATAAAAGCAGTTGCCg                                                                                    <  2:845311/141‑1 (MQ=255)
                                                         gcgcgctcGGGTGTGCTCGGATTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAAGTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGCGGGAACCTGGAATAAAAGCAGTTGCCGc                                                                                   <  2:140800/140‑1 (MQ=255)
                                                          ctcgctcTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCa                                                                                  <  1:897298/143‑1 (MQ=255)
                                                               tctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGTAACCTGGAATAAAAGCAGTTGCCGCAGTTaa                                                                             >  1:226212/1‑143 (MQ=255)
                                                               tctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTaa                                                                             >  1:223925/1‑143 (MQ=255)
                                                               tctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGCCTTTTTTTGCGAAATCAGTGGGAACCTGGAATACAAGCAGTTGCCGCAGTTaa                                                                             >  2:698267/1‑143 (MQ=255)
                                                               tctGGTGTGGTCGGACTTTCGCGCGGCGAGATGATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTaa                                                                             >  2:990405/1‑143 (MQ=255)
                                                                ctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAAt                                                                            <  2:876460/143‑1 (MQ=255)
                                                                  ggTGGGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAAttt                                                                          <  1:221324/143‑1 (MQ=255)
                                                                      tgGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTg                                                                      >  1:314680/1‑143 (MQ=255)
                                                                        ggcGGACTTTCGCGCGGCGATAAAATAAAGCGATGAGAATGAGCACAATGCGCAAGTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATGTTCTgcg                                                                    >  2:455965/3‑143 (MQ=255)
                                                                         tCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCAAACATGTCACGTTGGGCTTTTTTTGCGAAATCACTGGGAACCTGGAATAAAAGCAATTGCTGCAGTTAATTTTCTgcgc                                                                   >  2:295819/1‑143 (MQ=255)
                                                                          cGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCCACATGTCCCGTTGGGCTTTGTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCGGCAGTTAATTTTCTGCGCt                                                                  >  2:130585/1‑143 (MQ=255)
                                                                              cTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGa                                                              <  1:745433/143‑1 (MQ=255)
                                                                                ttCGCGCGGCGATAAAATAATGCGTTGAGAATGTTATCAATGCGCGATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCAGGAGTTAATTTTCTACGCTAAGATg                                                            >  2:1078576/1‑143 (MQ=255)
                                                                                 tCGCGCGGCGATAAAATAATACGTTGTGAGTGATCTCAATGCGCAATTGACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGt                                                           <  1:417399/143‑1 (MQ=255)
                                                                                     gcggccATAAAATAATGCGTTGAGAATGATCCCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAAt                                                       <  1:953094/143‑1 (MQ=255)
                                                                                       ggcgATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGCTGTTAATGaa                                                    >  1:649770/1‑143 (MQ=255)
                                                                                       agcgATAAAATAATGCGTTCAGAAGGATCTCAATGCGCAATTGACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGa                                                     <  1:563934/142‑1 (MQ=255)
                                                                                          gttAAAATAATGCGTTGGGGATGATCTCAAGGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAAtt                                                  <  2:158773/141‑1 (MQ=255)
                                                                                          gATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGTGGCCATTAATTTTCTGCGCTTAGATGTTAATGAAtt                                                  <  2:1110356/143‑1 (MQ=255)
                                                                                            tAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTa                                                >  2:214844/1‑143 (MQ=255)
                                                                                               aatgatGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAAcc                                             <  1:409094/143‑1 (MQ=255)
                                                                                                        ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGt                                    >  1:1032236/1‑143 (MQ=255)
                                                                                                        ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCAATACCAGt                                    >  1:143236/1‑143 (MQ=255)
                                                                                                        ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCGTAGATGTTGATGAATTTAACCTATACCAGt                                    >  1:876393/1‑143 (MQ=255)
                                                                                                         tGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAACAAAACCAGTTGCCGCCATTAATTTTCAGCGCTTCGATATCAATTAATGTAACCCTTACCAATa                                   >  2:849392/1‑143 (MQ=255)
                                                                                                          gagaATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTAc                                  >  2:1112004/1‑143 (MQ=255)
                                                                                                                 aTCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCACTACAATGGcc                           >  2:950370/1‑142 (MQ=255)
                                                                                                                    tcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTAt                         <  2:212072/143‑1 (MQ=255)
                                                                                                                      aaTGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCCCTTAGATGTTAATGAATTTAACCCATACCAGTACAAGTGCAATGGt                      >  2:51057/1‑143 (MQ=255)
                                                                                                                        tGCGCAATTTACAGCCGAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTAACCGCAGGTATTGTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttt                    <  2:187628/143‑1 (MQ=255)
                                                                                                                            cAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTa                 >  1:99743/1‑143 (MQ=255)
                                                                                                                               tttACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGGCGGAACGTTGAGAATGATGAGTTGCCGCAGTGAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACAtt             <  1:258077/143‑1 (MQ=255)
                                                                                                                               tggagaGCCGAAGATGTAACGTTGGGCTTGTTCTGAGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTGAATTTGCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACAtt             <  2:769774/138‑1 (MQ=255)
                                                                                                                                  aCAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTa          <  1:264877/143‑1 (MQ=255)
                                                                                                                                    aGCCCAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAAGTTTCTGCGCTGAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACg        >  1:685111/1‑143 (MQ=255)
                                                                                                                                    aGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGTATAAAAGCAGATACCGCAGTTAATTTTATGAGATTAGATGATAATGAATTTAACCCATACTAGTACAATGGCTATGGTATTTACATTTTACg        >  1:661123/1‑143 (MQ=255)
                                                                                                                                      cccAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCa      >  2:881198/1‑143 (MQ=255)
                                                                                                                                          aCATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAggg  >  1:541758/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                              
GGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGG  >  NC_000913/87741‑88021

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