Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,816 C→G 5.6% 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,8160CG 5.6% 394.2 / 12.2 125intergenic (‑48/+18)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base C (62/56);  new base G (3/4);  total (65/60)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 131 total aligned reads displayed.

AAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGC  >  NC_000913/3815668‑3815965
                                                                                                                                                    |                                                                                                                                                     
aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtca                                                                                                                                                      <  2:650658/150‑1 (MQ=255)
 aaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcat                                                                                                                                                     >  2:487926/1‑150 (MQ=255)
  aTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCTTCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                                                                                    <  2:154712/150‑1 (MQ=255)
    aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAg                                                                                                                                                  >  1:559971/1‑150 (MQ=255)
       ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCg                                                                                                                                               >  2:100819/1‑150 (MQ=255)
        gCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGc                                                                                                                                              <  1:80846/150‑1 (MQ=255)
         cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATAGGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                                                                             >  1:926082/1‑150 (MQ=255)
              gcgcCCGGTTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAACCTGGCGCTGATATGGTTTCATGCCTCCGCTCCTCATCTTACTTTTCAACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTaaa                                                                                                                                        <  1:1570727/150‑1 (MQ=255)
               cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                                                                       >  2:1844206/1‑150 (MQ=255)
                 cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                                                                     >  2:716160/1‑150 (MQ=255)
                  ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                                                                    <  1:1773558/150‑8 (MQ=255)
                  ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc                                                                                                                                    >  1:1103044/1‑150 (MQ=255)
                  ccGGATTTCAGCGTAAACTCGCCAAACTTAAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                                                                    <  2:1639040/150‑8 (MQ=255)
                   cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                                   <  1:1439309/150‑1 (MQ=255)
                    ggATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAg                                                                                                                                  >  1:1743818/1‑150 (MQ=255)
                    ggATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGATTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAg                                                                                                                                  <  1:811633/150‑1 (MQ=255)
                          aaGCGTAAAATCGCCAAACTTTAACACCTGCTTTCTAAGCGCAAATTCAATAAACTGTCGCTGATATGGTTTCATGCCTACGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                                            <  1:67791/149‑1 (MQ=255)
                                      gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTC                                                                                                                >  1:1018097‑M1/1‑150 (MQ=255)
                                       ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCT                                                                                                               >  2:643392‑M1/1‑150 (MQ=255)
                                        aaaacTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                                                                              <  2:953653‑M1/149‑1 (MQ=255)
                                          aaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                                            >  2:203564‑M1/1‑149 (MQ=255)
                                           aCTTTAACACCTGCTTGATAAGCGAAAATTCAATAAACTGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACCGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcg                                                                                                           <  1:606978‑M1/150‑3 (MQ=255)
                                            cTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgg                                                                                                          >  1:1439899‑M1/1‑147 (MQ=255)
                                            cTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCATTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                          >  1:96158‑M1/1‑147 (MQ=255)
                                              taaacaCCTGCTTGCTACGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcc                                                                                                        <  2:1526759‑M1/148‑6 (MQ=255)
                                               tAACACCTGCTTGCTAAGCGCAAAATCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAACAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                       <  1:112235‑M1/150‑7 (MQ=255)
                                                aaCACCTGCTTGCTAAGCGCAAATTCACTAAACTGGCGCGGATATGGTTTCATGCCTTAGCTCCTCATCTTACTATTCTACACACAAAAAAACGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgccccc                                                                                                      <  2:321158‑M1/150‑8 (MQ=255)
                                                aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccc                                                                                                      >  1:1968192‑M1/1‑143 (MQ=255)
                                                  cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                                                    >  2:204132‑M1/1‑141 (MQ=255)
                                                    ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                                  >  1:2003433‑M1/1‑139 (MQ=255)
                                                           gctACGCGAAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                           <  1:1205846‑M1/150‑19 (MQ=255)
                                                           gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                           <  1:719465‑M1/150‑19 (MQ=255)
                                                            ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                          >  1:1223264‑M1/1‑131 (MQ=255)
                                                               aGCGCAAAATCAATAAACTGGCGCTGAAATGGATTCATGCCTACCCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                       <  1:369706‑M1/150‑23 (MQ=255)
                                                                  gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                    >  2:361743‑M1/1‑125 (MQ=255)
                                                                  gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                    >  1:897251‑M1/1‑125 (MQ=255)
                                                                  gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                    >  1:1271325‑M1/1‑125 (MQ=255)
                                                                   cAAATTCAATAAACAGGCGCTGATATGGTCTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCcgtccccttcacttcaatgatgcgcc                                                                                   <  1:320527‑M1/150‑27 (MQ=255)
                                                                    aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgccc                                                                                  >  2:1600373‑M1/1‑123 (MQ=255)
                                                                          aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                            >  2:482214‑M1/1‑117 (MQ=255)
                                                                          aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctgc                                                                            >  2:1091145‑M1/1‑117 (MQ=255)
                                                                            tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgg                                                                          <  2:838783‑M1/150‑36 (MQ=255)
                                                                              aaCTGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATATAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                        <  2:1279577‑M1/150‑38 (MQ=255)
                                                                                 tGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCAACACACAAAAAAACGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                                     <  2:897414‑M1/150‑41 (MQ=255)
                                                                                  ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTAATTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                    <  2:1235483‑M1/150‑42 (MQ=255)
                                                                                  ggCGCTGATATGGTTGCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                    <  1:2001799‑M1/150‑42 (MQ=255)
                                                                                     gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                                                 >  2:896948‑M1/1‑106 (MQ=255)
                                                                                      cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGACAGCGCCACCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactac                                                                >  1:1317544‑M1/1‑105 (MQ=255)
                                                                                       tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                                               >  1:266454‑M1/1‑104 (MQ=255)
                                                                                        gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                              >  2:405971‑M1/1‑103 (MQ=255)
                                                                                        gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                              >  1:331319‑M1/1‑103 (MQ=255)
                                                                                        gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                              >  2:1220757‑M1/1‑103 (MQ=255)
                                                                                           atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttca                                                           >  2:1374764‑M1/1‑100 (MQ=255)
                                                                                             ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatt                                                         >  2:273863‑M1/1‑98 (MQ=255)
                                                                                             ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcggcccctgcactgcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                         >  2:1184856‑M1/1‑98 (MQ=255)
                                                                                             ggTTTCATGCCTTCGCTACTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCAGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                         <  2:252420‑M1/150‑53 (MQ=255)
                                                                                             ggTTTCATGACTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAACGCCGCCATCTGCagtcccctgcagttcaatgatgcgcccgtcatctgtcatctcgacgttcatt                                                         >  2:1251665‑M1/1‑98 (MQ=255)
                                                                                              gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                                        <  1:906620‑M1/150‑54 (MQ=255)
                                                                                              gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtgcatgg                                                        >  1:1604291‑M1/1‑97 (MQ=255)
                                                                                               tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtc                                                       <  2:226515‑M1/150‑55 (MQ=255)
                                                                                                 tCATGCGTTCGCTCTTCATCTTACTTTTCTACATACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                                     <  1:985141‑M1/150‑57 (MQ=255)
                                                                                                   aTGCTTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                                   <  1:724888‑M1/150‑59 (MQ=255)
                                                                                                    tGCCTTCGCTCCTCATCTTACTTTTCTACAGAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                  <  2:564794‑M1/150‑60 (MQ=255)
                                                                                                     gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                                 >  1:1201957‑M1/1‑90 (MQ=255)
                                                                                                         tCGCGCCTCATCTTAGATTTCTACATACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcg                                             <  1:380282‑M1/150‑65 (MQ=255)
                                                                                                          cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcgg                                            >  1:192752‑M1/1‑85 (MQ=255)
                                                                                                          cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCAGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcgg                                            >  1:955908‑M1/1‑85 (MQ=255)
                                                                                                            ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggca                                          >  1:1254544‑M1/1‑83 (MQ=255)
                                                                                                                tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagt                                      <  1:265355‑M1/150‑72 (MQ=255)
                                                                                                                 cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                                     >  1:1133766‑M1/1‑78 (MQ=255)
                                                                                                                  aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                    >  1:789585‑M1/1‑77 (MQ=255)
                                                                                                                  aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                    >  1:695880‑M1/1‑77 (MQ=255)
                                                                                                                  aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                    >  1:1881963‑M1/1‑77 (MQ=255)
                                                                                                                        cTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAGAAATCAGTTTGCAAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccctcgtcgatcatcactacgttcatgtcggtcactgcggcagagacgtcagcg                              >  1:447041‑M1/1‑71 (MQ=255)
                                                                                                                         ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgt                             <  2:1664681‑M1/150‑81 (MQ=255)
                                                                                                                             cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattc                         >  2:1495490‑M1/1‑66 (MQ=255)
                                                                                                                                cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccag                      >  2:411769‑M1/1‑63 (MQ=255)
                                                                                                                                cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCagtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccag                      >  2:605094‑M1/1‑63 (MQ=255)
                                                                                                                                cAGACAAAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccag                      <  1:80489‑M1/150‑88 (MQ=255)
                                                                                                                                cAGAAAAAACAAGGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGGCGCCTTCTGCcggcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccag                      <  1:1698711‑M1/150‑88 (MQ=255)
                                                                                                                                 aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacggattccagg                     >  1:703996‑M1/1‑62 (MQ=255)
                                                                                                                                  gACAAAAAAAAGGCGACTCATCATTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactaccttcatgtcggtcactgcggcacagtcttcaacgtattccagaa                    >  2:1660625‑M1/1‑61 (MQ=255)
                                                                                                                                  gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaaagatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttccacgtattccagat                    >  2:1361535‑M1/1‑61 (MQ=255)
                                                                                                                                   aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgctgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc                   <  2:843181‑M1/150‑91 (MQ=255)
                                                                                                                                   aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatc                   <  1:1374764‑M1/150‑91 (MQ=255)
                                                                                                                                    cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcagcgtattccagatcg                  >  2:749308‑M1/1‑59 (MQ=255)
                                                                                                                                    cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcg                  >  2:1073540‑M1/1‑59 (MQ=255)
                                                                                                                                    cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgggctctgcggcagagtcggcaacgtattcccgctcg                  >  2:827284‑M1/1‑59 (MQ=255)
                                                                                                                                     atttttaaggcgacTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgc                 >  1:1303712‑M1/7‑58 (MQ=255)
                                                                                                                                          ataggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgttctctgcggcagagtcttcaacgtattccagatcgcaaacc            <  1:1183170‑M1/148‑98 (MQ=255)
                                                                                                                                          aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaccgtattccagatcgcaaacc            >  1:644800‑M1/1‑53 (MQ=255)
                                                                                                                                            aGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacaacgttcctgtcggtctctgcggcagaggcgtcaacgtattcccgatcgcaaaccgc          >  1:376867‑M1/1‑51 (MQ=255)
                                                                                                                                            aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgc          <  1:932782‑M1/150‑100 (MQ=255)
                                                                                                                                             ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgttgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgct         <  2:1743818‑M1/150‑101 (MQ=255)
                                                                                                                                             ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgct         >  2:1031022‑M1/1‑50 (MQ=255)
                                                                                                                                               cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcggcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttc       <  2:2000062‑M1/150‑103 (MQ=255)
                                                                                                                                                gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcg      >  1:1913941‑M1/1‑47 (MQ=255)
                                                                                                                                                    catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatccctacgttcatgtcggtctctgcggcagagtcttcaacgtattccagattgcaaaccgcttcgccgg  >  2:1288502‑M1/1‑43 (MQ=255)
                                                                                                                                                    catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagcgtcttcaacgtattccagatcgcaaaccgcttcgccgt  >  2:321625‑M1/1‑43 (MQ=255)
                                                                                                                                                    catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgt  >  2:1565630‑M1/1‑43 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGC  >  NC_000913/3815668‑3815965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: