Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,222 Δ1 bp 6.0% intergenic (+30/‑111) udp → / → rmuC uridine phosphorylase/putative recombination limiting protein RmuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,017,2220T. 6.0% 1081.5 / 54.3 331intergenic (+30/‑111)udp/rmuCuridine phosphorylase/putative recombination limiting protein RmuC
Reads supporting (aligned to +/- strand):  ref base T (161/149);  new base . (11/9);  total (172/159)
Fisher's exact test for biased strand distribution p-value = 8.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 335 total aligned reads displayed.

TGGTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTG  >  NC_000913/4017074‑4017371
                                                                                                                                                     |                                                                                                                                                     
 ggTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTc                                                                                                                                                     <  1:770256/150‑1 (MQ=255)
  gTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                                                                                    >  1:297833/1‑150 (MQ=255)
  gTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                                                                                    <  2:578608/150‑1 (MQ=255)
   gaGCGGGCGTTCTCGTCACCCCCAGCCCGCCAGAGATCCCGAAGCCTGAGACCATCACACAAACCGAAAGCCATGCGGTGACAATCGGGGTGCAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCgg                                                                                                                                                   <  2:365429/149‑1 (MQ=255)
    aGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGc                                                                                                                                                  >  1:1741748/1‑150 (MQ=255)
          tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                                                                            >  2:1936734/1‑147 (MQ=255)
          tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                                                                            <  1:321106/150‑4 (MQ=255)
            ttATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcagac                                                                                                                                          <  1:831265/150‑6 (MQ=255)
              aTCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGCCGCGTTCGGCCTTTTGTAtt                                                                                                                                        >  1:1709505/1‑150 (MQ=255)
              aTCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAtt                                                                                                                                        >  2:1038609/1‑150 (MQ=255)
                cGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAtttt                                                                                                                                      <  2:1623540/150‑1 (MQ=255)
                 gTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAttttt                                                                                                                                     >  2:1340635/1‑150 (MQ=255)
                 gTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTAGGCCTTTTGTAttttt                                                                                                                                     >  1:1124626/1‑150 (MQ=255)
                  ttAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTg                                                                                                                                    <  1:851190/150‑1 (MQ=255)
                    aaCCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCg                                                                                                                                  >  2:61106/1‑150 (MQ=255)
                    aaCCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCg                                                                                                                                  <  1:1321902/150‑1 (MQ=255)
                     aCCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGt                                                                                                                                 >  1:747930/1‑150 (MQ=255)
                      ccGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTg                                                                                                                                    >  1:1654258/1‑146 (MQ=255)
                        gCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAgc                                                                                                                              <  2:1396549/150‑1 (MQ=255)
                         cACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAgcg                                                                                                                             <  1:95239/150‑1 (MQ=255)
                           cccAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGcc                                                                                                                           >  1:996750/1‑150 (MQ=255)
                                 aaGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                                                                     >  1:1051239/1‑150 (MQ=255)
                                    agaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaa                                                                                                                  >  2:1426914/1‑150 (MQ=255)
                                     gaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAgcg                                                                                                                             <  2:1233461/138‑1 (MQ=255)
                                     gaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCATGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                                                                                 >  2:1150593/1‑150 (MQ=255)
                                      aTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAAt                                                                                                                <  1:1035029/150‑1 (MQ=255)
                                       tCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATg                                                                                                               <  2:781211/150‑1 (MQ=255)
                                       tCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATg                                                                                                               >  2:640950/1‑150 (MQ=255)
                                             aTGCTGATACGATGAAAAAAACCGAAAGCCATGCGGTGAAAATCGTTGTGGAAGCTGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCTGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGTAAATGCCTTTc                                                                                                         <  1:1762239/150‑1 (MQ=255)
                                             aTGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                                                                         >  1:596575/1‑150 (MQ=255)
                                             aTGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                                                                         <  1:1926568/150‑1 (MQ=255)
                                             aTGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTc                                                                                                         <  2:1590583/150‑1 (MQ=255)
                                              tGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTcc                                                                                                        <  1:640431/150‑1 (MQ=255)
                                              tGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTcc                                                                                                        <  1:894903/150‑1 (MQ=255)
                                              tGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCAACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTcc                                                                                                        <  1:1619495/150‑1 (MQ=255)
                                                cTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCaa                                                                                                      >  2:901975/1‑150 (MQ=255)
                                                     aCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGa                                                                                                 <  1:340718/150‑1 (MQ=255)
                                                           aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                                                                           <  2:1106324/150‑1 (MQ=255)
                                                              caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                        <  2:669565/150‑1 (MQ=255)
                                                              caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                        >  1:1159636/1‑150 (MQ=255)
                                                              caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                        >  2:825333/1‑150 (MQ=255)
                                                                 acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                                                                                     <  1:1680724/150‑1 (MQ=255)
                                                                   cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcac                                                                                   <  2:1672105/150‑1 (MQ=255)
                                                                   cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcac                                                                                   <  1:208802/150‑1 (MQ=255)
                                                                    gAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcaca                                                                                  >  1:452573/1‑150 (MQ=255)
                                                                      aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAAt                                                                                >  2:1304981/1‑150 (MQ=255)
                                                                        gCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTc                                                                              >  2:1567240/1‑150 (MQ=255)
                                                                        gCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTc                                                                              <  1:1073054/150‑1 (MQ=255)
                                                                          cATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTa                                                                            >  2:1882768/1‑150 (MQ=255)
                                                                             gCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTg                                                                                                                                    <  2:863850/91‑1 (MQ=255)
                                                                              cGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtttt                                                                        <  1:1898208/150‑1 (MQ=255)
                                                                                 tGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtgtg                                                                     <  2:747930/150‑1 (MQ=255)
                                                                                   aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGc                                                                                                                       >  2:1872065/1‑98 (MQ=255)
                                                                                   aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg                                                                   <  1:538334/150‑1 (MQ=255)
                                                                                    aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                                                                  <  1:1340635/150‑1 (MQ=255)
                                                                                    aaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGAAAATCTCTTGTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATATTTGCGTAGCCCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTGTATTTTGTGCgg                                                                  >  2:772802/1‑150 (MQ=255)
                                                                                      aTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGGTCGGCCGTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGt                                                                <  2:18459/150‑1 (MQ=255)
                                                                                       tCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa                                                               <  2:221682/150‑1 (MQ=255)
                                                                                         gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCCCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAg                                                             >  1:1777722/1‑150 (MQ=255)
                                                                                            gtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg                                                          >  2:417462/1‑150 (MQ=255)
                                                                                            gtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg                                                          <  1:1421773/150‑1 (MQ=255)
                                                                                            gtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttg                                                          >  1:969733/1‑150 (MQ=255)
                                                                                                aaGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGc                                                      <  2:321178/150‑1 (MQ=255)
                                                                                                 aGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCg                                                     <  1:67584/150‑1 (MQ=255)
                                                                                                     gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCag                                                 >  1:284690/1‑150 (MQ=255)
                                                                                                       gcgTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGCCGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCGCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACCGCACAATTCGATTTTGGGCGGGTAAGTTGTTGCGTCagga                                               >  2:458293/1‑150 (MQ=255)
                                                                                                           cgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCg                                           <  2:119070/150‑1 (MQ=255)
                                                                                                              ctgctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGCAAGTTGTTGCGTCAGGAGGCGTtg                                        <  1:61378/150‑1 (MQ=255)
                                                                                                                gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgtg                                      >  1:1029007/1‑150 (MQ=255)
                                                                                                                gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgtg                                      <  2:596745/150‑1 (MQ=255)
                                                                                                                 ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTgg                                     >  2:949936/1‑150 (MQ=255)
                                                                                                                  tgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCagga                                               >  2:1656737/1‑139 (MQ=255)
                                                                                                                  tgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGa                                    >  1:1256442/1‑150 (MQ=255)
                                                                                                                   gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGAt                                   <  1:146080/150‑1 (MQ=255)
                                                                                                                    tAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGAtt                                  >  1:1519763/1‑150 (MQ=255)
                                                                                                                     aaTTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGAttt                                 >  2:1899010/1‑150 (MQ=255)
                                                                                                                     aaTTCTCTTCTCCTCTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGAttt                                 >  2:654380/1‑150 (MQ=255)
                                                                                                                      aTTCTCTTCTCCTTTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAAAGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGTGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTtc                                <  1:1151593/150‑1 (MQ=255)
                                                                                                                      aTTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTtc                                <  2:969733/150‑1 (MQ=255)
                                                                                                                      aTTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTtc                                <  2:1741748/150‑1 (MQ=255)
                                                                                                                         ctcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCa                             >  2:1079397/1‑150 (MQ=255)
                                                                                                                           cttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAAt                           >  1:740828/1‑150 (MQ=255)
                                                                                                                             tctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCa                         <  2:195080/150‑1 (MQ=255)
                                                                                                                             tctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCa                         <  1:1100307/150‑1 (MQ=255)
                                                                                                                             tctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCa                         <  1:804574/150‑1 (MQ=255)
                                                                                                                              ctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCAt                        >  1:1302006/1‑150 (MQ=255)
                                                                                                                              ctcCTGTCTGAAGGCCGACGCGTTCGGCCTGTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCGTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCAt                        >  1:1482015/1‑150 (MQ=255)
                                                                                                                                ccTGTCTGAAAGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                                                                                     >  2:1464343/1‑87 (MQ=255)
                                                                                                                                       aaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGc               >  1:1227846/2‑150 (MQ=255)
                                                                                                                                       aaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGc               >  1:1205853/2‑150 (MQ=255)
                                                                                                                                        aaGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCa              <  1:1567240/150‑1 (MQ=255)
                                                                                                                                         aGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAg             <  1:1871859/150‑1 (MQ=255)
                                                                                                                                         aGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAg             <  1:181031/150‑1 (MQ=255)
                                                                                                                                          ggCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGt            >  2:1691021/1‑150 (MQ=255)
                                                                                                                                           gCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGGCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGtt           <  1:391261/150‑1 (MQ=255)
                                                                                                                                            ccGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTa          >  2:258032/1‑150 (MQ=255)
                                                                                                                                            ccGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTa          >  1:1953415/1‑150 (MQ=255)
                                                                                                                                              gACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTAtt        >  2:144366/1‑150 (MQ=255)
                                                                                                                                                aCGCGTTCGGCGTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTTCTGTTAAGTTGTTGCGTCATGAGGCTTTGTTGATTTCTCAATCATGGTTTACGCAGTTATTg       <  1:1376947/150‑1 (MQ=255)
                                                                                                                                                   cgTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGt    >  1:721779/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTTG  >  NC_000913/4017074‑4017371

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