Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,007,485 A→G 5.7% K541R (AAA→AGA) ‡ recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,4850AG 5.7% 272.4 / 4.9 87K541R (AAA→AGA) ‡recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (36/46);  new base G (2/3);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACG  >  NC_000913/4007338‑4007632
                                                                                                                                                   |                                                                                                                                                   
tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAAAGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCaaa                                                                                                                                                    >  1:682886/1‑149 (MQ=255)
 aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAAc                                                                                                                                                   <  2:934862/149‑1 (MQ=255)
 aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGAAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAAc                                                                                                                                                   <  1:105275/149‑1 (MQ=255)
  cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACt                                                                                                                                                  <  2:76775/149‑1 (MQ=255)
  cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACt                                                                                                                                                  <  1:372368/149‑1 (MQ=255)
     cTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGAAAACTGCg                                                                                                                                               >  1:206191/1‑149 (MQ=255)
      tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGt                                                                                                                                              >  1:920968/1‑149 (MQ=255)
      tGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGt                                                                                                                                              <  2:1278103/149‑1 (MQ=255)
        aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTaa                                                                                                                                            >  2:1002235/1‑149 (MQ=255)
          cgagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAAt                                                                                                                                          <  1:1552734/148‑1 (MQ=255)
          agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGCTCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAAt                                                                                                                                          >  2:643556/1‑149 (MQ=255)
          agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAAt                                                                                                                                          <  1:1487845/149‑1 (MQ=255)
            cggCCGCGCGCCCGGTGCTGCGCGGCGTATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAATCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttg                                                                                                                                        <  2:634629/148‑5 (MQ=255)
             ggCCGCGCGCCCGGTGCTTCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttgg                                                                                                                                       <  2:334018/149‑6 (MQ=255)
             ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGAGCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttgg                                                                                                                                       <  1:279179/149‑6 (MQ=255)
               ccGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATa                                                                                                                                     <  2:1565126/149‑1 (MQ=255)
                     gcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGAt                                                                                                                               <  1:908041/149‑1 (MQ=255)
                        cGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATCCAGAAATCGTTCGGCGGGAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGaa                                                                                                                            >  1:676636/1‑149 (MQ=255)
                              tgcGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAAt                                                                                                                      <  2:5080/149‑1 (MQ=255)
                                   ggcgAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCATAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTccc                                                                                                                 <  2:821029/149‑1 (MQ=255)
                                     cgAATCCTGTTTGCAACTTGACGTGCCGCGTAGCTTGGAGCTCAAACCGAAAGCTATGCAGAAATCGTTCGGCGGCAACTATGATAGCAAACTGTTCGACAAATTACGCAAACTTATTAAATAGATAGCCGATGAAATTAATGTCccgc                                                                                                               <  2:333921/149‑1 (MQ=255)
                                       aaTCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCccgccg                                                                                                             <  2:644505/149‑1 (MQ=255)
                                         tCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTa                                                                                                           >  2:1134743/1‑149 (MQ=255)
                                           ctctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACg                                                                                                         <  1:206237/149‑1 (MQ=255)
                                           ctctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACg                                                                                                         >  1:1356974/1‑149 (MQ=255)
                                           ctctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGAGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACg                                                                                                         >  2:574400/1‑149 (MQ=255)
                                             ctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCTAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtg                                                                                                       <  1:654648/149‑1 (MQ=255)
                                             ctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtg                                                                                                       <  2:682677/149‑1 (MQ=255)
                                             ctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtg                                                                                                       <  2:339775/149‑1 (MQ=255)
                                             ctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtg                                                                                                       <  2:458146/149‑1 (MQ=255)
                                                tGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCGGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtggtg                                                                                                    <  1:959074/149‑1 (MQ=255)
                                                 gCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGgtgt                                                                                                   <  2:1330639/149‑1 (MQ=255)
                                                 gCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGgtgt                                                                                                   <  2:1558906/149‑1 (MQ=255)
                                                  cAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGAAATGTCCCGCCGTACGTGGTGtt                                                                                                  <  1:1047771/149‑1 (MQ=255)
                                                  cAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTACGCCAAATTACGCAAACTGCGTATATCGATAGCCGATCAAAGTAATGTCCCGCCGTACCTGGTGtt                                                                                                  >  1:268318/1‑149 (MQ=255)
                                                     cTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTaa                                                                                               <  1:532422/149‑1 (MQ=255)
                                                      ttGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAac                                                                                              >  1:1507459/1‑149 (MQ=255)
                                                        gccgtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCATAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAacga                                                                                            <  2:1122769/149‑1 (MQ=255)
                                                         ccgtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAacgac                                                                                           <  2:1551442/149‑1 (MQ=255)
                                                           gtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGc                                                                                         >  1:357386/1‑149 (MQ=255)
                                                           gtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGAAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATTAAAGTAATGTCCCGCCGTACGTGGGGTTTAACGACGc                                                                                         <  2:1545930/149‑1 (MQ=255)
                                                                  gTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttg                                                                                  <  1:696753/149‑1 (MQ=255)
                                                                   tATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttga                                                                                 <  1:742945/149‑1 (MQ=255)
                                                                    aTCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttgat                                                                                <  1:748365/149‑1 (MQ=255)
                                                                           cgcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATg                                                                         <  1:608684/149‑1 (MQ=255)
                                                                             cTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCTATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                       <  1:1084248/149‑1 (MQ=255)
                                                                              tCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCt                                                                      <  2:357386/149‑1 (MQ=255)
                                                                               cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTg                                                                     >  2:1188007/1‑149 (MQ=255)
                                                                                aaaCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGa                                                                    >  2:1561796/1‑149 (MQ=255)
                                                                                 aaCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGaa                                                                   <  2:923163/149‑1 (MQ=255)
                                                                                 aaCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGaa                                                                   <  2:619465/149‑1 (MQ=255)
                                                                                   ccGACAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACa                                                                 <  2:245160/149‑1 (MQ=255)
                                                                                    cGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGTAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAg                                                                >  1:1066405/1‑149 (MQ=255)
                                                                                      aaaGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGAt                                                              >  2:1084409/1‑149 (MQ=255)
                                                                                      aaaGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGAt                                                              >  1:1382672/1‑149 (MQ=255)
                                                                                                aGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCACGCCGTACGTTGTGTTTAACGACGCAACCTTGATTGAGATGGCTGACCAGATGCCGATCAcc                                                    >  1:352190/1‑149 (MQ=255)
                                                                                                 gAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCg                                                   >  1:531011/1‑149 (MQ=255)
                                                                                                 gAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCg                                                   >  1:671593/1‑149 (MQ=255)
                                                                                                  aaaTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGc                                                  <  1:1002235/149‑1 (MQ=255)
                                                                                                  aaaTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGc                                                  >  1:561563/1‑149 (MQ=255)
                                                                                                   aaTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGcc                                                 >  2:1511386/1‑149 (MQ=255)
                                                                                                   aaTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGcc                                                 >  1:388239/1‑149 (MQ=255)
                                                                                                    aTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAAcc                                                                                     >  1:43227/1‑112 (MQ=255)
                                                                                                    aTCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAAcc                                                                                     <  2:43227/112‑1 (MQ=255)
                                                                                                       gTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCg                                             >  1:1353948/1‑149 (MQ=255)
                                                                                                           ggcggcAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAAt                                         >  1:1349917/1‑149 (MQ=255)
                                                                                                            gcggcAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATg                                        >  1:106971/1‑149 (MQ=255)
                                                                                                                  aCTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGc                                  <  2:920968/149‑1 (MQ=255)
                                                                                                                   cTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCAGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCg                                 >  1:955976/1‑149 (MQ=255)
                                                                                                                    tATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATTGCCGATGAAAGTAATGTCCAGCCGTACGTGGTGTTTAACGACGCAACCTTTATAGAGATTGCTGAAAAGATTCCGATCACAGCCAGCGAAATGATCCGCGt                                >  2:1418728/1‑149 (MQ=255)
                                                                                                                    tATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGt                                <  2:1368503/149‑1 (MQ=255)
                                                                                                                    tATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGt                                >  1:646336/1‑149 (MQ=255)
                                                                                                                        aTCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAAc                            >  1:844149/1‑149 (MQ=255)
                                                                                                                         tCGCAAACTGATCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGACGTACGTGGTGTTTAACGAAGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACAGCCAGCGAAATGCTCAGCGTAAACg                           >  1:261226/1‑149 (MQ=255)
                                                                                                                          cGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACgg                          >  1:533122/1‑149 (MQ=255)
                                                                                                                                cTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTgg                    <  2:55320/149‑1 (MQ=255)
                                                                                                                                cTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTgg                    <  2:671593/149‑1 (MQ=255)
                                                                                                                                 tGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTACCGGCGTTggg                   >  1:751829/1‑149 (MQ=255)
                                                                                                                                   ttCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGAt                                                              <  1:1388097/104‑1 (MQ=255)
                                                                                                                                   ttCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGAt                                                              >  2:1388097/1‑104 (MQ=255)
                                                                                                                                   ttCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGAt                 >  2:188844/1‑149 (MQ=255)
                                                                                                                                       ccAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATgcgc             >  1:1408836/1‑149 (MQ=255)
                                                                                                                                          atttACGCAGTTTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAg          >  1:1126845/3‑149 (MQ=255)
                                                                                                                                          aaTTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAg          <  1:1199477/149‑1 (MQ=255)
                                                                                                                                           tttACGCAGTTTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGc         >  1:552204/2‑149 (MQ=255)
                                                                                                                                           aTTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGc         <  2:1080279/149‑1 (MQ=255)
                                                                                                                                           aTTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGc         <  1:657846/149‑1 (MQ=255)
                                                                                                                                                ccAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGaa    <  1:334357/148‑1 (MQ=255)
                                                                                                                                                  aaaCTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACg  >  2:382993/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACG  >  NC_000913/4007338‑4007632

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