Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,816 C→G 16.2% intergenic (‑48/+18) 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,8160CG16.2% 368.2 / 58.4 136intergenic (‑48/+18)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base C (66/48);  new base G (5/17);  total (71/65)
Fisher's exact test for biased strand distribution p-value = 4.32e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 147 total aligned reads displayed.

ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTT  >  NC_000913/3815717‑3815915
                                                                                                    |                                                                                                   
acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACtc                                                                                                     <  1:797121/100‑1 (MQ=255)
 cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTGCAGAC‑AAAAAAAAGGCGACtca                                                                                                    <  1:1683986/100‑1 (MQ=255)
   ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACtcatc                                                                                                  <  2:2475543/100‑1 (MQ=255)
   ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCACCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACtcatc                                                                                                  <  2:1672085/100‑1 (MQ=255)
      gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGt                                                                                               >  1:584117/1‑100 (MQ=255)
      gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGt                                                                                               >  2:629088/1‑100 (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)
            tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTAAAAAAAAAA‑AAAAAAAAGGCGACTGATGAGTCGCCtt                                                                                         <  1:1598172/100‑1 (MQ=255)
             aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttt                                                                                        <  2:2773999/100‑2 (MQ=255)
               ccgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCtttct                                                                                      <  2:1971248/99‑4 (MQ=255)
                cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                     >  2:235378/1‑100 (MQ=255)
                 gtataTTCAATAATCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  2:2525798/96‑6 (MQ=38)
                 gcAAATTCAATAAACTTGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGTC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  1:1932290/100‑6 (MQ=38)
                 gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  1:2024332/100‑6 (MQ=39)
                 gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  2:1899291/100‑6 (MQ=39)
                 gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  1:1644385/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)
                  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:634770/100‑1 (MQ=255)
                         aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                            <  2:822570/100‑1 (MQ=255)
                           tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                          >  2:2875387/1‑100 (MQ=255)
                           tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                          >  2:2875485/1‑100 (MQ=255)
                           tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                          <  2:264152/100‑1 (MQ=255)
                             aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAg                                                                        >  2:2625157/1‑100 (MQ=255)
                              aCTTGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAACAAAACGCTACTCATCCGTCGCCTTAAAAATCAGTTAGCCAgc                                                                       >  1:394068/1‑100 (MQ=255)
                                   cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGC                                                                  >  2:1298016‑M1/1‑100 (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)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           <  2:2070433‑M1/100‑1 (MQ=255)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           <  1:235376‑M1/100‑1 (MQ=255)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           <  2:530381‑M1/100‑1 (MQ=255)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                            <  1:954973‑M1/99‑1 (MQ=255)
                                            ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                         >  2:299993‑M1/1‑98 (MQ=255)
                                                 cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                    <  2:2702679‑M1/100‑8 (MQ=255)
                                                  aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                   >  1:1454720‑M1/1‑92 (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)
                                                     ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                >  1:229054‑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)
                                                         cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                            <  1:2857163‑M1/100‑16 (MQ=255)
                                                         cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                            <  2:1775312‑M1/100‑16 (MQ=255)
                                                           ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                          >  2:285415‑M1/1‑83 (MQ=255)
                                                           ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                          >  2:40937‑M1/1‑83 (MQ=255)
                                                            tcctcATCTTTCGTTTCTACAGAC‑AAAGAAAAGGCGACTCATCAGTTGCCTTAAAAATCTTTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                         <  1:912621‑M1/100‑19 (MQ=255)
                                                                   cTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                  >  2:1586925‑M1/1‑75 (MQ=255)
                                                                   cTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                  >  2:2695263‑M1/1‑75 (MQ=255)
                                                                   cTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                  >  2:2832402‑M1/1‑75 (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‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcacctgcacttcaatgatgcgcc                                 >  1:2148774‑M1/1‑71 (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                    >  1:658391‑M1/1‑61 (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)
                                                                                     aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                 >  2:496988‑M1/1‑58 (MQ=255)
                                                                                      aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                >  1:2051214‑M1/1‑57 (MQ=255)
                                                                                         aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg             >  2:240023‑M1/1‑54 (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)
                                                                                              gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat        <  1:1654295‑M1/100‑52 (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)
                                                                                                  cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg    <  2:2404031‑M1/100‑56 (MQ=255)
                                                                                                  cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg    >  2:2726388‑M1/1‑45 (MQ=255)
                                                                                                  cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc     <  2:2675778‑M1/99‑55 (MQ=255)
                                                                                                    catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  >  1:273252‑M1/1‑43 (MQ=255)
                                                                                                    catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  >  2:1042023‑M1/1‑43 (MQ=255)
                                                                                                    catcaGTCGCATTAAAAATCAGTTAGCCAGCGGCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcg    <  1:2732802‑M1/98‑56 (MQ=255)
                                                                                                    |                                                                                                   
ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTT  >  NC_000913/3815717‑3815915

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: ‑