Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 16,923 T→G 5.2% K13T (AAG→ACG)  mokC ← regulatory protein for HokC, overlaps CDS of hokC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091316,9230TG 5.2% 270.6 / 5.2 97K13T (AAG→ACG) mokCregulatory protein for HokC, overlaps CDS of hokC
Reads supporting (aligned to +/- strand):  ref base T (36/55);  new base G (2/3);  total (38/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

CAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCT  >  NC_000913/16777‑17070
                                                                                                                                                  |                                                                                                                                                   
cAGCACCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTa                                                                                                                                                   <  2:491984/149‑1 (MQ=255)
 aGCAACCTCGGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGGTGCTTCATTGCTCTCTTCTCCTTGACCTTat                                                                                                                                                  <  2:844231/149‑2 (MQ=255)
 aGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCttcc                                                                                                                                                  >  1:955119/1‑147 (MQ=255)
    aaCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGt                                                                                                                                               <  1:492861/149‑1 (MQ=255)
    aaCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGt                                                                                                                                               <  1:76505/149‑1 (MQ=255)
      ccTCCGTCTGGCCCGTTCGGATGTGAACCTCACAGAGGACTTATATAACAAAAAAAGCCACCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCa                                                                                                                                             <  2:1103257/149‑1 (MQ=255)
      ccTCCGTCTGGCCAGTTCGGCGGTGAACCGCACAGAGGTCTTCCCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCa                                                                                                                                             <  1:629590/149‑1 (MQ=255)
      ccTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCa                                                                                                                                             >  2:170238/1‑149 (MQ=255)
      ccTCCGTCTGGCCAGTTCGGATGTGAACCCCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACCGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGGTTTCTTCTCCTTGACCTTACGGTCa                                                                                                                                             >  1:705091/1‑149 (MQ=255)
      ccTCCGTCTGGCCAGTTCGGATGTCAACCTCACAGAGGTCTTTTCCCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTACTTCATTGCTCTCTTCTCCTTGACCTTACGGTCa                                                                                                                                             <  1:242270/149‑1 (MQ=255)
       cTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAg                                                                                                                                            >  1:756722/1‑149 (MQ=255)
        tCCGTCTGGCCAGTTCGGATGTGAACCTGACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCGTAAGGTCaag                                                                                                                                           <  1:492379/149‑3 (MQ=255)
        tCCGTCTGGCCAGTTCGGATGTGAACCTCACAGATGTCGTTTCTCGGTACCAGCGCCGCCACGACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCGTACGGTCaag                                                                                                                                           <  2:907483/149‑3 (MQ=37)
        tCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCGTAAGGTCaag                                                                                                                                           <  1:492388/149‑3 (MQ=255)
         cggTCGGGCCAGTTCGGATGGGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTGTGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTa                                                                                                                                          <  1:275050/147‑1 (MQ=255)
            tCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAaga                                                                                                                                       <  2:108171/149‑1 (MQ=255)
            tCTGGCCAGTTCGGACGGGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAaga                                                                                                                                       <  1:417184/149‑1 (MQ=255)
             cTGGCCAGTTCGGATGTGACCCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAagag                                                                                                                                      <  1:1113767/149‑1 (MQ=255)
             cTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAagag                                                                                                                                      <  2:664968/149‑1 (MQ=255)
             cTGGCCACTTCCGATGTGAACCTCACAGACGTCTTTTCTCGTTACCAGCGCCGCCACCACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAagag                                                                                                                                      <  2:373339/149‑1 (MQ=255)
               ggCCAGTTCGGATGTGACCCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGc                                                                                                                                    <  2:749287/149‑1 (MQ=255)
                gcgaGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCa                                                                                                                                   <  1:873859/146‑1 (MQ=255)
                gCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCa                                                                                                                                   >  1:808275/1‑149 (MQ=255)
                gCCAGTTCGGAGGGGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCa                                                                                                                                   <  1:873679/149‑1 (MQ=255)
                   agtagggATGAGAACGTCACCGGGGTCCGTGCTCGTTACCAGCGCCGCCATTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCCTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCActc                                                                                                                                <  1:547050/144‑1 (MQ=255)
                   aGTTCGGATGTGAACCTCACAGAGGCCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGCTCAGGGCGACAATCATCGCGTTATGCTGGTTCATTGCCCCCTTCTCCTTTACCTTCCGGTCAGCAAGAGGCActc                                                                                                                                >  2:630165/1‑149 (MQ=255)
                     ttCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTa                                                                                                                              <  2:494410/149‑1 (MQ=255)
                        ggATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATCCAGATGACGATAAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTCCGGTCAGTAAGAGGCACTCTACAt                                                                                                                           >  1:760760/1‑149 (MQ=255)
                        ggATGTCCACCTCACAGAGGTCTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGCTGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACAt                                                                                                                           <  2:1095167/149‑1 (MQ=255)
                             tgAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGtt                                                                                                                      >  2:1117302/1‑149 (MQ=255)
                               aaCCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTCATACAGATGACGATCATTGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCa                                                                                                                    <  2:672009/149‑1 (MQ=255)
                                    cacaGAGGTCTTCTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCata                                                                                                               <  1:397169/149‑1 (MQ=255)
                                              ttttCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCt                                                                                                     >  2:923303/1‑149 (MQ=255)
                                                ttCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGCTCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAATAAGACGCACTCTACATGTGTTAACGATATGGGAGGCCTCg                                                                                                   >  2:310777/1‑149 (MQ=255)
                                                 tctcGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCgg                                                                                                  <  2:629661/149‑1 (MQ=255)
                                                    cGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGtt                                                                                               >  2:1147291/1‑149 (MQ=255)
                                                     gTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTCCGGTCAGTAAGAGGCACTCTACATGTGTTCCGCATGTCGACGGCCTCGGGTTg                                                                                              >  1:1089014/1‑149 (MQ=255)
                                                        aCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATg                                                                                           <  1:905916/149‑1 (MQ=255)
                                                              gccgccACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTaaaa                                                                                     <  2:746772/149‑1 (MQ=255)
                                                                cggcACTACGGCGGTGATACCGATGACGATCAGGGCGACAATCAGCGCCTTATGCTGCTTCATTGCTCTCTGCTCCTTGACTTTACGGTCAGTAAGAGGCACTCTACATGTGTTTAGCATATAGGAGGCCTCGGGTTGATGGTAAAata                                                                                   <  1:299341/146‑1 (MQ=255)
                                                                cgccACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAata                                                                                   <  1:12225/149‑1 (MQ=255)
                                                                 gccACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAatat                                                                                  >  1:491807/1‑149 (MQ=255)
                                                                  ccACGACGGAGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATc                                                                                 <  1:89595/149‑1 (MQ=255)
                                                                    aCTACGGCGGCGATACAGATGACGATCAGGGCGAAAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCAc                                                                               <  1:328395/149‑1 (MQ=255)
                                                                       aCGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCg                                                                            <  1:877229/149‑1 (MQ=255)
                                                                       aCGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCg                                                                            >  2:816845/1‑149 (MQ=255)
                                                                       aCGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCg                                                                            <  2:1058851/149‑1 (MQ=255)
                                                                               gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCtttt                                                                    <  2:756722/149‑1 (MQ=255)
                                                                                 tACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtct                                                                  >  1:13499/1‑149 (MQ=255)
                                                                                 tACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtct                                                                  >  2:102619/1‑149 (MQ=255)
                                                                                 tACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCGCCTTGAGCTTCCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATAAAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGTGCTTTtct                                                                  >  2:231862/1‑149 (MQ=255)
                                                                                    aggtGACGAGCAGGGCCACAATCATCGCCTTATGCTGCTCCATGGCCCTCTCCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCTGGTTGATGGTCAAATATCACTCGGGGCTTTTCTCTa                                                               <  2:483339/146‑1 (MQ=255)
                                                                                    aGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTa                                                               <  2:698975/149‑1 (MQ=255)
                                                                                    aGATGAAGATCAGGGCGACAATCATCGCATTATGTTCCATCACTGCCATGTTCTCCTTGCCATTTTGGTCCGCACGAAGCACTCTACTCGCATTCAGTATCTAGTATGACTCGGGTTGATGGTAGATTATGAATCCGGGCtctttgctc                                                               >  2:373985/1‑144 (MQ=255)
                                                                                      aTGACGATCAGGGCGACAACCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTCATGGTAAAATATCACTCGGGGCTTTTCTCTATc                                                             <  1:841088/149‑1 (MQ=255)
                                                                                       tGACGATCAGGGCGACAAGCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCACCTAGGTGTCCGCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCt                                                            <  2:192589/149‑1 (MQ=255)
                                                                                          tgagcggggCGACAATCATCGCCTTATGCTGGTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTGCATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGcc                                                         <  2:286492/143‑1 (MQ=255)
                                                                                          cGACCAGGCCGACAACCATCGCCTTATGCCGCTTCATTGCTCTCTTCCCCGTGACCTTACGGTCAGTGGCAGGCACGCTACCTGTTGTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGcc                                                         <  1:843457/149‑1 (MQ=255)
                                                                                           gATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCg                                                        >  1:192560/1‑149 (MQ=255)
                                                                                                 ggCGACAATCATCGCCTTATGCTACTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCAGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGc                                                  >  2:791096/1‑149 (MQ=255)
                                                                                                     aaaaccatcGCCTGATGCTGCTTCTGTGCTCGATTCTCCTTGCGCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCAGATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAAt                                              <  2:342324/144‑1 (MQ=255)
                                                                                                     aCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAAt                                              <  1:1022594/149‑1 (MQ=255)
                                                                                                       aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGc                                            >  1:156941/1‑149 (MQ=255)
                                                                                                        atcatcGCCTTATGCTTCTTCATTGCTCTCTGCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATACCTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGcc                                           <  2:248959/149‑1 (MQ=255)
                                                                                                         tcatcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCt                                          <  1:651962/149‑1 (MQ=255)
                                                                                                         tcatcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCt                                          <  1:84239/149‑1 (MQ=255)
                                                                                                         tactgtccTTATGCTTCTTCATTGCTCTCTTCTCCTTGACCCTACGGTCAGGAAGAGGCACTCTACATGTGTTCCGCATATAGGAGGCCTCGTCTTGATGGTAAAAGTTCACTCGGGGCTTTTCTCTATCTGCCGTGCAGCGAATGCCt                                          <  2:174855/143‑1 (MQ=255)
                                                                                                         ataccGCCTTTTGCTGCTTCATTGCTCTCTTCTCCTTGACCTTCCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAAGATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCt                                          <  2:341923/145‑1 (MQ=255)
                                                                                                           atcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTga                                        <  2:20664/149‑1 (MQ=255)
                                                                                                           atcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTga                                        <  1:170238/149‑1 (MQ=255)
                                                                                                            tcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTgag                                       <  2:308265/149‑1 (MQ=255)
                                                                                                              gCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAAAGCCTGagac                                     >  2:37343/1‑149 (MQ=255)
                                                                                                                 ttATGCTGCTTCATTGCTCCCTCCTCCTTGACCTTCCGGTCGGTAAGGGGCACTCTACATGTGTTCAGCATATAGGCGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagacaga                                  <  2:470932/149‑1 (MQ=255)
                                                                                                                      ctgctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTACGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGcc                             <  2:969587/149‑1 (MQ=255)
                                                                                                                      ctgctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTAAGACAGACGGcc                             >  2:832507/1‑149 (MQ=255)
                                                                                                                         ctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCa                          <  2:491807/149‑1 (MQ=255)
                                                                                                                          ttCATTGCTCTCTTTTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGGTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCaa                         <  1:975507/149‑1 (MQ=255)
                                                                                                                          ttCATTGCTCTCTTCTCCTTGACCTTCCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCaa                         <  2:1088218/149‑1 (MQ=255)
                                                                                                                          ttCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCaa                         >  2:1035696/1‑149 (MQ=255)
                                                                                                                           tCAGTGCTCTCTTCTCCTTGACCTTACGGGCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAg                        <  2:210174/149‑1 (MQ=255)
                                                                                                                                  tctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccgc                 <  2:954729/149‑1 (MQ=255)
                                                                                                                                    tcttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccgccg               >  2:532029/1‑149 (MQ=255)
                                                                                                                                    tcttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccgccg               >  1:778087/1‑149 (MQ=255)
                                                                                                                                    tcttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCACCTACTCGCCGAGACAGACAGCCTGCACCAGccgccg               >  1:832842/1‑149 (MQ=255)
                                                                                                                                    tcttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCATGTCTCTATCTGCCGTTCAGCTAATCCCTGAGGCAGACAGCCTCAAGCACccgcca               >  2:696610/1‑148 (MQ=255)
                                                                                                                                    tcttctCCTTGACCTTACGGTCAGTAAGAGGCACATGACATGTGTTCAGCATATAGGAGGGCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccgccg               >  2:818398/1‑149 (MQ=255)
                                                                                                                                      ttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATGGGGGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCt             >  1:793362/1‑149 (MQ=255)
                                                                                                                                      ttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCAGGTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCt             >  1:890194/1‑149 (MQ=255)
                                                                                                                                        cttactGACCGTAAGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGCGGCCTCGGGTTGATGGTACAATATCACTCGGGGCTTTTCTCTCTCTGCCGTTCAGCTGATGCCTGCGGCAGAAAACCTCAAGCAACCGCCGCgat           >  2:633087/6‑149 (MQ=255)
                                                                                                                                        cttactGACCGTAAGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCAGCCGCCGCtat           >  2:1058606/6‑149 (MQ=255)
                                                                                                                                        ctcCTTGACCTACCGGTCACTAAGAGGTCCCATACATGTGTTCAGCATATAGGTTTCCTCGGGTTGATGGTACAATATCACTCGGGGCATTTCTCTATCTGCCGTTCAGCTCATGGCTGAGACAGACAGCCCCAAACAAccgccggcat           >  2:226513/1‑145 (MQ=255)
                                                                                                                                         tcCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCtatt          >  1:249799/1‑149 (MQ=255)
                                                                                                                                         tcCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCtatt          >  1:811562/1‑149 (MQ=255)
                                                                                                                                            ttGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATtata       >  2:870032/1‑149 (MQ=255)
                                                                                                                                             ggcccTTACGGTTAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTatat      <  2:718191/146‑1 (MQ=255)
                                                                                                                                               cccTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCATCTAATGCCTGAGATAGACAGCCTCAAGCACCCGCCGCTATTATATCg    <  1:743757/148‑1 (MQ=255)
                                                                                                                                                 cTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGct  >  1:482374/1‑149 (MQ=255)
                                                                                                                                                 cTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTGTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCCCAAGCACCCGCCGCCATTATATCGca  >  2:1016911/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCT  >  NC_000913/16777‑17070

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