Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,816 C→G 15.3% 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,8160CG15.3% 390.4 / 56.0 136intergenic (‑48/+18)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 = 1.00e+00

Only 100 of 151 total aligned reads displayed.

ACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTT  >  NC_000913/3815717‑3815915
                                                                                                   |                                                                                                   
 cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtca                                                                                                    >  1:2789725/1‑100 (MQ=255)
 cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtca                                                                                                    <  2:1746804/100‑1 (MQ=255)
   ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                  <  1:1784154/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                                                                                          >  2:3264255/1‑100 (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:1238807/1‑100 (MQ=255)
                cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                     <  2:622229/100‑1 (MQ=255)
                   aaaTTCAATAAACTGGCTCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                  <  1:2340074/100‑8 (MQ=38)
                   aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                  <  1:2934695/100‑8 (MQ=38)
                   aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc                                                                                  <  2:3403065/100‑1 (MQ=255)
                     aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttgtn                                                                                <  2:71870/100‑10 (MQ=38)
                       tCAATTAACTGGCGCTGAGAGGGTTACATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                  <  1:720064/96‑8 (MQ=25)
                       tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGtt                                                                              >  1:1553595/1‑100 (MQ=255)
                         aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                            <  2:1558550/100‑1 (MQ=255)
                          aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                             >  1:828202/1‑98 (MQ=255)
                          aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                           >  2:2311623/1‑100 (MQ=255)
                          aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                           >  2:3067192/1‑100 (MQ=255)
                          aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                           <  1:1107715/100‑1 (MQ=255)
                             aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAg                                                                        >  2:546840/1‑100 (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                                                                  >  2:111752‑M1/1‑100 (MQ=255)
                                   cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGC                                                                  >  1:8769‑M1/1‑100 (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:2424934‑M1/100‑1 (MQ=255)
                                         tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                            <  2:2425404‑M1/100‑1 (MQ=255)
                                          atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                           >  2:3584288‑M1/1‑100 (MQ=255)
                                            ggTTTAATGCCTTCGCTCCTCATCTTACTTAAAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGccg                                                         <  1:3481225/100‑3 (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)
                                                 cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                    <  2:1344129‑M1/100‑8 (MQ=255)
                                                 cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                    <  2:844711‑M1/100‑8 (MQ=255)
                                                 cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCAACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                    <  1:3873629‑M1/100‑8 (MQ=255)
                                                    gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgc                                                 <  2:212682‑M1/100‑11 (MQ=255)
                                                     ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                   >  1:2550557‑M1/1‑89 (MQ=255)
                                                      cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTCCcgtcccctgcac                                               >  1:2961704‑M1/1‑88 (MQ=255)
                                                      cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                               >  2:1370161‑M1/1‑88 (MQ=255)
                                                          gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttca                                           >  1:1292758‑M1/1‑84 (MQ=255)
                                                           ctcctcATCTTACTTTTCTACAGACAAAAAAAACGCGACTCATCAGTAGCCTTAAAAATCAATTTGCAAGCGCAGACTTCTGCcgtcccctgcacttcaa                                          >  1:370784‑M1/1‑83 (MQ=255)
                                                            tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccttgctcttcaat                                         <  1:2037506‑M1/100‑19 (MQ=255)
                                                            tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                         <  2:1248021‑M1/100‑19 (MQ=255)
                                                             ccccatcTTACTTATCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttc                                            <  1:3667214‑M1/93‑16 (MQ=255)
                                                              ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                       <  2:2558882‑M1/100‑21 (MQ=255)
                                                              ctcATCTAACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                       <  2:3152627‑M1/100‑21 (MQ=255)
                                                               tcATCTTCCTTTTCTACAGCCAAAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                        <  1:2639974‑M1/98‑20 (MQ=255)
                                                                cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                     <  2:679389‑M1/100‑23 (MQ=255)
                                                                cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                     >  2:7177‑M1/1‑78 (MQ=255)
                                                                cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                     >  1:3628905‑M1/1‑78 (MQ=255)
                                                                   cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                  <  2:2182486‑M1/100‑26 (MQ=255)
                                                                    ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                 <  2:276218‑M1/100‑27 (MQ=255)
                                                                    ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                 <  1:3410738‑M1/100‑27 (MQ=255)
                                                                     tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                >  2:3711533‑M1/1‑73 (MQ=255)
                                                                       cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                              >  2:1109298‑M1/1‑71 (MQ=255)
                                                                       cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgt                              >  2:1956510‑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)
                                                                            cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                         >  2:2757542‑M1/1‑66 (MQ=255)
                                                                            cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                         >  2:105297‑M1/1‑66 (MQ=255)
                                                                            cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                         >  2:254808‑M1/1‑66 (MQ=255)
                                                                              aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgttccctgcacttcaatgatgcgcccgtcttcggt                       >  2:3676640‑M1/1‑64 (MQ=255)
                                                                                aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                     <  2:1859382‑M1/100‑39 (MQ=255)
                                                                                aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtca                     <  2:426894‑M1/100‑39 (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:2267806‑M1/1‑59 (MQ=255)
                                                                                   cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                  >  2:1637682‑M1/1‑59 (MQ=255)
                                                                                   cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                   <  2:893612‑M1/99‑41 (MQ=255)
                                                                                    atttttaaggcgacTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtgacctgcacttcaatgatgcgcccgtcttcggtcatcac                 >  1:913921‑M1/7‑58 (MQ=255)
                                                                                    aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                 >  1:1276134‑M1/1‑58 (MQ=255)
                                                                                       aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactac              <  1:3083658‑M1/100‑46 (MQ=255)
                                                                                        ttaaggCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg             >  2:2655662‑M1/3‑54 (MQ=255)
                                                                                        aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg             >  1:28769‑M1/1‑54 (MQ=255)
                                                                                         taagGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt            >  2:2990811‑M1/2‑53 (MQ=255)
                                                                                         taagGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt            >  2:662303‑M1/2‑53 (MQ=255)
                                                                                          aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt           >  1:3817381‑M1/1‑52 (MQ=255)
                                                                                          aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt           >  2:2949233‑M1/1‑52 (MQ=255)
                                                                                           aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc          >  2:937382‑M1/1‑51 (MQ=255)
                                                                                           aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatacgcccgtcttcggtcatcactacgttc          <  1:3660980‑M1/100‑50 (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)
                                                                                                   catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  >  2:1772236‑M1/1‑43 (MQ=255)
                                                                                                   catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  >  1:466147‑M1/1‑43 (MQ=255)
                                                                                                   catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  <  1:1187558‑M1/100‑58 (MQ=255)
                                                                                                   catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtaccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggt  >  1:1483313‑M1/1‑43 (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: