Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,587,886 C→A 7.6% intergenic (+23/+23) mrr → / ← yjiA methylated adenine and cytosine restriction protein/P‑loop guanosine triphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,587,8860CA 7.6% 276.8 / 9.0 85intergenic (+23/+23)mrr/yjiAmethylated adenine and cytosine restriction protein/P‑loop guanosine triphosphatase
Reads supporting (aligned to +/- strand):  ref base C (43/35);  new base A (2/5);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

CAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGC  >  NC_000913/4587738‑4588035
                                                                                                                                                    |                                                                                                                                                     
cAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGACTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTgcc                                                                                                                                                      >  2:1427068/1‑149 (MQ=255)
   tCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCt                                                                                                                                                   <  2:818374/150‑1 (MQ=255)
   tCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCt                                                                                                                                                   <  1:1391281/150‑1 (MQ=255)
    ccGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCtt                                                                                                                                                  >  2:465730/1‑150 (MQ=255)
     cGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTa                                                                                                                                                 >  2:557370/1‑150 (MQ=255)
     cGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTa                                                                                                                                                 >  2:1710509/1‑150 (MQ=255)
      gTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTAt                                                                                                                                                >  1:1004570/1‑150 (MQ=255)
          aGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCgg                                                                                                                                            >  1:389767/1‑150 (MQ=255)
           gggTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATGAGGCCGGATTCGGCGGCGCCTTATCGGGc                                                                                                                                           <  2:919975/150‑1 (MQ=255)
            ggTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAACCGAAGTAGGGGGTTCTTCACGTTTGTTGACGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATAAGGCGCAGCCGCATCCGGcc                                                                                                                                          <  1:1356199/150‑1 (MQ=255)
            ggTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGcc                                                                                                                                          <  1:131576/150‑1 (MQ=255)
            ggTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGcc                                                                                                                                          <  2:405961/150‑1 (MQ=255)
            ggTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATAAGGCGCAGCCGCATCCGGcc                                                                                                                                          <  1:213636/150‑1 (MQ=255)
             gTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATAAGGCGCAGCCGCATCCGGcct                                                                                                                                         <  1:509212/150‑2 (MQ=255)
              tATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATAAGGCGCAGCCGCATCCGGcctg                                                                                                                                        <  1:408082/150‑3 (MQ=255)
               aTGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCCCGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTTCGCCTTATCCGGCCCAt                                                                                                                                       >  2:1408852/1‑150 (MQ=255)
               aTGGTGTAGGGTGATGGGGAACGCCTGTTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCAt                                                                                                                                       <  2:1796230/150‑1 (MQ=255)
                tGGTGTTGGTTGATGGGGACCGCCTGGTGCACTTATTGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTTCCGAAACTGGATATGGAGTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATACGGCCCATa                                                                                                                                      >  2:1333407/1‑150 (MQ=255)
                 ggTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATaa                                                                                                                                     >  2:1978375/1‑150 (MQ=255)
                    gttggttgATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAAccc                                                                                                                                  >  2:716554/1‑150 (MQ=255)
                         ttgATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTAc                                                                                                                             <  2:1340555/150‑1 (MQ=255)
                            aTGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTc                                                                                                                          >  2:10739/1‑150 (MQ=255)
                            aTGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATAAGGCGCAGCCGCATCCGGcctg                                                                                                                                        <  2:1930868/136‑3 (MQ=255)
                            aTGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATAAGGCGCAGCCGCATCCGGcctg                                                                                                                                        >  1:1930868/1‑134 (MQ=255)
                                      ccTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCg                                                                                                                >  2:816920/1‑150 (MQ=255)
                                         ggTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCaa                                                                                                             <  1:1950550/150‑1 (MQ=255)
                                           tGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCCCGTTTGTTGAAGGTGCCGAAACTGGATATGGACTAATTTGAGTGAAATATCATGCCGGATGCGACTGAGCCTTATCAGTCCCAGAACCCCTTACTTCCTCCACACCGCAAAc                                                                                                           >  2:343963/1‑150 (MQ=255)
                                                ttAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGc                                                                                                      <  2:469678/150‑1 (MQ=255)
                                                  taTGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAAc                                                                                                                    <  1:794773/133‑1 (MQ=255)
                                                  aaTGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAAc                                                                                                                    >  2:794773/1‑134 (MQ=255)
                                                    tGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGa                                                                                                  <  1:437294/150‑1 (MQ=255)
                                                     gATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGaa                                                                                                 >  1:294605/1‑150 (MQ=255)
                                                       tCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAAtc                                                                                               <  2:389767/150‑1 (MQ=255)
                                                          aaaaCGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATctct                                                                                            <  1:1410867/150‑1 (MQ=255)
                                                          aaaaCGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATctct                                                                                            >  1:802273/1‑150 (MQ=255)
                                                           aaaCGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCtt                                                                                           >  2:1347968/1‑150 (MQ=255)
                                                           aaaCGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCtt                                                                                           <  1:686632/150‑1 (MQ=255)
                                                           aaaCGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACACAGCCCGAATCTCtt                                                                                           >  1:206801/1‑150 (MQ=255)
                                                            aaCGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCCCAAACGCACCCCGAATCTCAtc                                                                                          >  2:1655109/1‑150 (MQ=255)
                                                                   tAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGc                                                                                   <  1:1710509/150‑1 (MQ=255)
                                                                    aGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCa                                                                                  <  1:557370/150‑1 (MQ=255)
                                                                     ggggTTTCTTCACGTTTGTTAAGGTGCCGAAACTGGATATGGACTATTTTGAGTCAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCACATAACCCCTTACTTCCTCAACCCCGCAAACGCAGACCGAATCTCTTCCTCCGGCAGc                                                                                >  2:600148/1‑150 (MQ=255)
                                                                       ggTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCt                                                                               <  1:1493040/150‑1 (MQ=255)
                                                                           tcttcACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGAt                                                                           <  1:1629301/150‑1 (MQ=255)
                                                                           tcttcACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGAt                                                                           >  2:1346825/1‑150 (MQ=255)
                                                                            cttcACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATc                                                                          <  2:743596/150‑1 (MQ=255)
                                                                             ttcACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATcc                                                                         <  1:943036/150‑1 (MQ=255)
                                                                              tcACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATccc                                                                        >  1:712387/1‑150 (MQ=255)
                                                                              tcACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATccc                                                                        <  1:10739/150‑1 (MQ=255)
                                                                               cACGTTTGTTGAAGGAGCAGAAAATGGATATGGACTATTTTGAGTGAAATATCACGCCGGATGCCGCTGCGCCTTATCCGGCCCATAACCCCATACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGCCACCTGGATccac                                                                       >  2:642426/1‑148 (MQ=255)
                                                                                     tgttgAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAac                                                                 <  1:1845278/150‑1 (MQ=255)
                                                                                      gttgAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAaca                                                                >  2:1898194/1‑150 (MQ=255)
                                                                                      gttgAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAaca                                                                >  2:658182/1‑150 (MQ=255)
                                                                                      gttgAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCAGAATCTCTTCCTACGGCAGCTGGATCCCGATAAaca                                                                >  1:522542/1‑150 (MQ=255)
                                                                                        tgAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACAcc                                                              >  1:1405187/1‑150 (MQ=255)
                                                                                             gTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGt                                                         >  2:620579/1‑150 (MQ=255)
                                                                                               gCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCTGCAGCTGGATCCCGATAAACACAATCGTga                                                       >  2:438998/1‑149 (MQ=255)
                                                                                                 cGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTa                                                     >  2:1445906/1‑150 (MQ=255)
                                                                                                  gAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTAt                                                    <  1:680607/150‑1 (MQ=255)
                                                                                                  gAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTAt                                                    >  1:139765/1‑150 (MQ=255)
                                                                                                       tGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGt                                               <  1:438998/150‑1 (MQ=255)
                                                                                                         gATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACCCCATCGTGCTAGGCGGttt                                             >  1:852378/1‑150 (MQ=255)
                                                                                                          atatGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGtttt                                            <  1:1978375/150‑1 (MQ=255)
                                                                                                              ggACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGc                                                                                   >  1:871056/1‑107 (MQ=255)
                                                                                                              ggACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGc                                                                                   <  2:871056/107‑1 (MQ=255)
                                                                                                              ggACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCAt                                        <  1:1445906/150‑1 (MQ=255)
                                                                                                              ggACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCAt                                        >  2:1668969/1‑150 (MQ=255)
                                                                                                               gACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATc                                       <  1:502976/150‑1 (MQ=255)
                                                                                                                 cTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGc                                     >  2:1231078/1‑150 (MQ=255)
                                                                                                                 cTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGc                                     >  2:465907/1‑150 (MQ=255)
                                                                                                                   aTTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGccc                                   <  2:139765/150‑1 (MQ=255)
                                                                                                                       tGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACg                               >  2:802027/1‑150 (MQ=255)
                                                                                                                         aGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGc                             <  1:1167018/150‑1 (MQ=255)
                                                                                                                          gTGAAATATCCGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCACCCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGGTTTTCATCGCCCCCCGGcc                            >  2:1979069/1‑150 (MQ=255)
                                                                                                                          gTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGcc                            >  1:1008246/1‑150 (MQ=255)
                                                                                                                          gTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGcc                            <  1:1820184/150‑1 (MQ=255)
                                                                                                                          gTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGcc                            <  1:1898194/150‑1 (MQ=255)
                                                                                                                            gAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTg                          >  2:1763643/1‑150 (MQ=255)
                                                                                                                            gAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTg                          >  1:879477/1‑150 (MQ=255)
                                                                                                                              aaTATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCACCCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTc                        >  2:1806539/1‑150 (MQ=255)
                                                                                                                                    aGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTc                  <  1:658182/150‑1 (MQ=255)
                                                                                                                                     ggCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCg                 >  2:1598806/1‑150 (MQ=255)
                                                                                                                                        cGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGgcg              <  1:343963/150‑1 (MQ=255)
                                                                                                                                          gATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCt            <  1:854283/150‑1 (MQ=255)
                                                                                                                                               ggCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTaga       >  2:1969364/1‑150 (MQ=255)
                                                                                                                                                 cTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAgg     >  2:720489/1‑150 (MQ=255)
                                                                                                                                                    cgcCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGgcgc  >  2:454151/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTTGAAGGTGCCGAAACTGGATATGGACTATTTTGAGTGAAATATCAGGCCGGATGCGGCTGCGCCTTATCCGGCCCATAACCCCTTACTTCCTCAACCCCGCAAACGCAGCCCGAATCTCTTCCTCCGGCAGCTGGATCCCGATAAACACCATCGTGCTATGCGGTTTTTCATCGCCCCACGGCCTGTCCCAGTCGGCGCTGTAGAGGCGC  >  NC_000913/4587738‑4588035

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