Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,809 Δ1 bp 38.3% intergenic (‑41/+25) 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,8090G.38.3% 769.9 / inf 586intergenic (‑41/+25)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (130/150);  new base . (193/31);  total (327/259)
Fisher's exact test for biased strand distribution p-value = 3.30e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 629 total aligned reads displayed.

CGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAAT  >  NC_000913/3815660‑3815958
                                                                                                                                                       |                                                                                                                                                     
cGTTGAAGAATTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGAAAATTC‑AATAAACTGGCGCTGATTTGTTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAg                                                                                                                                                       <  2:1329957/150‑1 (MQ=255)
  ttGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAgcg                                                                                                                                                     >  2:1489838/1‑148 (MQ=255)
    gaagaaATAGGGGCTTTTGCGCCCGGTTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAgcgac                                                                                                                                                   <  2:6197847/150‑5 (MQ=255)
       gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATCTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACtc                                                                                                                                                >  2:5323576/1‑150 (MQ=255)
       gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACtc                                                                                                                                                >  2:4905473/1‑150 (MQ=255)
        aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACtcat                                                                                                                                              >  2:564887/1‑150 (MQ=255)
               ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGc                                                                                                                                       >  1:559502/1‑150 (MQ=255)
                  ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCtt                                                                                                                                    >  1:857453/1‑150 (MQ=255)
                     tGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACGGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTaaa                                                                                                                                 >  2:8996305/1‑150 (MQ=255)
                          ccgtttttCAGCTAAAACCCGCAAAAGTTTAACCCGTGCTTGCAAAGCCAAATTTA‑AATAAACTGGCTCTTATTTGGTTTCTTGCCTTCGCTCCTCATCTTCATTTTCTACACA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                            <  2:1124726/145‑1 (MQ=255)
                          ccGGATTTCAGCGTAAACTCGCCAAACTTTAACCCCGGTTTGCTAAGCGCAAATTC‑AATACACGGGCGGTGATCTGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACACA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                            <  2:10862949/150‑1 (MQ=255)
                          ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc                                                                                                                             >  2:5837305/1‑150 (MQ=255)
                           cgtaTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                            <  1:8011677/147‑1 (MQ=255)
                           cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                            >  1:3802407/1‑150 (MQ=255)
                           cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAg                                                                                                                           >  2:179587/1‑150 (MQ=255)
                              ttttCCGCGTAAACGCCCAAAACTTACCCCCGTGGTTGAAAAGCG‑AAATTCAAATAAACGGGGGTTTTTTTGTTTTCAGGCGTTCGCCCCTCATCTTA‑TTTTC‑ACACACAAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                                                                        <  2:5046724/149‑1 (MQ=255)
                               tttCAGCTTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATTTGTTTACATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                                                                       <  1:355707/150‑1 (MQ=255)
                               tttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                                                                        >  1:6787969/1‑150 (MQ=255)
                                ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                                                                                                       >  1:463978/1‑150 (MQ=255)
                                ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                                                                      <  1:6546722/150‑1 (MQ=255)
                                    gCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgc                                                                                                                  >  2:10058414/1‑150 (MQ=255)
                                     cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                                                                 >  1:10045821/1‑150 (MQ=255)
                                       tAAACTCGCAAAACTTTAACACCTGCTTGATAAGCTAAAATTA‑ATTAAACTGGCGCTTATATGTTTTATTGCCTTCGCTCCTCATCTTACTTTTCTACATA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCC                                                                                                               <  2:7959331‑M1/150‑1 (MQ=255)
                                        aaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCG                                                                                                              >  2:10759808‑M1/1‑150 (MQ=255)
                                            tCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTT                                                                                                          >  2:8396166‑M1/1‑150 (MQ=255)
                                                cAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTG                                                                                                       >  1:7850385‑M1/1‑150 (MQ=255)
                                                  aaCTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGt                                                                                                    >  2:11087694/1‑150 (MQ=255)
                                                  aaCTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGt                                                                                                    >  1:9803797/1‑150 (MQ=255)
                                                    cTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg                                                                                                  >  2:8083532/1‑150 (MQ=255)
                                                     tttAACACCTGCTTGCTAAGCGAAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                                  >  2:7420846‑M1/1‑146 (MQ=255)
                                                      ttAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGc                                                                                                 >  1:3105579/1‑150 (MQ=255)
                                                      ttAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                >  2:5497188‑M1/1‑144 (MQ=255)
                                                       tAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtccc                                                                                                <  2:1640374‑M1/150‑7 (MQ=255)
                                                          cacCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                            >  1:6041701/1‑150 (MQ=255)
                                                           acCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACact                                                                                           >  2:10014726/1‑148 (MQ=255)
                                                           acCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACCAt                                                                                           >  1:8188348/1‑150 (MQ=255)
                                                           acCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAAACAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                                                                                           >  1:10263511/1‑150 (MQ=255)
                                                            ccTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgca                                                                                           <  1:2321681‑M1/150‑12 (MQ=255)
                                                              tgcttgctACGCTCAATTAC‑AAAAACGGGGGGCTTATTTGTTTCCTGGCCTGCGCCCCTCTACTAACTTTACAACACA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGga                                                                                        <  1:1733899/150‑1 (MQ=255)
                                                              tgcttgctAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGga                                                                                        <  1:2958800/150‑1 (MQ=255)
                                                                     tAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                                                                                 >  2:4105904/1‑150 (MQ=255)
                                                                          gcAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                             <  2:724755‑M1/150‑26 (MQ=255)
                                                                          gcAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCt                                                                            >  1:779583/1‑150 (MQ=255)
                                                                          gcAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCt                                                                            >  1:1448013/1‑150 (MQ=255)
                                                                             aaTTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                                         <  2:4287906/150‑1 (MQ=255)
                                                                                   aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTGTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGGCTTCTGCGTCGCTACAATGGATTCGATTCGCCTCGGGCCag                                                                    >  1:5834490/1‑150 (MQ=255)
                                                                                   aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCag                                                                    >  1:4974245/1‑150 (MQ=255)
                                                                                    aTAACCTGGCGGTTTTTTGGTTTCTTGCCTTCGCTCCTCACCTTAATTTTATCCAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCaga                                                                   <  2:2586563/150‑1 (MQ=255)
                                                                                       aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtc                                                                 <  1:9042832‑M1/150‑38 (MQ=255)
                                                                                             cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTACTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                                                           <  1:1595282‑M1/150‑44 (MQ=255)
                                                                                              gcTTATTTGTTTTCTTGCCTTCGCTCCTCATCTTATTTTAATACACA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGa                                                         <  2:359750/150‑1 (MQ=255)
                                                                                              gcTGATATGGTTTCTTGCCTTCGATCCACATCTTACTTTACTACACA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGa                                                         <  2:656778/150‑1 (MQ=255)
                                                                                               gtGATTTGTTTCCTTGCGTGCGCTCCTCATCTTTATTTTCAACAGA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt                                                        <  2:5191521/149‑1 (MQ=255)
                                                                                                 tttttGGTTACATGCCTTCGCTCATCATCTTACTTTACTACACA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                                      <  1:1077460/146‑1 (MQ=255)
                                                                                                   tttGGTTTCATGCCTGGGCTCCTCATCTTATTTTTCAACACA‑‑AAAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                                    <  2:21962/148‑1 (MQ=255)
                                                                                                   tatGGTTTCATGCCTTCTCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                     <  1:4333191‑M1/150‑50 (MQ=255)
                                                                                                   tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttc                                                     <  1:5691949‑M1/150‑50 (MQ=255)
                                                                                                   taggttttCTTGCTTTCGATCCTCATCTGAATTTTCTACATA‑‑AAAAAAAAAGCTACTCATCAGTCGACTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                                    <  2:5999011/145‑1 (MQ=255)
                                                                                                    atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCC‑ACAAGATGAGTa                                                   <  1:6382482/150‑1 (MQ=255)
                                                                                                     tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTCTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcat                                                   >  2:10584526‑M1/1‑99 (MQ=255)
                                                                                                     tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                  >  1:8719544/1‑150 (MQ=255)
                                                                                                      ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatg                                                  <  1:7399642‑M1/150‑53 (MQ=255)
                                                                                                      ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGc                                                 >  2:1895308/1‑150 (MQ=255)
                                                                                                       gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                                 <  2:3075301‑M1/150‑54 (MQ=255)
                                                                                                        tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGctc                                               >  2:6929502/1‑150 (MQ=255)
                                                                                                          tctggCCTGCGCGCCTCATTTTATTTTTCAACAGAC‑‑AAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                             <  2:28060/146‑1 (MQ=255)
                                                                                                          tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCCACACGATGAGTAGcgcgg                                             >  2:764085/1‑146 (MQ=255)
                                                                                                            tggctttCGCTCCTCCTGTTACTTTTCTACACA‑‑AAAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                                           <  1:409716/145‑1 (MQ=255)
                                                                                                            aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtc                                            >  2:7100358‑M1/1‑92 (MQ=255)
                                                                                                             tGCTTTCGCGCCTCTTATTACTTTACTACCCA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAt                                          <  1:5182003/150‑1 (MQ=255)
                                                                                                               ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                         >  2:9520894/1‑150 (MQ=255)
                                                                                                               ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTACTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctct                                         >  1:53331‑M1/1‑89 (MQ=255)
                                                                                                               ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                                        <  2:3154747/150‑1 (MQ=255)
                                                                                                                cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                                       >  2:8843290/1‑150 (MQ=255)
                                                                                                                 ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGt                                      >  1:9422648/1‑150 (MQ=255)
                                                                                                                 ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCACGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGt                                      >  2:9733927/1‑150 (MQ=255)
                                                                                                                    gCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                                   >  1:10310436/1‑150 (MQ=255)
                                                                                                                     ctcctcATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                                  >  1:1966667/1‑150 (MQ=255)
                                                                                                                      tcctcATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACg                                 >  2:24148/1‑150 (MQ=255)
                                                                                                                          cATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                              <  2:3436995‑M1/150‑73 (MQ=255)
                                                                                                                            tCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtctt                            >  1:8876403‑M1/1‑76 (MQ=255)
                                                                                                                            tCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGc                           <  1:9676794/150‑1 (MQ=255)
                                                                                                                              ttACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt                         >  2:4476167/1‑150 (MQ=255)
                                                                                                                              ttACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt                         >  1:11058239/1‑150 (MQ=255)
                                                                                                                                aCTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                       >  1:10480898/1‑150 (MQ=255)
                                                                                                                                 cTTTTCTCCACAC‑‑AAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                      <  2:10379954/150‑1 (MQ=255)
                                                                                                                                  ttttCTCCCGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                     <  2:6801783/150‑1 (MQ=255)
                                                                                                                                  ttttCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                     >  2:3804326/1‑150 (MQ=255)
                                                                                                                                   tttCAACACAC‑‑AAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATGCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                    <  1:3409079/150‑1 (MQ=255)
                                                                                                                                    ttCTACACA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                   <  1:262278/150‑1 (MQ=255)
                                                                                                                                       tACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattcc                 >  1:7044242‑M1/1‑65 (MQ=255)
                                                                                                                                       tACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                >  2:10253379/1‑150 (MQ=255)
                                                                                                                                           cagaaaaaaaaaAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc              <  2:9061205/146‑1 (MQ=255)
                                                                                                                                           cacacaaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc              <  1:9165652/145‑1 (MQ=255)
                                                                                                                                            ataaaaaaaaaAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc             <  2:8549687/147‑1 (MQ=255)
                                                                                                                                                aaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa         >  2:4617011/1‑150 (MQ=255)
                                                                                                                                                aaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa         >  2:9206112/1‑150 (MQ=255)
                                                                                                                                                aaaaaaaCGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa         <  2:9345443/150‑1 (MQ=255)
                                                                                                                                                      aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgc   >  1:1290803‑M1/1‑51 (MQ=255)
                                                                                                                                                       ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat  <  2:3964356/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTC‑AATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAAT  >  NC_000913/3815660‑3815958

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑