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.C88.3% 990.3 / 87.1 352coding (35/1818 nt)ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base . (39/2);  new base C (81/229);  total (121/231)
Fisher's exact test for biased strand distribution p-value = 2.56e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 369 total aligned reads displayed.

GAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGGCCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCAC  >  NC_000913/3951074‑3951248
                                                                                                |                                                                                   
gagaCAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCCCCACCCCTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTgg                                                                                      <  1:4851723/93‑1 (MQ=255)
 agacagacACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTggg                                                                                    >  2:7051257/1‑93 (MQ=255)
   acagacaCTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTggcgc                                                                                   <  1:412369/93‑4 (MQ=255)
             ggAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGaa                                                                         <  1:1569347/93‑1 (MQ=255)
                gTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGa                                                                      <  1:927478/93‑1 (MQ=255)
                gTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGa                                                                      <  2:4825269/93‑1 (MQ=255)
                gTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGa                                                                      <  2:5406808/93‑1 (MQ=255)
                 taaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGAc                                                                     <  1:103815/93‑1 (MQ=255)
                  aaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGAcc                                                                    <  1:741572/93‑1 (MQ=255)
                  aaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGAcc                                                                    <  1:353788/93‑1 (MQ=255)
                    ataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGa                                                                  <  1:845191/93‑1 (MQ=255)
                    ataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGa                                                                  <  1:3196846/93‑1 (MQ=255)
                    ataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGa                                                                  <  1:6065333/93‑1 (MQ=255)
                     taaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGAc                                                                 <  2:5333102/93‑1 (MQ=255)
                       aaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGc                                                               <  2:957574/93‑1 (MQ=255)
                        aGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGcc                                                              <  2:2663605/93‑1 (MQ=255)
                        aGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGCGTGCTCCTGCGC‑‑TGTGGCGCGCCACCGGGATGACCGACGcc                                                              >  1:5906495/1‑93 (MQ=255)
                          tATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCgc                                                            >  2:3662612/1‑92 (MQ=255)
                           aTGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGAt                                                           <  2:1412959/93‑1 (MQ=255)
                            tGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGAtt                                                          <  2:5986847/93‑1 (MQ=255)
                              ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTc                                                        <  1:6983791/93‑1 (MQ=255)
                              ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTc                                                        <  1:986790/93‑1 (MQ=255)
                                tAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCgg                                                      <  1:2046976/93‑1 (MQ=255)
                                 aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGGCGTGCGCCGGGGGGGGGCCCCGGGAAGACCGACGCAGACTTCGGt                                                     >  1:7516668/1‑93 (MQ=255)
                                  aGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTa                                                    <  1:3000627/93‑1 (MQ=255)
                                  aGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCCGAATGACCGACGCCGCTTTCggga                                                    >  2:1658526/1‑91 (MQ=255)
                                    tACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAg                                                  <  1:5229961/93‑1 (MQ=255)
                                     aCCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGc                                                 <  1:6343984/93‑1 (MQ=255)
                                     aCCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTGGGGCGCCACCGGGATGACCGACGCCGATTTTGGTAAGc                                                 >  2:2992113/1‑93 (MQ=255)
                                      ccGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGcc                                                <  2:3222798/93‑1 (MQ=255)
                                      ccGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGGGGGC‑‑TGTGGGGCGCCCCCGGGATGACCGACCCCGACTTCGGTAAGcc                                                >  2:126234/1‑93 (MQ=255)
                                       cGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCg                                               <  1:7289263/93‑1 (MQ=255)
                                       cGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGGGC‑‑CGTGGCGCGCCACCCGAAAGACCGACGCCCATTTCGGTAAGCCg                                               >  2:6868591/1‑93 (MQ=255)
                                        gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGa                                              <  1:368293/93‑1 (MQ=255)
                                        gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGa                                              <  1:2628252/93‑1 (MQ=255)
                                         ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAt                                             <  2:4374603/93‑1 (MQ=255)
                                         ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAt                                             <  1:5304233/93‑1 (MQ=255)
                                         ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTGGCGCGCCCCCGGAATGACCGACGCCGATTTTGGTAAGACGAt                                             >  2:6375495/1‑93 (MQ=255)
                                         ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTGGCGCGCCACCGGGATGACCGACGCCGCTTTCGGTAAGCCGAt                                             >  2:2059183/1‑93 (MQ=255)
                                          tCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAtt                                            <  1:1120083/93‑1 (MQ=255)
                                            cGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTAt                                          <  2:2032299/93‑1 (MQ=255)
                                             gCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                                         <  2:4886070/93‑1 (MQ=255)
                                             gCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                                         <  2:5446489/93‑1 (MQ=255)
                                             gCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTGGGGCGCCCCCCGGATGAACGACCCCGGTTTTGGGAAGCCGATTTTc                                         >  2:5154059/1‑93 (MQ=255)
                                              ccaccaccaCTCATGGTCGTAATATGGCGGGGGGTCGTGCGC‑‑TGTGGCGCGCCACCCGGATGACCGACCCCCCCTTCCGTTAACCCATTATcc                                        >  1:5481670/1‑92 (MQ=255)
                                               caccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgc                                       <  1:3866044/93‑1 (MQ=255)
                                               caccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgc                                       <  2:4228915/93‑1 (MQ=255)
                                               caccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgc                                       <  1:4587403/93‑1 (MQ=255)
                                                accaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgcg                                      <  1:2753093/93‑1 (MQ=255)
                                                 ccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCgg                                     <  2:2169291/93‑1 (MQ=255)
                                                  caccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGt                                    <  2:2838864/93‑1 (MQ=255)
                                                  caccacTCATGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGt                                    >  2:853079/1‑93 (MQ=255)
                                                   accacTCATGGTCGTAATATGGCGGGGGGTCGTGCGC‑‑TGTGGGGCGCCACCCGGCTGAACGGCGCCGAATTCGCTAAGCCCATTATtgcccgt                                   >  2:3464763/1‑89 (MQ=255)
                                                   accacTCATGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTGGGGCGCCCCCGGGATGACCGACGCCGATTTTGGTAAGGCGATTATCGCGGtt                                   >  2:1712497/1‑93 (MQ=255)
                                                    ccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TCTGGCGCGCCACCGGCAAGACCGACCCCGCATTTGGGAAGCCGGTTATCCCGGAtg                                  >  2:2270777/1‑93 (MQ=255)
                                                    ccacTCATGGTCGTAATATGGCGGGGGCTCGCGGGC‑‑TGTGGCGCGCCACCGGAAAGACCCAAGCCGATTTTCGTAAGCCGAATATCGCGGTtg                                  >  2:4027117/1‑93 (MQ=255)
                                                     cacTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgt                                 <  2:7275722/93‑1 (MQ=255)
                                                       cTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGa                               <  1:18116/93‑1 (MQ=255)
                                                       cTCATGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTGGCGCGCCACCCGGATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGGTGTGa                               >  2:6668332/1‑93 (MQ=255)
                                                          aTGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACt                            <  1:1931831/93‑1 (MQ=255)
                                                           tGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc                           <  2:6217645/93‑1 (MQ=255)
                                                           tGGTCGTAATATGGCGGGGGCTCGTGCGC‑‑TGTGGCGCGCCCCCGGAATGACCGACCCCGACTTTGGTAAGCCGATTATTGCGGTTGTGGACTc                           >  1:5815698/1‑93 (MQ=255)
                                                            ggTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCg                          <  2:2124657/93‑1 (MQ=255)
                                                             gTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGt                         <  2:5820597/93‑1 (MQ=255)
                                                             gTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGt                         <  1:6453782/93‑1 (MQ=255)
                                                             gTCGTAATATGGCGGGTGCTCGTGCCC‑‑CGTGGCGCGCCACCCGGATGACCGACCCCCATTTCTGGAAGCCGAATATCGCGGGTGTGAACTCGt                         >  2:1835021/1‑93 (MQ=255)
                                                              tCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGtt                        <  2:4020692/93‑1 (MQ=255)
                                                              tCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGtt                        <  2:7165877/93‑1 (MQ=255)
                                                               cGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc                       <  1:2415418/93‑1 (MQ=255)
                                                                 tAATATGGCGGGGGCTCGTGCGC‑‑TGTGGCGCGCCACCGGGATGAACGACGCCGATTTCGGTAAGCCGATTATCGCGGGTGTGAACTCGTTCAc                     >  2:7379058/1‑93 (MQ=255)
                                                                  aaTATGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAcc                    <  1:229251/93‑1 (MQ=255)
                                                                   atatGGCGGGGGCTCGTGCGC‑‑TGTGGGGCGCCCCCGGAATGACCGACGCCGATTTCTGGAAGCCGATTATTGCGGGTGTGAACTCGTTTAccc                   >  2:3819537/1‑93 (MQ=255)
                                                                     atGGCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCaa                 <  2:1102658/93‑1 (MQ=255)
                                                                       ggCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAAtt               >  1:3451728/1‑93 (MQ=255)
                                                                        gCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAAttt              <  2:4552445/93‑1 (MQ=255)
                                                                        gCGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAAttt              <  2:4965476/93‑1 (MQ=255)
                                                                         cGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTg             <  2:2505553/93‑1 (MQ=255)
                                                                         cGGGTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTg             <  2:5625899/93‑1 (MQ=255)
                                                                           ggTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTa           >  1:1304150/1‑93 (MQ=255)
                                                                            gTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc          <  2:3586425/93‑1 (MQ=255)
                                                                            gTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc          <  2:6603125/93‑1 (MQ=255)
                                                                            gTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc          <  2:6512653/93‑1 (MQ=255)
                                                                            gTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc          <  2:2733269/93‑1 (MQ=255)
                                                                            gTGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc          <  1:213888/93‑1 (MQ=255)
                                                                             tGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc         <  1:3128340/93‑1 (MQ=255)
                                                                             tGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc         <  1:1467639/93‑1 (MQ=255)
                                                                             tGCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc         <  2:6279802/93‑1 (MQ=255)
                                                                              gCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg        <  2:3161799/93‑1 (MQ=255)
                                                                              gCTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg        >  1:444801/1‑93 (MQ=255)
                                                                               cTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCgg       <  1:4013517/93‑1 (MQ=255)
                                                                               cTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCgg       >  2:1899041/1‑93 (MQ=255)
                                                                               cTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCgg       <  2:1331193/93‑1 (MQ=255)
                                                                               cTCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATATCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCgg       >  2:2932866/1‑93 (MQ=255)
                                                                                tCGTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCggg      <  2:3290640/93‑1 (MQ=255)
                                                                                  gTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTc    <  2:1206985/93‑1 (MQ=255)
                                                                                  gTGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTc    >  1:6735586/1‑93 (MQ=255)
                                                                                   tGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCa   >  2:4733837/1‑93 (MQ=255)
                                                                                   tGCGC‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCa   >  2:1065487/1‑93 (MQ=255)
                                                                                    gcgc‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCAc  >  2:6808347/1‑93 (MQ=255)
                                                                                    gcgc‑‑TGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCAc  >  2:1587199/1‑93 (MQ=255)
                                                                                                |                                                                                   
GAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGC‑‑TGTGGGCCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCAC  >  NC_000913/3951074‑3951248

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