Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,816 C→G 17.8% 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,8160CG17.8% 295.1 / 59.4 119intergenic (‑48/+18)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base C (57/40);  new base G (9/12);  total (67/52)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

Only 100 of 128 total aligned reads displayed.

CACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTT  >  NC_000913/3815718‑3815915
                                                                                                  |                                                                                                   
cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtca                                                                                                    >  2:347767/1‑100 (MQ=255)
  ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                  <  2:1623470/100‑1 (MQ=255)
   cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatca                                                                                                 <  2:1713477/100‑1 (MQ=255)
   cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCAGCTTACTTTTCTACAGACAAAAAAAAGGCGCCtcatca                                                                                                 >  1:563356/1‑100 (MQ=255)
     gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCCTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGt                                                                                               >  2:125557/1‑100 (MQ=255)
      cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc                                                                                              >  2:1687773/1‑100 (MQ=255)
         gctAAGCGCAAATTCATTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                           <  1:850202/100‑1 (MQ=255)
          ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCt                                                                                          >  2:2315131/1‑100 (MQ=255)
          ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCATCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCt                                                                                          <  2:132145/100‑1 (MQ=255)
           tAAGCGCAAATTTAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCtt                                                                                         <  1:2173925/100‑1 (MQ=255)
            aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTa                                                                                        >  2:943316/1‑100 (MQ=255)
            aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTa                                                                                        <  1:1705082/100‑1 (MQ=255)
             aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCtttt                                                                                       <  2:2228651/100‑3 (MQ=255)
             aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaa                                                                                       >  2:1417942/1‑100 (MQ=255)
              gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaa                                                                                      >  2:1562464/1‑100 (MQ=255)
               cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                     <  2:560321/100‑1 (MQ=255)
               cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTAATTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                     <  1:2568191/100‑1 (MQ=255)
               cgcAAATTCAATAAACTCGCGCTTATATGGTTTAATGCCTTCGCTCCTCATCTTACTTCGCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttt                                                                                        <  1:1446836/97‑2 (MQ=255)
                gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                    <  1:2132745/100‑6 (MQ=39)
                gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaaa                                                                                    <  2:2186062/100‑1 (MQ=255)
                gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaaa                                                                                    >  2:1591328/1‑100 (MQ=255)
                 cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCtttttttt                                                                                   <  2:2198314/100‑7 (MQ=38)
                  aaaTTCAATAAACTGGCTCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCtttttttt                                                                                   <  1:1346171/99‑7 (MQ=38)
                  aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                  <  2:1554597/100‑8 (MQ=38)
                  aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc                                                                                  <  2:2893197/100‑1 (MQ=255)
                   aaTTCAACAAACTGGCGCTGATATGGTTTCATGCATTCGCTCCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCtttttttt                                                                                   <  1:1952859/98‑7 (MQ=37)
                    aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAg                                                                                <  2:728746/100‑1 (MQ=255)
                     ttCAATAAGCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                  <  1:2500499/97‑8 (MQ=37)
                       cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                             >  1:1093096/1‑100 (MQ=255)
                       cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATACCGACAAAAAAAAGGCCACTCATCAGTCGCCTTAAAAACCAGttt                                                                             >  1:1664508/1‑100 (MQ=255)
                        aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                            >  1:611113/1‑100 (MQ=255)
                         aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                           >  2:893158/1‑100 (MQ=255)
                         aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                           >  2:2644153/1‑100 (MQ=255)
                               tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcgc                                                                     <  2:2540762/100‑1 (MQ=255)
                                   gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCC                                                                 <  1:1812895‑M1/100‑1 (MQ=255)
                                    cTGATATGGTTTCATGCCTTCGCTGCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCTCCTTAAAAATCAGTGTGCCAGCGCCGC                                                                  <  1:836092‑M1/98‑1 (MQ=255)
                                    cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCT                                                                >  2:2859320‑M1/1‑100 (MQ=255)
                                     tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAAACGACTTAAACCTCAGTTTGCAAGCGCCGCCTT                                                               >  2:1696802‑M1/1‑100 (MQ=255)
                                      gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTC                                                              >  2:58754‑M1/1‑100 (MQ=255)
                                      gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTC                                                              >  2:2490059‑M1/1‑100 (MQ=255)
                                        tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                            >  2:2813788‑M1/1‑100 (MQ=255)
                                        tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                            <  1:2644137‑M1/100‑1 (MQ=255)
                                        tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTAATCAGTCGCCTTAAAAATCAGTTTGCAAGCGCCGCCTTCTG                                                            >  2:1696834‑M1/1‑100 (MQ=255)
                                         atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATACAGACAAAAAAAAGGCGACTCATCAGTCGAATTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           <  1:369504‑M1/100‑1 (MQ=255)
                                          tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACACCAAAAACGCGACTCATCAGTCCCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCtccc                                                          >  1:892688/1‑97 (MQ=255)
                                          tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                          >  2:2289197‑M1/1‑99 (MQ=255)
                                          tGGTTTCATGCCTTCGCTCCTCAGCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           <  1:543648‑M1/99‑1 (MQ=255)
                                          tGGTTGGATGCCGTCGCTCCACATCTAAATTAACAAAAGAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                           <  2:2006657/99‑1 (MQ=255)
                                                cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                    <  2:1864097‑M1/100‑8 (MQ=255)
                                                cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgagatc                                                    >  1:432373‑M1/1‑93 (MQ=255)
                                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                   >  1:2695996‑M1/1‑92 (MQ=255)
                                                  tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                  >  2:1611595‑M1/1‑91 (MQ=255)
                                                  tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                   >  1:686257‑M1/1‑91 (MQ=255)
                                                  tGCCTTCGCTCCTCATCTTACTTTTCTAAAGAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                  <  1:2834394‑M1/100‑10 (MQ=255)
                                                   gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                 <  2:313338‑M1/100‑11 (MQ=255)
                                                   gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                 <  2:611116‑M1/100‑11 (MQ=255)
                                                    ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCGGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                >  2:2912343‑M1/1‑89 (MQ=255)
                                                    ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                <  2:497275‑M1/100‑12 (MQ=255)
                                                    ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                >  2:949251‑M1/1‑89 (MQ=255)
                                                    ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                >  2:1250724‑M1/1‑89 (MQ=255)
                                                     cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcttcccctgcac                                               >  1:2435933‑M1/1‑88 (MQ=255)
                                                     cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                               >  2:410562‑M1/1‑88 (MQ=255)
                                                                aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgc                                    >  2:1861114‑M1/1‑77 (MQ=255)
                                                                     aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccg                               >  2:541012‑M1/1‑72 (MQ=255)
                                                                        tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGGCTTCTGCcgtcccctgcacttcaatgatgcgtccgtct                            >  1:1587913‑M1/1‑69 (MQ=255)
                                                                           cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                         >  2:448978‑M1/1‑66 (MQ=255)
                                                                           cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                         <  1:963573‑M1/100‑35 (MQ=255)
                                                                           cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                          <  2:132434‑M1/99‑34 (MQ=255)
                                                                             aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                       <  1:1820701‑M1/100‑37 (MQ=255)
                                                                             aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                        <  2:1343938‑M1/99‑36 (MQ=255)
                                                                                gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                    >  2:2446586‑M1/1‑61 (MQ=255)
                                                                                gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                    >  2:44307‑M1/1‑61 (MQ=255)
                                                                                gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                      <  1:1076508‑M1/98‑38 (MQ=255)
                                                                                 aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCGGTCTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcagc                   <  1:844309‑M1/100‑41 (MQ=255)
                                                                                 aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                   <  1:1361033‑M1/100‑41 (MQ=255)
                                                                                  cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                  >  2:1573117‑M1/1‑59 (MQ=255)
                                                                                   aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                 >  1:937774‑M1/1‑58 (MQ=255)
                                                                                   aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                 <  1:1865255‑M1/100‑43 (MQ=255)
                                                                                   aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                 >  2:2399801‑M1/1‑58 (MQ=255)
                                                                                    aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                >  1:745096‑M1/1‑57 (MQ=255)
                                                                                    aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                >  1:1200761‑M1/1‑57 (MQ=255)
                                                                                    aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                 >  1:721880‑M1/1‑57 (MQ=255)
                                                                                     aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta               >  2:587204‑M1/1‑56 (MQ=255)
                                                                                      tttaaggcGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactac              >  1:869817‑M1/4‑55 (MQ=255)
                                                                                         aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcccttcattgatgcgcccgtcttcggtcatcactacgtt           >  1:1963197‑M1/1‑52 (MQ=255)
                                                                                         aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt           >  1:1035362‑M1/1‑52 (MQ=255)
                                                                                         aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt           >  2:1900302‑M1/1‑52 (MQ=255)
                                                                                         aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt           >  2:67263‑M1/1‑52 (MQ=255)
                                                                                         aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCACCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg             >  1:197073‑M1/1‑52 (MQ=255)
                                                                                         aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttctatgatgcgcccgtcttcggtcatcactacgtt           <  2:2767354‑M1/100‑49 (MQ=255)
                                                                                           ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCGTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca         <  1:1865541‑M1/100‑51 (MQ=255)
                                                                                           ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcactgctgcgcccgtcttcggtcatcactacgttca         <  1:1913683‑M1/100‑51 (MQ=255)
                                                                                           ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca         <  1:73753‑M1/100‑51 (MQ=255)
                                                                                             cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg       <  2:2545614‑M1/100‑53 (MQ=255)
                                                                                              gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt      <  2:118654‑M1/100‑54 (MQ=255)
                                                                                                 tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg   <  2:289743‑M1/100‑57 (MQ=255)
                                                                                                 tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg   <  1:1003947‑M1/100‑57 (MQ=255)
                                                                                                  catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  >  2:2749548‑M1/1‑43 (MQ=255)
                                                                                                  catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  >  2:767422‑M1/1‑43 (MQ=255)
                                                                                                  catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  >  2:1124279‑M1/1‑43 (MQ=255)
                                                                                                  |                                                                                                   
CACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTT  >  NC_000913/3815718‑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: