Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,819 C→G 5.6% 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,8190CG 5.6% 396.3 / 12.6 125intergenic (‑51/+15)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 = 2.54e-01

Only 100 of 137 total aligned reads displayed.

ATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCG  >  NC_000913/3815670‑3815968
                                                                                                                                                     |                                                                                                                                                     
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)
       cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                                                                                <  1:1284806/150‑1 (MQ=255)
       cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATAGGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGcc                                                                                                                                                >  1:926082/1‑150 (MQ=255)
          ttGCTCCCGGATTTCAGCGTAAACTGGCCAAACTTTAACACCTGCTTGGTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCTTGCCTTCGCTCCTCTTCTTAGTTTTCTTCAGACAAAAAAATGGCGACTGATGAGTCGCCttt                                                                                                                                             <  2:650498/150‑2 (MQ=255)
            gcgcCCGGTTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAACCTGGCGCTGATATGGTTTCATGCCTCCGCTCCTCATCTTACTTTTCAACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTaaa                                                                                                                                           <  1:1570727/150‑1 (MQ=255)
             cgcaaggATATAAGCGTAAACTCGCCTAACTTTAACACCTGCTTGCTAATCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAACGGCGACTCATCAGTCGCCTTaaaa                                                                                                                                          <  1:1209565/145‑1 (MQ=255)
             cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                                                                          >  2:1844206/1‑150 (MQ=255)
               cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGAAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTAACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAt                                                                                                                                        <  1:771204/150‑1 (MQ=255)
                ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                                                                       <  1:1773558/150‑8 (MQ=255)
                ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc                                                                                                                                       >  1:1103044/1‑150 (MQ=255)
                ccGGATTTCAGCGTAAACTCGCCAAACTTAAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                                                                       <  2:1639040/150‑8 (MQ=255)
                ccGGATTTCAGCGTAAACTCGCCAAACTTAAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTGATGAGTCGCCttttttttg                                                                                                                                       <  2:1639044/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)
                    aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGtt                                                                                                                                   <  2:36818/150‑1 (MQ=255)
                        cAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                                               <  1:701818/150‑1 (MQ=255)
                         aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCa                                                                                                                              >  2:1464152/1‑150 (MQ=255)
                                    gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTC                                                                                                                   >  1:1018097‑M1/1‑150 (MQ=255)
                                    gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTACCTTTTCAACAGACAAAAAAAAGGAGAGTCATCAGTCGCCTTAAAAAACCGTTACCCAGCCCCGCCTTa                                                                                                                   >  2:823524‑M1/1‑149 (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)
                                          cTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCATTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                             >  1:96158‑M1/1‑147 (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)
                                                          ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                             >  1:1223264‑M1/1‑131 (MQ=255)
                                                             aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                          >  2:1560627‑M1/1‑128 (MQ=255)
                                                             aGCGCAAAATCAATAAACTGGCGCTGAAATGGATTCATGCCTACCCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatg                                                                                          <  1:369706‑M1/150‑23 (MQ=255)
                                                                gcAAATTCAATAACCTGGCGCTGATATGGTTTCATGCCTTCGATCCTCATCTTACTCTTCTACCGACAACAAAAACGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCACCGCCGCCTTCAGCcgtcccctgcacttcaatcatgcgc                                                                                       <  2:847479‑M1/150‑26 (MQ=255)
                                                                gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                       >  1:897251‑M1/1‑125 (MQ=255)
                                                                 cAAATTCAATAAACAGGCGCTGATATGGTCTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTCCcgtccccttcacttcaatgatgcgcc                                                                                      <  1:320527‑M1/150‑27 (MQ=255)
                                                                       cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctt                                                                                >  2:562091‑M1/1‑118 (MQ=255)
                                                                        aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttc                                                                               >  2:482214‑M1/1‑117 (MQ=255)
                                                                        aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtctgc                                                                               >  2:1091145‑M1/1‑117 (MQ=255)
                                                                            aaCTGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATATAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                           <  2:1279577‑M1/150‑38 (MQ=255)
                                                                            aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTAATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                           <  2:1234427‑M1/150‑38 (MQ=255)
                                                                               tGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCAACACACAAAAAAACGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatc                                                                        <  2:897414‑M1/150‑41 (MQ=255)
                                                                                ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                       <  1:1403796‑M1/150‑42 (MQ=255)
                                                                                ggCGCTGATATGGTTGCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatca                                                                       <  1:2001799‑M1/150‑42 (MQ=255)
                                                                                   gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                                                    >  2:896948‑M1/1‑106 (MQ=255)
                                                                                     tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                                                  >  1:266454‑M1/1‑104 (MQ=255)
                                                                                      gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                                 >  1:331319‑M1/1‑103 (MQ=255)
                                                                                      gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                                 >  2:1220757‑M1/1‑103 (MQ=255)
                                                                                      gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgt                                                                 >  2:405971‑M1/1‑103 (MQ=255)
                                                                                           ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatt                                                            >  2:273863‑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)
                                                                                               tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgg                                                        <  1:1660318‑M1/150‑57 (MQ=255)
                                                                                                 aTGCTTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                                      <  1:724888‑M1/150‑59 (MQ=255)
                                                                                                  tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                     >  1:620847‑M1/1‑91 (MQ=255)
                                                                                                  tGCCTTCGCTCCTCATCTTACTTTTCTACAGAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtct                                                     <  2:564794‑M1/150‑60 (MQ=255)
                                                                                                   gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctc                                                    >  1:1201957‑M1/1‑90 (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)
                                                                                                                aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                       >  1:789585‑M1/1‑77 (MQ=255)
                                                                                                                aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                       >  1:695880‑M1/1‑77 (MQ=255)
                                                                                                                aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtct                                       >  1:1881963‑M1/1‑77 (MQ=255)
                                                                                                                  cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttc                                     >  1:1170154‑M1/1‑75 (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)
                                                                                                                               aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacggattccagg                        >  1:703996‑M1/1‑62 (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)
                                                                                                                                         aaGGCGACTGATGAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCagtcccctgcacttcaatgatgagcccgtcttcggtcatcactacgttcatgtcggtctatgcggcagagtcttcaacgtattacagatcgcaaaccg              >  1:432258‑M1/1‑52 (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)
                                                                                                                                                    tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc   <  2:1505450‑M1/150‑110 (MQ=255)
                                                                                                                                                    tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc   <  2:1271325‑M1/150‑110 (MQ=255)
                                                                                                                                                    tcaGTCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTGCTGCcgtcgcctgcacttcaatgatgcgcccgtcttcggtcatcactactttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc   <  1:362807‑M1/150‑110 (MQ=255)
                                                                                                                                                    tcaGGCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTGCTGCcgtccgctgcacttgaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttc   <  1:208582‑M1/150‑110 (MQ=255)
                                                                                                                                                     caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttca  <  2:1223264‑M1/150‑111 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCG  >  NC_000913/3815670‑3815968

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: