Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,170 Δ1 bp coding (39/1818 nt) ilvD → dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,1680C.99.2% 1077.7 / ‑0.5 242coding (37/1818 nt)ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (27/213);  total (28/214)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

Only 100 of 308 total aligned reads displayed.

TGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCG‑‑CTGTG‑‑G‑‑‑GCCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCC  >  NC_000913/3951085‑3951252
                                                                                           |                                                                                     
gatactgggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGcc                                                                                     <  2:1196687‑M2/3‑1 (MQ=255)
gatactgggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGcc                                                                                     <  2:3146262‑M2/3‑1 (MQ=255)
  cactgggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcgccac                                                                                   <  1:7001471‑M2/5‑3 (MQ=255)
      gggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGa                                                                              <  2:7178046‑M2/9‑1 (MQ=255)
       ggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGaa                                                                             <  1:1769531‑M2/10‑1 (MQ=255)
        gagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAAt                                                                            <  1:5117007‑M2/11‑1 (MQ=255)
        gagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAAt                                                                            <  2:9153869‑M2/11‑1 (MQ=255)
        gagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAAt                                                                            <  2:8838398‑M2/11‑1 (MQ=255)
          gtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcggggggtcgtgcgctgtgggggGCCCCCcggatga                                                                           >  2:5763518‑M2/81‑90 (MQ=255)
           taaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGAc                                                                         <  1:6727178‑M2/14‑1 (MQ=255)
           taaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGAc                                                                         <  2:8447498‑M2/14‑1 (MQ=255)
           taaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcggggggtcgtgcgctgtggcgcGCCCCccggatgac                                                                          >  1:4837596‑M2/80‑89 (MQ=255)
           taaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcggggggtcctgcgccgtggggcGCCCCccggatgaa                                                                          >  2:6755443‑M2/80‑89 (MQ=255)
            aaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggggctcgcgcgctgtggcgcGCCACCCGGATGAcc                                                                        >  2:5219051‑M2/79‑93 (MQ=255)
             aataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCg                                                                       <  2:5081308‑M2/16‑1 (MQ=255)
             aataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCg                                                                       <  2:3012093‑M2/16‑1 (MQ=255)
              ataaagtttgccttaggtcccttcccccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGa                                                                      <  1:7244760‑M2/17‑1 (MQ=255)
               taaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGAc                                                                     <  1:2645545‑M2/18‑1 (MQ=255)
                  agtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGcc                                                                  <  2:4201615‑M2/21‑1 (MQ=255)
                      aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTCTGGG‑‑G‑‑‑GGGCCCCGGCATGAACGACGCCGCtttttgt                                                         >  2:9512961/1‑89 (MQ=255)
                       aGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCG‑‑CCGCG‑‑GC‑‑GCGCCCCCGGGATGAACCACGCCGGTTTCGGTa                                                        >  2:3891894/1‑93 (MQ=255)
                        cctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTc                                                            <  1:7406808‑M2/27‑1 (MQ=255)
                        cctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggggctcgtgcgctgtggggcGGCCCCCGACAGACCGACGCCGacttt                                                            >  1:4431351‑M2/67‑92 (MQ=255)
                        cctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggggctcctgcgctgtggggcGCCACCGGAATGACCGACCCCCatttt                                                            >  2:24217‑M2/67‑92 (MQ=255)
                          taagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggggcGGCACCGGAAGGACCGACCCCGATTTCgg                                                          >  2:3494940‑M2/65‑93 (MQ=255)
                          caagtacccttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCgg                                                          <  2:3549762‑M2/29‑1 (MQ=255)
                           ccGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCCGGCG‑‑CTGTGG‑‑G‑‑‑GCGCCCCCGGGATGAACTACCCCCATTTCGGTAAAcc                                                    >  2:8245668/1‑93 (MQ=255)
                           aagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGt                                                         <  1:9520080‑M2/30‑1 (MQ=255)
                           aagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGt                                                         <  2:1657489‑M2/30‑1 (MQ=255)
                            agtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTa                                                        <  2:741345‑M2/31‑1 (MQ=255)
                            agtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTa                                                        <  2:8375335‑M2/31‑1 (MQ=255)
                            agtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTa                                                        <  1:8759809‑M2/31‑1 (MQ=255)
                            agtaccgttccgccaccaccactcatggtcgtaatatggcgggggctcgtgcgccgtggggcGCCCCCCGGAAGAACGACGCCGCTTttggga                                                        >  2:7999245‑M2/63‑91 (MQ=255)
                             gtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTaa                                                       <  1:6917251‑M2/32‑1 (MQ=255)
                             gTTCCGCCACCGCCACTCATGGTCGTAATAGGGCGGGCGCTCGTCCG‑‑CTGTGCC‑‑‑GCGCCACCGAAATGGCCGACGCCGATTTCGGTAAGCCGa                                                  <  2:6042073/93‑1 (MQ=255)
                             gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCCTGCG‑‑CCGTG‑‑GGGCGCCCCCCGGAAGAACGACGCCGCCTTCGGGAAGGCgg                                                  >  2:4941164/1‑92 (MQ=255)
                              ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCG‑‑CCGTGGCG‑‑‑GGCCACCGGGATGACCGACCCCGCGTTCGGTAAAcccgc                                                 >  1:2051575/1‑90 (MQ=255)
                              taccgttccgccaccaccactcatggtcgtaatatggcgggggctcgggcgctgtggcgcGCCACCGGAATGACCGACGCCGAGTTCGGTCAg                                                      >  2:55644‑M2/61‑93 (MQ=255)
                               accgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGc                                                     <  2:4495654‑M2/34‑1 (MQ=255)
                                  cttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGa                                                  <  2:1529503‑M2/37‑1 (MQ=255)
                                   ccaccaccaCTCATGGTCGTAATATGGCGGGGGGTCGTGCG‑‑CTGTG‑‑GCGCGCCCCCCGGAAGAACGACGCCGACTTCCGTAAGCCGATTAttc                                            >  2:2285802/1‑91 (MQ=255)
                                    tccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAtt                                                <  1:1979522‑M2/39‑1 (MQ=255)
                                     ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                               <  2:436758‑M2/40‑1 (MQ=255)
                                     ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                               <  2:2068417‑M2/40‑1 (MQ=255)
                                     ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                               <  1:1373520‑M2/40‑1 (MQ=255)
                                     ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                               <  1:6515661‑M2/40‑1 (MQ=255)
                                     ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                               <  1:6870954‑M2/40‑1 (MQ=255)
                                     ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                               <  1:8372525‑M2/40‑1 (MQ=255)
                                      cgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTAt                                              <  1:754664‑M2/41‑1 (MQ=255)
                                       gccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                                             <  1:4012462‑M2/42‑1 (MQ=255)
                                       gccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                                             <  1:1000378‑M2/42‑1 (MQ=255)
                                       cccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                                             <  1:7842296‑M2/42‑1 (MQ=255)
                                       cccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                                             <  1:2254575‑M2/42‑1 (MQ=255)
                                        ccaccagcactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGCCGCCGATTTCGGTAAGCCGATTATcg                                            <  2:6946608‑M2/43‑1 (MQ=255)
                                         caccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgc                                           <  2:5343164‑M2/44‑1 (MQ=255)
                                         caccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgc                                           <  2:4302623‑M2/44‑1 (MQ=255)
                                           ccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCgg                                         <  1:8900437‑M2/46‑1 (MQ=255)
                                            caccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGt                                        <  1:2971577‑M2/47‑1 (MQ=255)
                                             cccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGtt                                       <  1:9914012‑M2/48‑1 (MQ=255)
                                              ccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtg                                      <  1:8409719‑M2/49‑1 (MQ=255)
                                               cactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgt                                     <  2:2364421‑M2/50‑1 (MQ=255)
                                                tGGTCGTAATATGGCGGGGGCTCGTGCG‑‑CTGTGG‑‑G‑‑‑GCGCCCCCCGGATGAACCACCCCGAATTTGGTAAGCCGATTATTGCGGGTGTGGACTc                               >  2:3712728/1‑93 (MQ=255)
                                                actcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgtg                                    <  1:8132132‑M2/51‑1 (MQ=255)
                                                 ctcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGa                                   <  1:1442584‑M2/52‑1 (MQ=255)
                                                  tcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGaa                                  <  2:5676805‑M2/53‑1 (MQ=255)
                                                  tcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGaa                                  <  1:9082085‑M2/53‑1 (MQ=255)
                                                  tcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGaa                                  <  1:6120856‑M2/53‑1 (MQ=255)
                                                  tcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGaa                                  <  1:6202106‑M2/53‑1 (MQ=255)
                                                    atggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACt                                <  1:7285268‑M2/55‑1 (MQ=255)
                                                    atggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACt                                <  1:3681351‑M2/55‑1 (MQ=255)
                                                    atggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAATATCGAGGTTGAGAaca                                <  2:8986892‑M2/55‑2 (MQ=255)
                                                        tcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGtt                            <  1:3603922‑M2/59‑1 (MQ=255)
                                                         tatGGCGGGTGCTCGTGCG‑‑CTGGG‑‑GC‑‑GCGCCCCCGGAATGAACGACGCCGCTTTCCGTAAGCCGGTTAACCAGGTTGTGAACTCGTTCACCCa                      >  2:7726714/1‑93 (MQ=255)
                                                          gtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCa                          <  1:354722‑M2/61‑1 (MQ=255)
                                                          gtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCa                          <  2:9661990‑M2/61‑1 (MQ=255)
                                                            aatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAcc                        <  1:8562970‑M2/63‑1 (MQ=255)
                                                             atatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAccc                       <  2:9382451‑M2/64‑1 (MQ=255)
                                                              tatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCa                      <  1:5585093‑M2/65‑1 (MQ=255)
                                                              tatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCa                      >  1:80842‑M2/29‑93 (MQ=255)
                                                              tatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGGTGTGAAATCGTTCACCCa                      >  2:9879865‑M2/29‑93 (MQ=255)
                                                                 ggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCCAtt                   >  2:4819434‑M2/26‑93 (MQ=255)
                                                                 ggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAAtt                   <  1:2713315‑M2/68‑1 (MQ=255)
                                                                 ggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAAtt                   <  1:3164634‑M2/68‑1 (MQ=255)
                                                                   cgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTg                 <  1:6755444‑M2/70‑1 (MQ=255)
                                                                    gggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGt                <  1:3911372‑M2/71‑1 (MQ=255)
                                                                     ggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTa               >  1:5402462‑M2/22‑93 (MQ=255)
                                                                     ggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTa               <  2:1756244‑M2/72‑1 (MQ=255)
                                                                      gtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc              <  1:3302602‑M2/73‑1 (MQ=255)
                                                                       tgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc             <  1:4137568‑M2/74‑1 (MQ=255)
                                                                       tgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc             <  1:1622792‑M2/74‑1 (MQ=255)
                                                                       tgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc             <  2:8934620‑M2/74‑1 (MQ=255)
                                                                       tgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc             >  2:1966533‑M2/20‑93 (MQ=255)
                                                                       tgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc             <  1:9754887‑M2/74‑1 (MQ=255)
                                                                        gctcgtgcgctgtggcgcGCCACCAGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg            <  2:8051060‑M2/75‑1 (MQ=255)
                                                                          tcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCggg          >  2:5704715‑M2/17‑93 (MQ=255)
                                                                          tcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCggg          >  2:9242456‑M2/17‑93 (MQ=255)
                                                                             tgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCa       >  2:2622365‑M2/14‑93 (MQ=255)
                                                                             tgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCa       >  2:1284284‑M2/14‑93 (MQ=255)
                                                                             tgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCa       <  1:5698594‑M2/80‑1 (MQ=255)
                                                                                 ctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTc   <  1:9624313‑M2/84‑1 (MQ=255)
                                                                                           |                                                                                     
TGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCG‑‑CTGTG‑‑G‑‑‑GCCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCC  >  NC_000913/3951085‑3951252

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