Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,754,310 C→A R382M (AGG→ATG)  yiaY ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,754,3100CA99.7% 1312.9 / NA 339R382M (AGG→ATG) yiaYputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (172/166);  minor base G (0/1);  total (172/167)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

GACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTGCATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAA  >  NC_000913/3754220‑3754398
                                                                                          |                                                                                        
gACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATc                                                                                        >  2:2264226/1‑93 (MQ=255)
gACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATc                                                                                        <  1:96605/93‑1 (MQ=255)
 aCAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCg                                                                                       >  2:4926602/1‑93 (MQ=255)
 aCAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCg                                                                                       >  2:6871107/1‑93 (MQ=255)
  cAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGc                                                                                      <  1:622887/93‑1 (MQ=255)
     cTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGc                                                                                   <  1:5912152/93‑1 (MQ=255)
      tATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCg                                                                                  >  2:4942528‑M1/1‑92 (MQ=255)
       aTCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgc                                                                                 <  1:2245839‑M1/93‑3 (MQ=255)
        tCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcg                                                                                <  1:2771983‑M1/93‑4 (MQ=255)
        tCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcg                                                                                >  2:3603655‑M1/1‑90 (MQ=255)
         cACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcga                                                                               >  1:2480256‑M1/1‑89 (MQ=255)
         cACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcga                                                                               <  2:2748038‑M1/93‑5 (MQ=255)
          aCGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgat                                                                              <  2:666693‑M1/93‑6 (MQ=255)
          aCGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgat                                                                              >  1:1796193‑M1/1‑88 (MQ=255)
             aaTTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaa                                                                           <  1:4323074‑M1/93‑9 (MQ=255)
              aTTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaat                                                                          <  2:3887438‑M1/93‑10 (MQ=255)
              aTTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaat                                                                          >  2:1618618‑M1/1‑84 (MQ=255)
               tttACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                         <  1:1044999‑M1/93‑11 (MQ=255)
                 tACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgc                                                                       <  1:5368644‑M1/93‑13 (MQ=255)
                  aCGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgcc                                                                      <  1:2957143‑M1/93‑14 (MQ=255)
                  aCGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgcc                                                                      <  1:1953926‑M1/93‑14 (MQ=255)
                   cGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgcca                                                                     >  1:1337853‑M1/1‑79 (MQ=255)
                    gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccac                                                                    <  2:172324‑M1/93‑16 (MQ=255)
                    gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccac                                                                    <  1:6293992‑M1/93‑16 (MQ=255)
                    gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccac                                                                    >  1:7226572‑M1/1‑78 (MQ=255)
                      gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaa                                                                  <  2:4843160‑M1/93‑18 (MQ=255)
                      gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaa                                                                  >  2:466160‑M1/1‑76 (MQ=255)
                        aaGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatt                                                                >  2:6014116‑M1/1‑74 (MQ=255)
                        aaGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatt                                                                <  1:5562346‑M1/93‑20 (MQ=255)
                         aGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaattt                                                               >  1:406363‑M1/1‑73 (MQ=255)
                          gggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttc                                                              >  2:4524589‑M1/1‑72 (MQ=255)
                           ggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttct                                                             >  1:314185‑M1/1‑71 (MQ=255)
                            gTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttctt                                                            <  1:2722765‑M1/93‑24 (MQ=255)
                             ttACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttc                                                           >  1:5152962‑M1/1‑69 (MQ=255)
                                cAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtg                                                        >  2:1945824‑M1/1‑66 (MQ=255)
                                 agagATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtga                                                       >  2:312351‑M1/1‑65 (MQ=255)
                                 agagATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtga                                                       >  1:5700264‑M1/1‑65 (MQ=255)
                                 agagATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtga                                                       <  2:4728430‑M1/93‑29 (MQ=255)
                                  gagaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgag                                                      >  2:7362737‑M1/1‑64 (MQ=255)
                                   agaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagt                                                     <  1:5247717‑M1/93‑31 (MQ=255)
                                    gaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagtt                                                    >  2:3556974‑M1/1‑62 (MQ=255)
                                       aGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcc                                                 <  1:2108539‑M1/93‑35 (MQ=255)
                                       aGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcc                                                 >  1:740975‑M1/1‑59 (MQ=255)
                                         ttGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctg                                               >  2:3413399‑M1/1‑57 (MQ=255)
                                          tGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctgg                                              <  2:855991‑M1/93‑38 (MQ=255)
                                          tGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctgg                                              <  2:4408388‑M1/93‑38 (MQ=255)
                                           gAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctgga                                             >  2:121526‑M1/1‑55 (MQ=255)
                                           gAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctgga                                             <  1:5288135‑M1/93‑39 (MQ=255)
                                              aaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcg                                          <  2:134725‑M1/93‑42 (MQ=255)
                                               aaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgg                                         <  1:112637‑M1/93‑43 (MQ=255)
                                                 gCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggt                                       <  1:3796775‑M1/93‑45 (MQ=255)
                                                    gCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttag                                    <  1:1995780‑M1/93‑48 (MQ=255)
                                                    gCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttag                                    >  2:5389228‑M1/1‑46 (MQ=255)
                                                      tAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgnnnnnaatcgccacaatttcttcgtgagttgcctggatcgggttagta                                  >  2:4010946‑M1/1‑44 (MQ=255)
                                                      tAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagta                                  <  1:7520917‑M1/93‑50 (MQ=255)
                                                      tAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagta                                  >  1:2192201‑M1/1‑44 (MQ=255)
                                                       aaCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaa                                 >  1:332777‑M1/1‑43 (MQ=255)
                                                       aaCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaa                                 >  2:4323216‑M1/1‑43 (MQ=255)
                                                          aaTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgctataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagc                              <  1:1524361‑M1/93‑54 (MQ=255)
                                                          aaTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagc                              <  1:7468456‑M1/93‑54 (MQ=255)
                                                           aTTTGTCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagcc                             >  1:4698213‑M1/1‑39 (MQ=255)
                                                           aTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagcc                             >  2:5191370‑M1/1‑39 (MQ=255)
                                                            tttGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagcca                            <  1:6835048‑M1/93‑56 (MQ=255)
                                                              tGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccaca                          <  2:5736141‑M1/93‑58 (MQ=255)
                                                               gCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacag                         <  1:234482‑M1/93‑59 (MQ=255)
                                                                ccAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacagg                        <  1:7129421‑M1/93‑60 (MQ=255)
                                                                ccAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacagg                        >  2:6973949‑M1/1‑34 (MQ=255)
                                                                ccAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacagg                        >  1:4787892‑M1/1‑34 (MQ=255)
                                                                 cAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggc                       <  2:5584334‑M1/93‑61 (MQ=255)
                                                                 cAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggc                       >  2:5304704‑M1/1‑33 (MQ=255)
                                                                 cAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggc                       <  1:3897556‑M1/93‑61 (MQ=255)
                                                                  aGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggca                      <  1:565987‑M1/93‑62 (MQ=255)
                                                                    ccGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatc                    <  2:2684795‑M1/93‑64 (MQ=255)
                                                                     cGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatct                   <  2:371314‑M1/93‑65 (MQ=255)
                                                                      gTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctt                  >  2:549742‑M1/1‑28 (MQ=255)
                                                                      gTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctt                  >  1:7257531‑M1/1‑28 (MQ=255)
                                                                        tgtgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttc                <  2:1432403‑M1/93‑68 (MQ=255)
                                                                         gtgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttca               <  2:2558639‑M1/93‑69 (MQ=255)
                                                                         gtgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttca               >  1:6994749‑M1/1‑25 (MQ=255)
                                                                         gtgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttca               <  2:1804134‑M1/93‑69 (MQ=255)
                                                                          tgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcag              >  1:6899187‑M1/1‑24 (MQ=255)
                                                                           ggAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagg             <  1:818432‑M1/93‑71 (MQ=255)
                                                                              aaTGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggca          >  2:3013520‑M1/1‑20 (MQ=255)
                                                                                tgatgaTTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt        >  2:1577114‑M1/1‑18 (MQ=255)
                                                                                tgatgaTTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt        <  1:2584265‑M1/93‑76 (MQ=255)
                                                                                tgatgaTTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt        >  2:4125025‑M1/1‑18 (MQ=255)
                                                                                tgatgaTTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt        >  2:7194733‑M1/1‑18 (MQ=255)
                                                                                tgatgaTTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt        >  2:6385440‑M1/1‑18 (MQ=255)
                                                                                 gatgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatta       >  2:4087030‑M1/1‑17 (MQ=255)
                                                                                 gatgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatta       <  2:5931050‑M1/93‑77 (MQ=255)
                                                                                 gatgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatta       <  2:2881296‑M1/93‑77 (MQ=255)
                                                                                  atgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattag      <  1:490894‑M1/93‑78 (MQ=255)
                                                                                   tgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagt     <  1:6628931‑M1/93‑79 (MQ=255)
                                                                                   tgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagt     >  1:1629074‑M1/1‑15 (MQ=255)
                                                                                   tgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagt     >  2:6177989‑M1/1‑15 (MQ=255)
                                                                                    gatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtc    >  1:4522809‑M1/1‑14 (MQ=255)
                                                                                     atTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcg   >  2:1057336‑M1/1‑13 (MQ=255)
                                                                                     atTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcg   <  1:3393081‑M1/93‑81 (MQ=255)
                                                                                     atTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattaagcg   >  1:3632227‑M1/1‑13 (MQ=255)
                                                                                      ttACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcgc  >  2:6092871‑M1/1‑12 (MQ=255)
                                                                                          |                                                                                        
GACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTGCATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAA  >  NC_000913/3754220‑3754398

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