Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,819 C→G 17.4% intergenic (‑51/+15) pyrE ← / ← rph orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,8190CG17.4% 357.9 / 65.6 137intergenic (‑51/+15)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base C (65/48);  new base G (7/17);  total (72/65)
Fisher's exact test for biased strand distribution p-value = 1.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 148 total aligned reads displayed.

CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCA  >  NC_000913/3815720‑3815917
                                                                                                    |                                                                                                  
ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                    <  2:2475543/100‑1 (MQ=255)
ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCACCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                    <  2:1672085/100‑1 (MQ=255)
  tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAg                                                                                                  >  2:2395193/1‑100 (MQ=255)
   gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGt                                                                                                 >  1:584117/1‑100 (MQ=255)
   gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGt                                                                                                 >  2:629088/1‑100 (MQ=255)
    cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc                                                                                                <  2:2595790/100‑1 (MQ=255)
      tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGc                                                                                              >  1:114628/1‑100 (MQ=255)
       tccgatctcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                             <  2:2773637/92‑1 (MQ=255)
       gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                             >  2:2074506/1‑100 (MQ=255)
       gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                             >  2:2183198/1‑100 (MQ=255)
        ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCACTCCCCt                                                                                            >  1:729800/1‑100 (MQ=255)
            ccgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCtttct                                                                                        <  2:1971248/99‑4 (MQ=255)
              gtataTTCAATAATCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  2:2525798/96‑6 (MQ=38)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  2:1899291/100‑6 (MQ=39)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  1:1644385/100‑6 (MQ=39)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  1:2024332/100‑6 (MQ=39)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  1:1881841/100‑6 (MQ=38)
               cAAATTCAATAAACTGGCGCTGATATGTTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACATACAAAAAAAAGGCGACTGATGAGTCGCCtttttttt                                                                                     <  1:1190426/100‑7 (MQ=38)
               cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGccagatcgga                                                                                     >  1:2773627/1‑92 (MQ=255)
               cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                     >  1:21965/1‑100 (MQ=255)
               cAAATTCAATAAACTGGCGCTAATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACACCACAACGGCGACTGACGAGTCGCCtttttttt                                                                                     <  1:1883973/100‑7 (MQ=255)
                aaaTTCAATAAACTGTCGCTCATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAAACAAAAAAAAGGCGACTCATCAGCCGCCTTAAAAATc                                                                                    <  1:1793392/100‑1 (MQ=255)
                aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                    <  2:1945230/100‑8 (MQ=38)
                aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                    <  2:757781/100‑8 (MQ=38)
                 aaTTCAAGAAACAGGCGCTGATATGGTCTCATGCCTTAGCTCCTCATCTTAATAAAAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                     <  1:2240115/98‑1 (MQ=255)
                  aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAg                                                                                  <  2:761524/100‑1 (MQ=255)
                   ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGt                                                                                 >  1:896174/1‑100 (MQ=255)
                     cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                               >  2:1877823/1‑100 (MQ=255)
                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                              <  2:822570/100‑1 (MQ=255)
                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                              <  2:634770/100‑1 (MQ=255)
                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                            <  2:264152/100‑1 (MQ=255)
                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                            >  2:2875387/1‑100 (MQ=255)
                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                            >  2:491730/1‑100 (MQ=255)
                          aaCTGGCGCTGATATGGTTTTATGCCGTCGCTCCACGTATTGCTTTTCTACAGACAAAAAAAAGGCGGCTCATGAATCCCCATAAAAACCAGTTTGcc                                                                            >  1:1938911/1‑98 (MQ=255)
                            cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTAGCCTTAAAAATCAGTTTGCCAgcg                                                                        >  2:31517/1‑100 (MQ=255)
                                  cTGATATGGTTTCATGCCTTCGCTCCTCCTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCT                                                                  <  1:2763345‑M1/100‑1 (MQ=255)
                                  cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTCCTACAGACAAAAAAAAGGCTACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCT                                                                  <  1:1584877‑M1/100‑1 (MQ=255)
                                   ttatatggTTTCATGCCTTCGCTCCTCATCTTACTTTACTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTT                                                                 <  1:2405554‑M1/98‑1 (MQ=255)
                                     atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCT                                                               >  1:1009015‑M1/1‑100 (MQ=255)
                                       atGGTTTCATGCGTTCTCCGCTCCTCTTACTCTTCTACAGCCAAAAAAAAGGCGACTCATCAGCCGCCTCACAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                             >  1:1574197/1‑100 (MQ=255)
                                       atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                             <  2:2070433‑M1/100‑1 (MQ=255)
                                       atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                             <  2:530381‑M1/100‑1 (MQ=255)
                                       atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                              <  1:954973‑M1/99‑1 (MQ=255)
                                           tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                             >  1:527795‑M1/1‑96 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAGCCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                      <  2:2755753‑M1/100‑8 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                      <  1:2893866‑M1/100‑8 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                      <  2:2702679‑M1/100‑8 (MQ=255)
                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                     >  1:1454720‑M1/1‑92 (MQ=255)
                                                 gCCTTCGCTTCTCATATTACGTTAAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGccgtcccctgc                                                   <  1:2393323/100‑11 (MQ=255)
                                                 gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                   >  1:2408869‑M1/1‑90 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAACAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCcgtctcctgca                                                  >  1:2599518‑M1/1‑89 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  >  1:935748‑M1/1‑89 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCGACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  <  2:1181313‑M1/100‑12 (MQ=255)
                                                  ccTTCGCTCCTCAGCTTACGTTTCGACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  <  1:1445235‑M1/100‑12 (MQ=255)
                                                   cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                 >  2:1863874‑M1/1‑88 (MQ=255)
                                                   cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                 >  2:2643796‑M1/1‑88 (MQ=255)
                                                      cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                              <  2:1775312‑M1/100‑16 (MQ=255)
                                                        ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                            >  1:164270‑M1/1‑83 (MQ=255)
                                                         tcctcATCTTTCGTTTCTACAGACAAAGAAAAGGCGACTCATCAGTTGCCTTAAAAATCTTTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                           <  1:912621‑M1/100‑19 (MQ=255)
                                                          cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                          <  2:2456259‑M1/100‑20 (MQ=255)
                                                                cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                    >  2:2832402‑M1/1‑75 (MQ=255)
                                                                cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                    >  2:2695263‑M1/1‑75 (MQ=255)
                                                                cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                    >  2:1586925‑M1/1‑75 (MQ=255)
                                                                 ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                   <  2:2905474‑M1/100‑27 (MQ=255)
                                                                  tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                  >  2:1723380‑M1/1‑73 (MQ=255)
                                                                   aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAGATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccg                                 <  1:1334844‑M1/100‑29 (MQ=255)
                                                                   aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccg                                 <  1:155180‑M1/100‑29 (MQ=255)
                                                                    cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGAcgtcccctgcacttcaatgatgcgcccg                                 >  1:2522680‑M1/1‑71 (MQ=255)
                                                                     ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                               >  2:1198873‑M1/1‑70 (MQ=255)
                                                                       ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                             <  2:584118‑M1/100‑33 (MQ=255)
                                                                       ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtatt                             >  2:307805‑M1/1‑68 (MQ=255)
                                                                        tccgatctcaaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                            <  1:91243‑M1/92‑34 (MQ=255)
                                                                         cTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                            <  2:1504860‑M1/100‑34 (MQ=255)
                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                            <  1:2322082‑M1/99‑34 (MQ=255)
                                                                           ncaGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                         >  1:108769‑M1/2‑64 (MQ=255)
                                                                           aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                         <  2:487171‑M1/100‑37 (MQ=255)
                                                                            cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                        >  2:2256599‑M1/1‑63 (MQ=255)
                                                                             aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                        >  1:2698834‑M1/1‑62 (MQ=255)
                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                      >  1:1539039‑M1/1‑61 (MQ=255)
                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                      >  2:2875357‑M1/1‑61 (MQ=255)
                                                                               aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                     >  2:2697983‑M1/1‑60 (MQ=255)
                                                                                cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcagatcgga                    >  2:91244‑M1/1‑59 (MQ=255)
                                                                                  aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                   >  2:1414735‑M1/1‑58 (MQ=255)
                                                                                  aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                   >  1:2642170‑M1/1‑58 (MQ=255)
                                                                                   aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGAGCCGCCTTCTGCcgtcacctgcacttcaatgatgcgcccgtcttcggtcatcact                  >  1:2073412‑M1/1‑57 (MQ=255)
                                                                                    ttttaaggcgACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                 >  1:2317124‑M1/5‑56 (MQ=255)
                                                                                      aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg               >  2:240023‑M1/1‑54 (MQ=255)
                                                                                        aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt             >  2:1332447‑M1/1‑52 (MQ=255)
                                                                                        aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt              >  2:1907434‑M1/1‑52 (MQ=255)
                                                                                        aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt              >  1:2517018‑M1/1‑52 (MQ=255)
                                                                                          ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca           <  2:2894591‑M1/100‑51 (MQ=255)
                                                                                           gCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat          >  2:651171‑M1/1‑49 (MQ=255)
                                                                                            cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg         <  1:2087346‑M1/100‑53 (MQ=255)
                                                                                             gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt        <  2:1009017‑M1/100‑54 (MQ=255)
                                                                                               cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGTCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt        <  1:293364‑M1/98‑54 (MQ=255)
                                                                                                tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg     <  2:2816005‑M1/100‑57 (MQ=255)
                                                                                                tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg     <  1:165611‑M1/100‑57 (MQ=255)
                                                                                                 catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt    >  1:273252‑M1/1‑43 (MQ=255)
                                                                                                 catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt    >  2:1042023‑M1/1‑43 (MQ=255)
                                                                                                   tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct  >  2:50368‑M1/1‑41 (MQ=255)
                                                                                                    |                                                                                                  
CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCA  >  NC_000913/3815720‑3815917

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