Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 86,547 T→A D306E (GAT→GAA ilvI → acetolactate synthase/acetohydroxybutanoate synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091386,5470TA99.7% 1209.1 / NA 282D306E (GAT→GAAilvIacetolactate synthase/acetohydroxybutanoate synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (150/131);  minor base C (0/1);  total (150/132)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

GATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAAT  >  NC_000913/86449‑86644
                                                                                                  |                                                                                                 
 aTGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCt                                                                                                >  1:131515/1‑101 (MQ=255)
    tgATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTAc                                                                                              <  1:131316/100‑1 (MQ=255)
    tgATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACt                                                                                             >  1:131523/1‑101 (MQ=255)
    tgATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACt                                                                                             >  2:131522/1‑101 (MQ=255)
      aTTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTc                                                                                           <  2:131081/101‑1 (MQ=255)
          tCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCAtt                                                                                       <  1:131002/101‑1 (MQ=255)
          tCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCAtt                                                                                       <  1:131152/101‑1 (MQ=255)
           cGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCAttt                                                                                      >  2:131528/1‑101 (MQ=255)
           cGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCAttt                                                                                      <  1:131093/101‑1 (MQ=255)
            gCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTc                                                                                     >  1:131529/1‑101 (MQ=255)
             ccGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCt                                                                                    >  2:131533/1‑101 (MQ=255)
              cGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTa                                                                                   <  2:131195/101‑1 (MQ=255)
               gTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTaa                                                                                  >  1:131535/1‑101 (MQ=255)
               gTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTaa                                                                                  <  2:131207/101‑1 (MQ=255)
                tCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTaaa                                                                                 >  2:131538/1‑101 (MQ=255)
                 cGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTaaaa                                                                                >  2:131541/1‑101 (MQ=255)
                  ggggTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAAc                                                                               >  2:131546/1‑101 (MQ=255)
                  ggggTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAAc                                                                               >  1:131544/1‑101 (MQ=255)
                    ggTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCg                                                                             >  1:131548/1‑101 (MQ=255)
                    ggTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCg                                                                             <  2:131219/101‑1 (MQ=255)
                      tACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTg                                                                           <  1:131076/101‑1 (MQ=255)
                       aCGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTaaaa                                                                                >  2:131552/1‑95 (MQ=255)
                       aCGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGa                                                                          <  2:131297/101‑1 (MQ=255)
                         gATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACt                                                                        >  1:131554/1‑101 (MQ=255)
                          aTTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTg                                                                       <  1:131085/101‑1 (MQ=255)
                          aTTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTg                                                                       >  2:131558/1‑101 (MQ=255)
                            ttGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCg                                                                     <  2:131251/101‑1 (MQ=255)
                            ttGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCg                                                                     >  1:131561/1‑101 (MQ=255)
                             tgatgaCCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCgg                                                                    <  1:131222/101‑1 (MQ=255)
                              gatgaCCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTaaaa                                                                                >  2:131562/1‑88 (MQ=255)
                               atgaCCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGat                                                                  <  2:131290/101‑1 (MQ=255)
                                tgaCCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGata                                                                 >  1:131567/1‑101 (MQ=255)
                                    cGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATccc                                                             <  2:131182/101‑1 (MQ=255)
                                     gAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCg                                                            >  2:131570/1‑101 (MQ=255)
                                     gAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCg                                                            <  1:130776/101‑1 (MQ=255)
                                       acgacgAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAt                                                          <  1:131186/101‑1 (MQ=255)
                                       acgacgAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAt                                                          >  2:131572/1‑101 (MQ=255)
                                       acgacgAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAt                                                          <  1:131261/101‑1 (MQ=255)
                                       acgacgAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAt                                                          <  1:131300/101‑1 (MQ=255)
                                       acgacgAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAt                                                          <  2:131136/101‑1 (MQ=255)
                                       acgacgAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAt                                                          >  1:131574/1‑101 (MQ=255)
                                        cgacgaACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAtt                                                         <  1:131048/101‑1 (MQ=255)
                                        cgacgaACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAtt                                                         >  1:131577/1‑101 (MQ=255)
                                         gacgaACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATtg                                                        >  2:131580/1‑101 (MQ=255)
                                         gacgaACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATtg                                                        <  1:131112/101‑1 (MQ=255)
                                          acgaACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATtgt                                                       >  1:131582/1‑101 (MQ=255)
                                           cgaACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATtgtg                                                      <  2:131475/101‑1 (MQ=255)
                                            gaACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTgg                                                     >  1:131585/1‑101 (MQ=255)
                                               cAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTggggg                                                  >  1:131590/1‑101 (MQ=255)
                                                aaTCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGa                                                 >  2:131591/1‑101 (MQ=255)
                                                 aTCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAtt                                                         <  1:131484/92‑1 (MQ=255)
                                                   cTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGc                                              >  1:131593/1‑101 (MQ=255)
                                                     ggCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTc                                            <  1:131389/101‑1 (MQ=255)
                                                     ggCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTc                                            <  1:131203/101‑1 (MQ=255)
                                                      gCAAAGTACTGCCCAAATGACACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCg                                           >  1:131596/1‑101 (MQ=255)
                                                       cAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGc                                          >  2:131599/1‑101 (MQ=255)
                                                         aaGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCa                                        <  2:131351/101‑1 (MQ=255)
                                                          aGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAg                                       >  1:131603/1‑101 (MQ=255)
                                                           gTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAgg                                      <  2:131373/101‑1 (MQ=255)
                                                            tACTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGc                                          >  1:131607/1‑96 (MQ=255)
                                                             aCTGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTc                                    >  2:131609/1‑101 (MQ=255)
                                                               tGCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCt                                  >  2:131611/1‑101 (MQ=255)
                                                                gCCCAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTc                                 >  1:131614/1‑101 (MQ=255)
                                                                 cccAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCg                                >  2:131617/1‑101 (MQ=255)
                                                                  ccAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGAt                                                          <  2:131107/74‑1 (MQ=255)
                                                                  ccAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGa                               <  2:131242/101‑1 (MQ=255)
                                                                   cAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGaa                              >  2:131620/1‑101 (MQ=255)
                                                                   cAAATGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGaa                              <  1:131130/101‑1 (MQ=255)
                                                                    aaaTGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAAc                             <  1:131273/101‑1 (MQ=255)
                                                                    aaaTGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAAc                             <  2:131037/101‑1 (MQ=255)
                                                                     aaTGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACa                            <  2:131191/101‑1 (MQ=255)
                                                                      aTGCCACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACaa                           >  2:131624/1‑101 (MQ=255)
                                                                         ccACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATg                        >  2:131627/1‑101 (MQ=255)
                                                                         ccACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATg                        >  1:131628/1‑101 (MQ=255)
                                                                          cACTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATcc                                                              >  1:131630/1‑62 (MQ=255)
                                                                           aCTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCt                      <  2:131424/101‑1 (MQ=255)
                                                                           aCTGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCt                      <  1:131361/101‑1 (MQ=255)
                                                                             tGTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTg                    >  2:131634/1‑101 (MQ=255)
                                                                              gTTCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGa                   <  1:131267/101‑1 (MQ=255)
                                                                               ttCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGaa                  >  1:131637/1‑101 (MQ=255)
                                                                                tCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAAc                 <  1:131237/101‑1 (MQ=255)
                                                                                tCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAAc                 >  2:131638/1‑101 (MQ=255)
                                                                                tCTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATAGCCCGATTGTGGGGGATg                                               <  2:131063/71‑1 (MQ=255)
                                                                                 cTGCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAAct                >  2:131641/1‑101 (MQ=255)
                                                                                   gCATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAActct              <  2:131228/101‑1 (MQ=255)
                                                                                    cATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCtt             >  1:131645/1‑101 (MQ=255)
                                                                                    cATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCtt             <  1:131247/101‑1 (MQ=255)
                                                                                    cATATCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCtt             <  1:131172/101‑1 (MQ=255)
                                                                                     atatCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTg            <  2:131308/101‑1 (MQ=255)
                                                                                     atatCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTg            >  1:131648/1‑101 (MQ=255)
                                                                                     atatCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTg            <  1:131101/101‑1 (MQ=255)
                                                                                      tatCCATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGt           >  2:131650/1‑101 (MQ=255)
                                                                                        tCGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCg         <  2:131260/101‑1 (MQ=255)
                                                                                         cGATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGt                                     >  1:131654/1‑72 (MQ=255)
                                                                                          gATATTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCa       <  1:131381/101‑1 (MQ=255)
                                                                                            tatTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAg     >  2:131658/1‑101 (MQ=255)
                                                                                            tatTGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAg     >  1:131659/1‑101 (MQ=255)
                                                                                              ttGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGaa   <  2:131404/101‑1 (MQ=255)
                                                                                              ttGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGaa   >  2:131663/1‑101 (MQ=255)
                                                                                               tGAACCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAAt  >  1:131665/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCCGATTGTGGGGGATGCTCGCCAGGTCCTCGAACAAATGCTTGAACTCTTGTCGCAAGAAT  >  NC_000913/86449‑86644

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