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.

CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCA  >  NC_000913/3815720‑3815917
                                                                                                    |                                                                                                  
ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACtcatc                                                                                                    <  2:2475543/100‑1 (MQ=255)
ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCACCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACtcatc                                                                                                    <  2:1672085/100‑1 (MQ=255)
  tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAg                                                                                                  >  2:2395193/1‑100 (MQ=255)
   gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGt                                                                                                 >  1:584117/1‑100 (MQ=255)
   gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGt                                                                                                 >  2:629088/1‑100 (MQ=255)
    cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTc                                                                                                <  2:2595790/100‑1 (MQ=255)
      tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGc                                                                                              >  1:114628/1‑100 (MQ=255)
       tccgatctcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGcc                                                                                             <  2:2773637/92‑1 (MQ=255)
       gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGcc                                                                                             >  2:2074506/1‑100 (MQ=255)
       gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGcc                                                                                             >  2:2183198/1‑100 (MQ=255)
        ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCACTCCCCt                                                                                            >  1:729800/1‑100 (MQ=255)
            ccgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCtttct                                                                                        <  2:1971248/99‑4 (MQ=255)
              gtataTTCAATAATCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  2:2525798/96‑6 (MQ=38)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  2:1899291/100‑6 (MQ=39)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  1:1644385/100‑6 (MQ=39)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  1:2024332/100‑6 (MQ=39)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                      <  1:1881841/100‑6 (MQ=38)
               cAAATTCAATAAACTGGCGCTGATATGTTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACATAC‑AAAAAAAAGGCGACTGATGAGTCGCCtttttttt                                                                                     <  1:1190426/100‑7 (MQ=38)
               cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGccagatcgga                                                                                     >  1:2773627/1‑92 (MQ=255)
               cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                     >  1:21965/1‑100 (MQ=255)
               cAAATTCAATAAACTGGCGCTAATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑ACCACAACGGCGACTGACGAGTCGCCtttttttt                                                                                     <  1:1883973/100‑7 (MQ=255)
                aaaTTCAATAAACTGTCGCTCATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAAAC‑AAAAAAAAGGCGACTCATCAGCCGCCTTAAAAATc                                                                                    <  1:1793392/100‑1 (MQ=255)
                aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                    <  2:1945230/100‑8 (MQ=38)
                aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                    <  2:757781/100‑8 (MQ=38)
                 aaTTCAAGAAACAGGCGCTGATATGGTCTCATGCCTTAGCTCCTCATCTTAATAAAAAAAAAAA‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                     <  1:2240115/98‑1 (MQ=255)
                  aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAg                                                                                  <  2:761524/100‑1 (MQ=255)
                   ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGt                                                                                 >  1:896174/1‑100 (MQ=255)
                     cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                               >  2:1877823/1‑100 (MQ=255)
                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                              <  2:822570/100‑1 (MQ=255)
                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                              <  2:634770/100‑1 (MQ=255)
                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                            <  2:264152/100‑1 (MQ=255)
                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                            >  2:2875387/1‑100 (MQ=255)
                        tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                            >  2:491730/1‑100 (MQ=255)
                          aaCTGGCGCTGATATGGTTTTATGCCGTCGCTCCACGTATTGCTTTTCTACAGAC‑AAAAAAAAGGCGGCTCATGAATCCCCATAAAAACCAGTTTGcc                                                                            >  1:1938911/1‑98 (MQ=255)
                            cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTAGCCTTAAAAATCAGTTTGCCAgcg                                                                        >  2:31517/1‑100 (MQ=255)
                                  cTGATATGGTTTCATGCCTTCGCTCCTCCTCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCT                                                                  <  1:2763345‑M1/100‑1 (MQ=255)
                                  cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTCCTACAGAC‑AAAAAAAAGGCTACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCT                                                                  <  1:1584877‑M1/100‑1 (MQ=255)
                                   ttatatggTTTCATGCCTTCGCTCCTCATCTTACTTTACTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTT                                                                 <  1:2405554‑M1/98‑1 (MQ=255)
                                     atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCT                                                               >  1:1009015‑M1/1‑100 (MQ=255)
                                       atGGTTTCATGCGTTCTCCGCTCCTCTTACTCTTCTACAGCC‑AAAAAAAAGGCGACTCATCAGCCGCCTCACAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                             >  1:1574197/1‑100 (MQ=255)
                                       atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                             <  2:2070433‑M1/100‑1 (MQ=255)
                                       atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                             <  2:530381‑M1/100‑1 (MQ=255)
                                       atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                              <  1:954973‑M1/99‑1 (MQ=255)
                                           tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                             >  1:527795‑M1/1‑96 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAGCCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                      <  2:2755753‑M1/100‑8 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                      <  1:2893866‑M1/100‑8 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                      <  2:2702679‑M1/100‑8 (MQ=255)
                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                     >  1:1454720‑M1/1‑92 (MQ=255)
                                                 gCCTTCGCTTCTCATATTACGTTAAAAAAAAA‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGccgtcccctgc                                                   <  1:2393323/100‑11 (MQ=255)
                                                 gCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                   >  1:2408869‑M1/1‑90 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAACAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCcgtctcctgca                                                  >  1:2599518‑M1/1‑89 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  >  1:935748‑M1/1‑89 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCGACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  <  2:1181313‑M1/100‑12 (MQ=255)
                                                  ccTTCGCTCCTCAGCTTACGTTTCGACAGAC‑AAAAAAAAGGCGGCTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  <  1:1445235‑M1/100‑12 (MQ=255)
                                                   cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                 >  2:1863874‑M1/1‑88 (MQ=255)
                                                   cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                 >  2:2643796‑M1/1‑88 (MQ=255)
                                                      cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                              <  2:1775312‑M1/100‑16 (MQ=255)
                                                        ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                            >  1:164270‑M1/1‑83 (MQ=255)
                                                         tcctcATCTTTCGTTTCTACAGAC‑AAAGAAAAGGCGACTCATCAGTTGCCTTAAAAATCTTTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                           <  1:912621‑M1/100‑19 (MQ=255)
                                                          cctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                          <  2:2456259‑M1/100‑20 (MQ=255)
                                                                cTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                    >  2:2832402‑M1/1‑75 (MQ=255)
                                                                cTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                    >  2:2695263‑M1/1‑75 (MQ=255)
                                                                cTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                    >  2:1586925‑M1/1‑75 (MQ=255)
                                                                 ttACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                   <  2:2905474‑M1/100‑27 (MQ=255)
                                                                  tACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                  >  2:1723380‑M1/1‑73 (MQ=255)
                                                                   aCTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAGATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccg                                 <  1:1334844‑M1/100‑29 (MQ=255)
                                                                   aCTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccg                                 <  1:155180‑M1/100‑29 (MQ=255)
                                                                    cTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGAcgtcccctgcacttcaatgatgcgcccg                                 >  1:2522680‑M1/1‑71 (MQ=255)
                                                                     ttttCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                               >  2:1198873‑M1/1‑70 (MQ=255)
                                                                       ttCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                             <  2:584118‑M1/100‑33 (MQ=255)
                                                                       ttCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtatt                             >  2:307805‑M1/1‑68 (MQ=255)
                                                                        tccgatctc‑aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                            <  1:91243‑M1/92‑34 (MQ=255)
                                                                         cTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                            <  2:1504860‑M1/100‑34 (MQ=255)
                                                                         cTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                            <  1:2322082‑M1/99‑34 (MQ=255)
                                                                           ncaGAC‑AAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                         >  1:108769‑M1/2‑64 (MQ=255)
                                                                           aCAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                         <  2:487171‑M1/100‑37 (MQ=255)
                                                                            cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                        >  2:2256599‑M1/1‑63 (MQ=255)
                                                                             aGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                        >  1:2698834‑M1/1‑62 (MQ=255)
                                                                              gAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                      >  1:1539039‑M1/1‑61 (MQ=255)
                                                                              gAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                      >  2:2875357‑M1/1‑61 (MQ=255)
                                                                               aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                     >  2:2697983‑M1/1‑60 (MQ=255)
                                                                                c‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcagatcgga                    >  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)
                                                                                                    |                                                                                                  
CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCA  >  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: ‑