Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,166 +C coding (35/1818 nt) ilvD → dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,1661.C86.8% 211.2 / 32.3 92coding (35/1818 nt)ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base . (11/1);  new base C (25/54);  total (37/55)
Fisher's exact test for biased strand distribution p-value = 1.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑G‑‑‑GCC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC  >  NC_000913/3951093‑3951231
                                                                             |                                                                      
aataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGTg‑‑g‑‑‑gg                                                                    >  1:246450/1‑75 (MQ=255)
 ataacGTAGGCCTAAGTCCCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTg‑‑gcgc                                                                      <  2:2837871/76‑4 (MQ=255)
 ataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTg‑‑g‑‑‑ggc                                                                   >  1:3338266/1‑74 (MQ=255)
 ataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCCCGTGCGCCGTg‑‑gcgc                                                                      >  2:4547052/1‑73 (MQ=255)
 ataaaGGAGGCCTAGGGGGCGGTCCGCCACCACCACTCGTGGTGGGAATATGGCGGGTGCTCGTGCGCTGTg‑‑gcgc                                                                      <  2:4243/76‑4 (MQ=255)
        aTGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCCTGCCGCGTG‑‑G‑‑‑Gcc‑‑ccccccc                                                         >  2:907197/1‑73 (MQ=255)
          gCCTAAGTACCGTGCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGa                                                         <  2:2577079/76‑1 (MQ=255)
           ccTAAGTACCGTCCCGCCCCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGaa                                                        <  2:2275969/76‑1 (MQ=255)
           ccTAAGTACCGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGaa                                                        <  1:3348697/76‑1 (MQ=255)
              aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGa                                                     <  2:1370389/76‑1 (MQ=255)
               aGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGAc                                                    <  1:4064661/76‑1 (MQ=255)
                gtcccGTCCCCCCCCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGAcc                                                   <  2:1777903/73‑1 (MQ=255)
                gccccGTTCCCCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGAcc                                                   <  2:5831407/73‑1 (MQ=255)
                 tcccGTCCCGCCCCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCg                                                  <  2:1203785/74‑1 (MQ=255)
                 tACCGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGG‑‑CACCGGAATGACCg                                                  <  2:4192007/76‑1 (MQ=255)
                  ccggTTCCGCCCCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGa                                                 <  2:5526582/73‑1 (MQ=255)
                  cccGTTCCGCCCCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGa                                                 <  2:5299247/75‑1 (MQ=255)
                  cccGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGa                                                 <  1:1166889/75‑1 (MQ=255)
                   ccGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGG‑‑G‑‑‑GCGC‑‑CACCGGGATGAACGAc                                                >  1:3536070/1‑76 (MQ=255)
                    cGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGG‑‑G‑‑‑GCGC‑‑CCCCCGAATGACCgccc                                               >  1:1835039/1‑73 (MQ=255)
                     gtcccGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGc                                              <  2:411285/73‑1 (MQ=255)
                     gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGTGG‑‑G‑‑‑GCGC‑‑CCCCCGGATGACCGACGc                                              >  1:2925591/1‑76 (MQ=255)
                     gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGG‑‑G‑‑‑GCGC‑‑CCCCGGGATGACCGACGc                                              >  1:2340145/1‑76 (MQ=255)
                     gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCCGTG‑‑GC‑‑GCGC‑‑CCCCGGCCTGACCGACGc                                              >  2:2874703/1‑76 (MQ=255)
                      tcccGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGcc                                             <  1:4478849/74‑1 (MQ=255)
                       tCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCg                                            <  2:4084214/76‑1 (MQ=255)
                        ccGCCCCCACCTCTCATGGTCGTAATAGGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGa                                           <  1:3466486/76‑1 (MQ=255)
                          gCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGAtt                                         <  2:1698447/76‑1 (MQ=255)
                          gCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGG‑‑G‑‑‑GCGC‑‑CCCCCGGATGACCGCCGCCGCtt                                         >  2:4364198/1‑76 (MQ=255)
                           ccaccaccaCTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGAttt                                        <  1:1274035/76‑1 (MQ=255)
                           ccaccaccaCTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGAttt                                        <  1:4310875/76‑1 (MQ=255)
                            ccccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTc                                       <  1:5510897/74‑1 (MQ=255)
                            caccaccacTCATGGTCGTAATATGGCGGGTGGTCGGGCGCTGTGG‑‑G‑‑‑GCGC‑‑CCCCGGGAAGAACGACGCCGAnttt                                       >  1:4348252/1‑75 (MQ=255)
                            caccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CCCCGGGATGACCCCCCCCCATTTc                                       >  1:2065393/1‑76 (MQ=255)
                             accaccacTCATGGTCGTAATGTGGCGGGGGCTCGTGCGCTGTGG‑‑G‑‑‑GCGC‑‑CCCCGGGATGACCGACGCCGCTTTCg                                      >  1:4091545/1‑76 (MQ=255)
                             accaccacTCATGGTCGTAATATGGCGGGTGGTCGTGGGCTGTGGCG‑‑‑GGC‑‑CCCCCGAATGACCGACGGCGATTTCg                                      >  1:1982334/1‑76 (MQ=255)
                             accaccacTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CCCCGGGATGACCGACGCCGATTTCg                                      >  1:3672105/1‑76 (MQ=255)
                              ccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCgg                                     <  1:2422018/76‑1 (MQ=255)
                               caccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGt                                    <  2:3119328/76‑1 (MQ=255)
                                cccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTa                                   <  2:5572862/75‑1 (MQ=255)
                                accacTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTG‑‑G‑‑‑GGCGGCCCCGGAATGACCCACCCCCGTTTCGGTa                                   >  1:2649028/1‑76 (MQ=255)
                                  cgcTCATTGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CCCCGGAATGACCGACGCCGATTTCGGTAAg                                 <  1:5101290/74‑1 (MQ=255)
                                  cacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAg                                 <  1:2429851/76‑1 (MQ=255)
                                   acTCATGGTCGTAATATGGCGTGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGc                                <  2:2909650/76‑1 (MQ=255)
                                   acTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGc                                <  2:1025977/76‑1 (MQ=255)
                                   acTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTG‑‑GGGCGCC‑‑CCCCGGAAGACCCACGCCGCTTTTGGGAAgg                                >  1:3615311/1‑75 (MQ=255)
                                     tCATGGTCGTAATATGGCGGGTGGTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CCCCGGGATGACCGACGCCGACTTCGGTAAGCCg                              >  2:4010536/1‑76 (MQ=255)
                                     tCATGGTCGTAATATGGCGGGTGCTCGGGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCg                              <  2:3671179/76‑1 (MQ=255)
                                      cATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGa                             <  1:677167/76‑1 (MQ=255)
                                      cATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGa                             <  1:455680/76‑1 (MQ=255)
                                       aTGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGAt                            <  1:2586050/76‑1 (MQ=255)
                                       aTGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGAt                            <  1:5082729/76‑1 (MQ=255)
                                        tGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGAtt                           <  2:782532/76‑1 (MQ=255)
                                          gTCGTAATATGGCGGGGGCTCGTGCGCTGTGG‑‑G‑‑‑GCGC‑‑CCCCGGGATGACCGACGCCGGTTTTGGTGAGGCGGTTAt                         >  1:5229260/1‑76 (MQ=255)
                                           tCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                        <  1:4364199/76‑1 (MQ=255)
                                           tCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                        <  1:5222443/76‑1 (MQ=255)
                                           tCGTAATATGGCGGGGGGTCGTGCGCTGGG‑‑GC‑‑GCGC‑‑CCCCGGGATGACCGACCCCCGTTTTGGTAAAGCGATTAtt                        >  1:3930809/1‑75 (MQ=255)
                                             gTAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTTTGTAAGCCGATTATcgc                      >  1:949511/1‑76 (MQ=255)
                                             gTAATAGGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgc                      <  2:5244440/76‑1 (MQ=255)
                                             ctAATATGGCGGTTGCTCGTGGGCTGTGTC‑‑‑GCGTCACCGGAATGCCCGACGCCGATTTCGGTAAGGCGATTATcgc                      <  2:2219313/75‑1 (MQ=255)
                                              tAATATGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgcg                     <  1:1055146/76‑1 (MQ=255)
                                               aaTATGGCGGGGGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCgg                    <  2:4973827/76‑1 (MQ=255)
                                                atatGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGt                   >  1:4960406/1‑76 (MQ=255)
                                                atatGGCGGGTGCTCGTGAGCTGTG‑‑GC‑‑GCGC‑‑CACAGGAATGACCGACGCCGATTTCGCGAAGCCGATTATAGCGGt                   >  1:30131/1‑76 (MQ=255)
                                                 tatGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑TCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGtt                  <  1:380530/76‑1 (MQ=255)
                                                  atGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtg                 <  2:3264656/76‑1 (MQ=255)
                                                   tGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgt                <  2:1928397/76‑1 (MQ=255)
                                                   tGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgt                <  2:5229259/76‑1 (MQ=255)
                                                   tGGCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGGtgg                >  1:629795/1‑75 (MQ=255)
                                                    ggCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgtg               >  2:2644272/1‑76 (MQ=255)
                                                    ggCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgtg               >  2:1320452/1‑76 (MQ=255)
                                                    ggCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTAACGCGGTtgtg               <  1:3701123/76‑1 (MQ=255)
                                                     gCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGa              >  1:2378324/1‑76 (MQ=255)
                                                     gCGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGa              <  1:1162960/76‑1 (MQ=255)
                                                      ggggTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGaa             <  2:2600434/75‑1 (MQ=255)
                                                      cGGGTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGaa             >  2:5206173/1‑76 (MQ=255)
                                                       gggTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAAc            <  2:3537152/76‑1 (MQ=255)
                                                       gggTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAAc            <  1:1005863/76‑1 (MQ=255)
                                                       gggTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAAc            >  1:3198377/1‑76 (MQ=255)
                                                       gggTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGGTGTGAAc            >  2:1657870/1‑76 (MQ=255)
                                                       gggTGCTCGTGCGCTGTG‑‑GC‑‑GCAC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGacc            >  2:505895/1‑74 (MQ=255)
                                                       gggTGCTCGTGCGCTGTG‑‑GC‑‑GCAC‑‑CACCGGAATGACCGACGCAGATTTATGAACGCCGATTATCGCGGTtgtgccc            >  2:511618/1‑73 (MQ=255)
                                                       gggTGATCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGACGATTTCGGTAAGCCGATTATCGCGGTTGTGAAc            >  1:5281988/1‑76 (MQ=255)
                                                         gTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          >  1:5051100/1‑76 (MQ=255)
                                                         gTGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          <  2:3543493/76‑1 (MQ=255)
                                                          tGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCg         >  2:4816072/1‑76 (MQ=255)
                                                          tGCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCg         <  1:2193243/76‑1 (MQ=255)
                                                           gCTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGt        <  1:4483767/76‑1 (MQ=255)
                                                            cTCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGtt       <  1:796741/76‑1 (MQ=255)
                                                             tCGTTCGCTGTG‑‑GC‑‑GCTC‑‑CTCCGGAATTGCCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCTTTc      <  1:2975304/76‑1 (MQ=255)
                                                             tCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGGACTCGTTc      >  2:361189/1‑76 (MQ=255)
                                                             tCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc      <  1:2560324/76‑1 (MQ=255)
                                                             tCGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc      >  2:1820191/1‑76 (MQ=255)
                                                              cGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCa     <  1:40876/76‑1 (MQ=255)
                                                              cGTGCGCTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGGAACCCGATTATCGCGGTTGTGAACCCCTTCa     >  2:4446849/1‑76 (MQ=255)
                                                               gTGCGCTGTG‑‑GC‑‑GCGC‑‑CCCCGGAATGACCGACGCCGATTTCGGTAAACCGATTATCGCGGGTGGGAAAGCGTTCAc    >  2:4395995/1‑76 (MQ=255)
                                                                 gcgcTGTG‑‑GC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAccc  >  2:2685544/1‑76 (MQ=255)
                                                                             |                                                                      
AATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTG‑‑G‑‑‑GCC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC  >  NC_000913/3951093‑3951231

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