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.

ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTT  >  NC_000913/3815717‑3815915
                                                                                                    |                                                                                                   
acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtc                                                                                                     <  1:797121/100‑1 (MQ=255)
 cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTGCAGACAAAAAAAAGGCGACtca                                                                                                    <  1:1683986/100‑1 (MQ=255)
   ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                  <  2:2475543/100‑1 (MQ=255)
   ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCACCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                  <  2:1672085/100‑1 (MQ=255)
      gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGt                                                                                               >  1:584117/1‑100 (MQ=255)
      gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGt                                                                                               >  2:629088/1‑100 (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)
            tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTAAAAAAAAAAAAAAAAAAGGCGACTGATGAGTCGCCtt                                                                                         <  1:1598172/100‑1 (MQ=255)
             aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttt                                                                                        <  2:2773999/100‑2 (MQ=255)
               ccgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCtttct                                                                                      <  2:1971248/99‑4 (MQ=255)
                cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                     >  2:235378/1‑100 (MQ=255)
                 gtataTTCAATAATCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  2:2525798/96‑6 (MQ=38)
                 gcAAATTCAATAAACTTGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGTCAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  1:1932290/100‑6 (MQ=38)
                 gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  1:2024332/100‑6 (MQ=39)
                 gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  2:1899291/100‑6 (MQ=39)
                 gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  1:1644385/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)
                  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:634770/100‑1 (MQ=255)
                         aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                            <  2:822570/100‑1 (MQ=255)
                           tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                          >  2:2875387/1‑100 (MQ=255)
                           tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                          >  2:2875485/1‑100 (MQ=255)
                           tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                          <  2:264152/100‑1 (MQ=255)
                             aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAg                                                                        >  2:2625157/1‑100 (MQ=255)
                              aCTTGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAACAAAACGCTACTCATCCGTCGCCTTAAAAATCAGTTAGCCAgc                                                                       >  1:394068/1‑100 (MQ=255)
                                   cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGC                                                                  >  2:1298016‑M1/1‑100 (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)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           <  2:2070433‑M1/100‑1 (MQ=255)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           <  1:235376‑M1/100‑1 (MQ=255)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           <  2:530381‑M1/100‑1 (MQ=255)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                            <  1:954973‑M1/99‑1 (MQ=255)
                                            ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                         >  2:299993‑M1/1‑98 (MQ=255)
                                                 cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                    <  2:2702679‑M1/100‑8 (MQ=255)
                                                  aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                   >  1:1454720‑M1/1‑92 (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)
                                                     ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                >  1:229054‑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)
                                                         cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                            <  1:2857163‑M1/100‑16 (MQ=255)
                                                         cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                            <  2:1775312‑M1/100‑16 (MQ=255)
                                                           ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                          >  2:285415‑M1/1‑83 (MQ=255)
                                                           ctcctcATCTTACTTTTCTACAGACAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                          >  2:40937‑M1/1‑83 (MQ=255)
                                                            tcctcATCTTTCGTTTCTACAGACAAAGAAAAGGCGACTCATCAGTTGCCTTAAAAATCTTTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                         <  1:912621‑M1/100‑19 (MQ=255)
                                                                   cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                  >  2:1586925‑M1/1‑75 (MQ=255)
                                                                   cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                  >  2:2695263‑M1/1‑75 (MQ=255)
                                                                   cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                  >  2:2832402‑M1/1‑75 (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)
                                                                       cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcacctgcacttcaatgatgcgcc                                 >  1:2148774‑M1/1‑71 (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                    >  1:658391‑M1/1‑61 (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)
                                                                                     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)
                                                                                                    |                                                                                                   
ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTT  >  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: