Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,819 C→G 15.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,8190CG15.4% 384.1 / 57.0 136intergenic (‑51/+15)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base C (60/55);  new base G (15/6);  total (75/61)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 151 total aligned reads displayed.

CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAT  >  NC_000913/3815720‑3815918
                                                                                                   |                                                                                                   
ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                     <  1:1784154/100‑1 (MQ=255)
ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTAGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                     <  1:3678984/100‑1 (MQ=255)
ccTGCTTGCTAAGCGCAAATGCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                     >  1:483272/1‑100 (MQ=255)
 cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatca                                                                                                    >  2:3253421/1‑100 (MQ=255)
  tgcttgctAAGCGAAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAAACAAAAAAAAGGCGACtcatc                                                                                                     <  1:3488851/98‑1 (MQ=255)
    cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTc                                                                                                 <  1:157681/100‑1 (MQ=255)
        ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCt                                                                                             >  1:744157/1‑100 (MQ=255)
          aaGCGAAAATTCAATAAACTGGCGCTGATATGGTTTAATTCCTTCGCTCCTCATCTTAATAAAAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGcc                                                                                              <  1:3763855/97‑1 (MQ=255)
            gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttt                                                                                         <  2:2074615/100‑4 (MQ=255)
            gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaa                                                                                         <  1:2628577/100‑1 (MQ=255)
             cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                        <  2:622229/100‑1 (MQ=255)
             cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                        >  2:1238807/1‑100 (MQ=255)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaaa                                                                                       <  2:34114/100‑1 (MQ=255)
               cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                      >  2:1311393/1‑100 (MQ=255)
               cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                      >  1:745337/1‑100 (MQ=255)
                aaaTTCAATAAACTGGCTCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                     <  1:2340074/100‑8 (MQ=38)
                 aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                                    <  2:165066/100‑1 (MQ=255)
                  aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttgtn                                                                                   <  2:71870/100‑10 (MQ=38)
                    tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGtt                                                                                 >  1:1553595/1‑100 (MQ=255)
                     cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                                >  1:2933804/1‑100 (MQ=255)
                      aagaatcTGGCGCTGATATGGTTTAAGGCCTTCGATCCTCATCTTACTTAAAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                               <  1:573920/94‑1 (MQ=255)
                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                               <  2:1558550/100‑1 (MQ=255)
                      aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                               >  1:3548216/1‑100 (MQ=255)
                       aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                              >  2:2311623/1‑100 (MQ=255)
                       aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                              >  2:3067192/1‑100 (MQ=255)
                       aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                              <  1:1107715/100‑1 (MQ=255)
                         aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCa                                                                            <  2:1176992/100‑1 (MQ=255)
                            cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAACAAAAGGCGACTCATCAGTAGCCTTAAAAATCAGATTGCCAgc                                                                          >  1:963510/1‑99 (MQ=255)
                              ggcggagATATGGTTTCATGCCTTCGCTCCCCATCTTACCTTGCTACAGACAAAAAAAAGGCGGCTCATCAAAAAACAAAAAAATCAGTTTGCCAgcg                                                                         <  1:1779437/92‑1 (MQ=255)
                              ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTGACTTTTCTACAGACAAAACAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCC                                                                       >  1:1416684‑M1/1‑100 (MQ=255)
                                cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCac                                                                     >  1:1683225‑M1/1‑98 (MQ=255)
                                cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGC                                                                     >  1:8769‑M1/1‑100 (MQ=255)
                                cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCG                                                                      >  1:893604‑M1/1‑99 (MQ=255)
                                cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTAAACAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                         <  1:16418/96‑1 (MQ=255)
                                    gctatggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTC                                                                 >  2:2619814‑M1/3‑100 (MQ=255)
                                    gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTC                                                                 <  1:3449792‑M1/100‑1 (MQ=255)
                                     atatGGTTTCATGCCTTCGCGCCTCACCTTACTTTTCTACAGACAAAAAAAAGGCGACGCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCT                                                                <  1:3639444‑M1/100‑1 (MQ=255)
                                      tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                               <  2:2425404‑M1/100‑1 (MQ=255)
                                      tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                               <  2:2424934‑M1/100‑1 (MQ=255)
                                       atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                              >  2:3584288‑M1/1‑100 (MQ=255)
                                       atGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                              <  1:3764925‑M1/100‑1 (MQ=255)
                                        tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                             >  2:626599‑M1/1‑99 (MQ=255)
                                          gTTTCATGCGTTCGCTCCTCTTCTTACTATACAAAAAACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                            <  1:2534408‑M1/99‑3 (MQ=255)
                                          gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                           >  1:1980802‑M1/1‑97 (MQ=255)
                                          ccgatctgccttcGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                           <  2:3384031‑M1/94‑4 (MQ=255)
                                           tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                          >  1:2992294‑M1/1‑96 (MQ=255)
                                             tCATGCCGTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                         <  1:3643532‑M1/99‑6 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                       <  2:1344129‑M1/100‑8 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCAACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                       <  1:3873629‑M1/100‑8 (MQ=255)
                                                 gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                    <  2:212682‑M1/100‑11 (MQ=255)
                                                 gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtagatcgg                                                    >  1:3384004‑M1/1‑90 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                      >  1:2550557‑M1/1‑89 (MQ=255)
                                                   cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTCCcgtcccctgcac                                                  >  1:2961704‑M1/1‑88 (MQ=255)
                                                   cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                  <  2:515353‑M1/100‑13 (MQ=255)
                                                       gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttca                                              >  1:1292758‑M1/1‑84 (MQ=255)
                                                         tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccttgctcttcaat                                            <  1:2037506‑M1/100‑19 (MQ=255)
                                                          ccccatcTTACTTATCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                               <  1:3667214‑M1/93‑16 (MQ=255)
                                                           ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                          <  2:2558882‑M1/100‑21 (MQ=255)
                                                            tcATCTTCCTTTTCTACAGCCAAAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                           <  1:2639974‑M1/98‑20 (MQ=255)
                                                             cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                        <  2:679389‑M1/100‑23 (MQ=255)
                                                                cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                     <  2:2182486‑M1/100‑26 (MQ=255)
                                                                 ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                    <  2:3701789‑M1/100‑27 (MQ=255)
                                                                 ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                    <  2:276218‑M1/100‑27 (MQ=255)
                                                                    cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                 >  2:1956510‑M1/1‑71 (MQ=255)
                                                                    cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                                 >  2:1109298‑M1/1‑71 (MQ=255)
                                                                     ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                                >  1:2797972‑M1/1‑70 (MQ=255)
                                                                      tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtct                               >  2:1940382‑M1/1‑69 (MQ=255)
                                                                       ttCTACAGACAAAATAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                              <  1:2438087‑M1/100‑33 (MQ=255)
                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcgcctgcacttcaatgaggcgcccgtcttcg                            <  1:392676‑M1/100‑35 (MQ=255)
                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                            >  2:105297‑M1/1‑66 (MQ=255)
                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                            >  2:254808‑M1/1‑66 (MQ=255)
                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                            >  2:2757542‑M1/1‑66 (MQ=255)
                                                                           aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgttccctgcacttcaatgatgcgcccgtcttcggt                          >  2:3676640‑M1/1‑64 (MQ=255)
                                                                            cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                         <  2:779215‑M1/100‑38 (MQ=255)
                                                                             aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                        <  2:1859382‑M1/100‑39 (MQ=255)
                                                                             aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                        <  2:426894‑M1/100‑39 (MQ=255)
                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                       >  1:2808598‑M1/1‑61 (MQ=255)
                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                        >  2:1837723‑M1/1‑61 (MQ=255)
                                                                               aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                      >  1:2231005‑M1/1‑60 (MQ=255)
                                                                               aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                      <  1:3192064‑M1/100‑41 (MQ=255)
                                                                                cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                     >  2:303787‑M1/1‑59 (MQ=255)
                                                                                cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                     >  2:3373137‑M1/1‑59 (MQ=255)
                                                                                cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                     >  2:2267806‑M1/1‑59 (MQ=255)
                                                                                cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                      <  2:893612‑M1/99‑41 (MQ=255)
                                                                                 aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                    >  1:1276134‑M1/1‑58 (MQ=255)
                                                                                    tttaaggcGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactac                 >  1:3453469‑M1/4‑55 (MQ=255)
                                                                                     ttaaggCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                >  2:2655662‑M1/3‑54 (MQ=255)
                                                                                      taagGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt               >  2:1447605‑M1/2‑53 (MQ=255)
                                                                                       aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt              >  2:2949233‑M1/1‑52 (MQ=255)
                                                                                       aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt              >  2:900704‑M1/1‑52 (MQ=255)
                                                                                        aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc             >  2:697870‑M1/1‑51 (MQ=255)
                                                                                        aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc             >  2:937382‑M1/1‑51 (MQ=255)
                                                                                         ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca            >  2:2288471‑M1/1‑50 (MQ=255)
                                                                                         ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca            <  1:2888343‑M1/100‑51 (MQ=255)
                                                                                              cTCATCAGTCGCCTTAAAAATCAGTTTGCCATCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc        <  1:2455513‑M1/99‑55 (MQ=255)
                                                                                               tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg      <  2:3225573‑M1/100‑57 (MQ=255)
                                                                                                catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt     >  2:1772236‑M1/1‑43 (MQ=255)
                                                                                                catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt     <  1:1187558‑M1/100‑58 (MQ=255)
                                                                                                catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtaccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt     >  1:1483313‑M1/1‑43 (MQ=255)
                                                                                                   caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc  <  2:28769‑M1/100‑61 (MQ=255)
                                                                                                   |                                                                                                   
CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAT  >  NC_000913/3815720‑3815918

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: