Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,726,045 C→A 5.6% intergenic (‑299/+24) kgtP ← / ← rrfG alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,726,0450CA 5.6% 318.3 / 2.7 92intergenic (‑299/+24)kgtP/rrfGalpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base C (43/42);  new base A (2/3);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAA  >  NC_000913/2725898‑2726190
                                                                                                                                                   |                                                                                                                                                   
gATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTc                                                                                                                                                     <  1:865682/148‑1 (MQ=255)
gATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCa                                                                                                                                                    >  4:23167/1‑149 (MQ=255)
  tAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAgg                                                                                                                                                 <  3:931540/149‑1 (MQ=255)
    aaTTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGa                                                                                                                                                <  2:126144/148‑1 (MQ=255)
     aTTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGAt                                                                                                                                               <  1:844896/148‑1 (MQ=255)
      tttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATgg                                                                                                                                             >  4:868928/1‑149 (MQ=255)
       ttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATgg                                                                                                                                             <  1:393698/148‑1 (MQ=255)
       ttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATgg                                                                                                                                             >  2:550223/1‑148 (MQ=255)
        tCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGcc                                                                                                                                           <  3:516252/149‑1 (MQ=255)
         cTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGACGGATGGcc                                                                                                                                           <  2:14998/148‑1 (MQ=255)
          tGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGACGGATGGCCtt                                                                                                                                         <  3:71729/149‑1 (MQ=255)
          tGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCtt                                                                                                                                         <  3:724478/149‑1 (MQ=255)
           ggTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGACGGATGGCCttt                                                                                                                                        <  4:348647/149‑1 (MQ=38)
            gTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTc                                                                                                                                       >  4:399622/1‑149 (MQ=255)
             tAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGACGGATGGCCtttt                                                                                                                                       <  1:809497/148‑2 (MQ=255)
             tAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGACGGATGGCCtttt                                                                                                                                       <  2:67017/148‑2 (MQ=38)
               aggaggACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCa                                                                                                                                   >  3:905145/1‑149 (MQ=255)
                   ggaCACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATaa                                                                                                                                >  1:315653/1‑148 (MQ=255)
                       acGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTc                                                                                                                                       <  3:163906/138‑1 (MQ=255)
                       acGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTc                                                                                                                                       >  4:163906/1‑138 (MQ=255)
                        cGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGa                                                                                                                           >  1:448569/1‑148 (MQ=255)
                         gTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATGTAGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGTATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATg                                                                                                                         >  4:793626/1‑149 (MQ=255)
                           aTGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGc                                                                                                                        >  1:640766/1‑148 (MQ=255)
                               aaGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGc                                                                                                                   >  3:65346/1‑149 (MQ=255)
                                aGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCa                                                                                                                  <  4:883322/149‑1 (MQ=255)
                                   gggCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGt                                                                                                                >  2:122787/1‑148 (MQ=255)
                                      cAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTcc                                                                                                             <  2:448569/148‑1 (MQ=255)
                                       aaGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCt                                                                                                           <  4:477811/149‑1 (MQ=255)
                                        aGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTa                                                                                                          <  4:905145/149‑1 (MQ=255)
                                           tGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        >  1:143476/1‑148 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTc                                                                                                              >  2:293210/1‑140 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTc                                                                                                              <  1:293210/140‑1 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        >  1:784094/1‑146 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        <  2:784094/146‑1 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACtctc                                                                                                     <  4:758292/149‑1 (MQ=255)
                                                 aGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCAt                                                                                                 <  3:18369/149‑1 (MQ=255)
                                                  gCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                <  3:865425/149‑1 (MQ=255)
                                                   ccGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                <  1:654985/148‑1 (MQ=255)
                                                   ccGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                >  2:807211/1‑148 (MQ=255)
                                                    cGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATggg                                                                                              <  4:64995/149‑1 (MQ=255)
                                                           gTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                       <  4:398639/149‑1 (MQ=255)
                                                           gTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                       >  3:746065/1‑149 (MQ=255)
                                                              gCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcac                                                                                    >  4:102692/1‑149 (MQ=255)
                                                               cTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                                   >  4:868907/1‑149 (MQ=255)
                                                                tGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                                   >  2:135260/1‑148 (MQ=255)
                                                                     cTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                                              >  1:466627/1‑148 (MQ=255)
                                                                     cTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                                              >  1:233611/1‑148 (MQ=255)
                                                                      tGGCATATGTGGGAGTATAAGACGCGCATCGTCGCATCAGGCATTTTCTTCTGCGCCACTGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTTTCGCATGGGGAGACCCCACACTACCAt                                                                            >  4:719709/1‑149 (MQ=255)
                                                                      tGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCAt                                                                            >  4:247031/1‑149 (MQ=255)
                                                                        gCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCt                                                                                                                      <  3:446783/105‑1 (MQ=255)
                                                                        gCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCt                                                                                                                      >  4:446783/1‑105 (MQ=255)
                                                                                gggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc                                                                   >  1:200701/1‑148 (MQ=255)
                                                                                gggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACg                                                                  >  4:833769/1‑149 (MQ=255)
                                                                                 ggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTc                                                                                                                                       >  3:936567/1‑80 (MQ=255)
                                                                                 ggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTc                                                                                                                                       <  4:936567/80‑1 (MQ=255)
                                                                                       taAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                            <  1:655599/148‑1 (MQ=255)
                                                                                        aaGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                          >  3:100757/1‑149 (MQ=255)
                                                                                          gACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCa                                                                                                                                   >  1:595912/1‑74 (MQ=255)
                                                                                          gACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCa                                                                                                                                   <  2:595912/74‑1 (MQ=255)
                                                                                           aCGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                        >  1:208397/1‑148 (MQ=255)
                                                                                           aCGCGCAGCGTCGCATCAGGCATGTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                        >  1:823203/1‑148 (MQ=255)
                                                                                             gcgcAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                     <  4:772860/149‑1 (MQ=255)
                                                                                              cgcAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATACAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                     >  1:793553/1‑148 (MQ=255)
                                                                                              cgcAGCGTCGCATCAGGCATTTTTTTCTGCACCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                    <  4:901080/149‑1 (MQ=255)
                                                                                                cAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                   <  2:793553/148‑1 (MQ=255)
                                                                                                   cGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               <  4:811325/149‑1 (MQ=255)
                                                                                                      cGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                             >  2:168023/1‑148 (MQ=255)
                                                                                                       gCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                            <  2:200701/148‑1 (MQ=255)
                                                                                                          tCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                        <  4:100757/149‑1 (MQ=255)
                                                                                                              gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               <  2:631810/138‑1 (MQ=255)
                                                                                                              gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               >  1:631810/1‑138 (MQ=255)
                                                                                                              gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc                                    >  4:178477/1‑149 (MQ=255)
                                                                                                               cATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTTTGAGTTCGGCATGGGGTc                                    >  2:371415/1‑148 (MQ=255)
                                                                                                                aTTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                                  <  3:532851/149‑1 (MQ=255)
                                                                                                                  ttttttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  4:388581/149‑1 (MQ=255)
                                                                                                                  ttttttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  3:571874/149‑1 (MQ=255)
                                                                                                                   tttttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  2:937245/148‑1 (MQ=255)
                                                                                                                    ttttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg                               <  2:143476/148‑1 (MQ=255)
                                                                                                                    ttttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg                               <  1:135260/148‑1 (MQ=255)
                                                                                                                     tttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg                             <  3:178477/149‑1 (MQ=255)
                                                                                                                             ccAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg                      <  2:709463/148‑1 (MQ=255)
                                                                                                                                aTGCAAAAAGGCCATCCGTCAAGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCgcgc                   >  1:570816/1‑148 (MQ=255)
                                                                                                                                    aaaaaGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACgg              >  4:104239/1‑149 (MQ=255)
                                                                                                                                     gaaaGGCCATCCTGACGGATGGCCTTTTTGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACgg              >  3:511113/2‑149 (MQ=11)
                                                                                                                                       aaGGCCATCCTGACGGATGGCCTTTCGCATAATTTGATGCTTGGGAATTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACTACCGCGCTACGgccg           >  3:808017/1‑149 (MQ=14)
                                                                                                                                         ggCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgccgc          >  1:465812/1‑148 (MQ=255)
                                                                                                                                          gCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCa        <  4:341814/149‑1 (MQ=255)
                                                                                                                                           ccATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCa        >  2:268794/1‑148 (MQ=255)
                                                                                                                                           ccATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAg       <  4:509982/149‑1 (MQ=255)
                                                                                                                                             aTCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAgg      >  2:814500/1‑148 (MQ=255)
                                                                                                                                              tCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgc             >  1:950213/1‑140 (MQ=255)
                                                                                                                                              tCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgc             <  2:950213/140‑1 (MQ=255)
                                                                                                                                                cGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCaaa  <  3:102692/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAA  >  NC_000913/2725898‑2726190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: