Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,819 C→G 17.8% 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,8190CG17.8% 302.3 / 58.2 119intergenic (‑51/+15)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 = 8.85e-01

Only 100 of 128 total aligned reads displayed.

CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCA  >  NC_000913/3815720‑3815917
                                                                                                   |                                                                                                  
ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                    <  2:1623470/100‑1 (MQ=255)
 cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCAGCTTACTTTTCTACAGACAAAAAAAAGGCGCCtcatca                                                                                                   >  1:563356/1‑100 (MQ=255)
    cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc                                                                                                >  2:1688502/1‑100 (MQ=255)
       gctAAGCGCAAATTCATTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                             <  1:850202/100‑1 (MQ=255)
        ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCt                                                                                            <  2:1814722/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)
           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:1591328/1‑100 (MQ=255)
              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaaa                                                                                      <  2:2186062/100‑1 (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)
                          aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAg                                                                          >  2:1281371/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                                                              <  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)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                      <  2:1864097‑M1/100‑8 (MQ=255)
                                              cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgagatc                                                      >  1:432373‑M1/1‑93 (MQ=255)
                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccct                                                     >  2:2595091‑M1/1‑92 (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:2913516‑M1/1‑89 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  >  2:1250724‑M1/1‑89 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  >  2:949251‑M1/1‑89 (MQ=255)
                                                  ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                  <  2:497275‑M1/100‑12 (MQ=255)
                                                   cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcttcccctgcac                                                 >  1:2435933‑M1/1‑88 (MQ=255)
                                                   cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                 >  2:410562‑M1/1‑88 (MQ=255)
                                                        ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                            <  2:139993‑M1/100‑18 (MQ=255)
                                                          cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                          >  2:2261446‑M1/1‑81 (MQ=255)
                                                              aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgc                                      >  2:1861114‑M1/1‑77 (MQ=255)
                                                                   aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccg                                 >  2:541012‑M1/1‑72 (MQ=255)
                                                                     ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtc                               >  2:1591240‑M1/1‑70 (MQ=255)
                                                                      tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGGCTTCTGCcgtcccctgcacttcaatgatgcgtccgtct                              >  1:1587913‑M1/1‑69 (MQ=255)
                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                           >  2:2067234‑M1/1‑66 (MQ=255)
                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcg                           >  2:448978‑M1/1‑66 (MQ=255)
                                                                         cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                            <  2:132434‑M1/99‑34 (MQ=255)
                                                                           aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggt                         <  1:1820701‑M1/100‑37 (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)
                                                                                gatttttaaggcgactGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                    >  2:1024269‑M1/8‑59 (MQ=255)
                                                                                cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                    >  2:2475304‑M1/1‑59 (MQ=255)
                                                                                 aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                   >  1:937774‑M1/1‑58 (MQ=255)
                                                                                 aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                   >  2:2399801‑M1/1‑58 (MQ=255)
                                                                                 aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcac                   <  1:1865255‑M1/100‑43 (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             >  2:67263‑M1/1‑52 (MQ=255)
                                                                                       aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt             >  2:1900302‑M1/1‑52 (MQ=255)
                                                                                       aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt             >  1:1035362‑M1/1‑52 (MQ=255)
                                                                                       aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCACCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg               >  1:197073‑M1/1‑52 (MQ=255)
                                                                                       aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttctatgatgcgcccgtcttcggtcatcactacgtt             <  2:2767354‑M1/100‑49 (MQ=255)
                                                                                        aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc            >  2:2834944‑M1/1‑51 (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)
                                                                                            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)
                                                                                                  tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct  >  1:90281‑M1/1‑41 (MQ=255)
                                                                                                   |                                                                                                  
CCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCA  >  NC_000913/3815720‑3815917

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: