Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,878 A→G 7.1% W148R (TGG→CGG)  glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,8780AG 7.1% 234.6 / 8.5 99W148R (TGG→CGG) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base A (42/46);  new base G (4/3);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

CCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTG  >  NC_000913/3570731‑3571025
                                                                                                                                                   |                                                                                                                                                   
ccAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCAc                                                                                                                                                    >  1:32217/1‑149 (MQ=255)
 cAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACg                                                                                                                                                   <  2:1521626/149‑1 (MQ=255)
   gcgcGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGc                                                                                                                                                 <  1:791451/149‑1 (MQ=255)
      cgGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTg                                                                                                                                              >  1:823112/1‑149 (MQ=255)
       ggTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTAAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCCTgc                                                                                                                                             >  2:108730/1‑149 (MQ=255)
            tGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCgcggc                                                                                                                                        >  2:1473166/1‑149 (MQ=255)
              cccAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCCCGGCGTGCGcggcgg                                                                                                                                      >  2:131750/1‑149 (MQ=255)
               ccAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCAGCCCCACGGCGTGCGCGGCggg                                                                                                                                     <  1:1223574/149‑1 (MQ=255)
                 aaTTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGc                                                                                                                                   >  1:466589/1‑149 (MQ=255)
                 aaTTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGAGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGTTGCCCCACGGCGTGCGCGGCGGGGc                                                                                                                                   <  1:453866/149‑1 (MQ=255)
                  atagTTTCAAATCGTTGATCATCACCGGAGGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATATCCGGGTGCAAGAACGTTAATCCTTTGACATGCGCTTCATAAATGATCGTGCTGCCCCGCCGCGCACGCCGTGGGGca                                                                                                                                  <  1:1344704/146‑1 (MQ=255)
                    tGTTTCAAATAGTTGATCATCCCCGGAGGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                                                <  2:611623/149‑1 (MQ=255)
                    tGTTTCAAATAGCTGATCATCCCCGGATGCCCGAGGGCTTTCTAAGTGCCACGGATCTCGACCGGGATCTCCGGGGGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                                                <  2:618495/149‑1 (MQ=255)
                    tGTTTCAAAAAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGGTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                                                <  1:299213/149‑1 (MQ=255)
                     gTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCGCCGCGCACGCCGTGGGgcagca                                                                                                                               <  2:951239/149‑1 (MQ=255)
                      tttCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGATGCCCCGCCGCGCACGCCGTGGGgcagcac                                                                                                                              <  1:637326/149‑2 (MQ=255)
                      tttCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGATGCCCCGCCGCGCACGCCGTGGGgcagcac                                                                                                                              >  2:637326/1‑148 (MQ=255)
                       ttCAACTCGTTCATCATCACCGTAGGCCCGAGGGCGTTATAAGTTCGACGGATCTCGACCGGGATATACTGGTGCAAGTACGTTAATCCTTTGACATGGGATTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCatcatc                                                                                                                             <  2:693698/149‑1 (MQ=255)
                       ttCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGCCTTGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGGGGGGCatgatc                                                                                                                             >  1:1506054/1‑149 (MQ=255)
                         cAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCtt                                                                                                                           >  2:771825/1‑149 (MQ=255)
                         cAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTGTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCtt                                                                                                                           >  2:1257891/1‑149 (MQ=255)
                          aaaTAGTTGATCCTCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTc                                                                                                                          <  2:25652/149‑1 (MQ=255)
                          aaaTAGTTGATCATCACCGGATGCCCGAGGGCGTTTTAAGGGCCACGGATCGCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGCCATGGGCTTGATAAATGATGGTGCGGCCCACCGGCGTGCGCGGCGGGGCATCATCTTc                                                                                                                          <  2:1057570/149‑1 (MQ=255)
                           aaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                         <  1:1247075/149‑1 (MQ=255)
                           aaTAGTGGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                                                         <  2:1037367/149‑1 (MQ=255)
                            aTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAGCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTccc                                                                                                                        <  1:1278902/149‑1 (MQ=255)
                             tAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCa                                                                                                                       <  1:612644/149‑1 (MQ=255)
                             tAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGAGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCa                                                                                                                       <  2:171359/149‑1 (MQ=255)
                                ttGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCCCGGCGTGCGCGGCGGGGAATCATGTACCCAGTc                                                                                                                    >  1:1068700/1‑149 (MQ=255)
                                  gATCATCACCGGATGCACGAGGGCTTTATAAGTGCCACGGATCTCCACCAGCATCTCCGGGGGCAAGTACGTTAATCCTTTGACATGGGCTTCATTAATGATGGTGCTGACCCACTGCGTGCGCGGCTGGGCATCAACTTACCAGTCAt                                                                                                                  >  1:848187/1‑149 (MQ=255)
                                   atcatcACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATa                                                                                                                 <  2:186373/149‑1 (MQ=255)
                                   atcatcACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATa                                                                                                                 <  1:1473166/149‑1 (MQ=255)
                                     catcaCCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGt                                                                                                               >  2:961657/1‑149 (MQ=255)
                                        caCCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGAt                                                                                                            >  1:1123332/1‑149 (MQ=255)
                                         aCCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTTATCCTTTGACATGGGCTTCATTAATGATGGTGCTGCCCCACGGCGTGCGCGGAGGGGGATCATCTTCCCAGTCATAGTGATc                                                                                                           >  2:1203915/1‑149 (MQ=255)
                                          ccGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCa                                                                                                          <  1:463767/149‑1 (MQ=255)
                                           cGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCaa                                                                                                         <  2:949572/149‑1 (MQ=255)
                                            ggATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAac                                                                                                        <  1:690704/149‑1 (MQ=255)
                                             gATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATATTCACAGTCATAGTGATCAacc                                                                                                       >  1:554858/1‑149 (MQ=255)
                                             gATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTACGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAacc                                                                                                       <  2:288722/149‑1 (MQ=255)
                                                  ccGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACcaccac                                                                                                  <  2:1135371/149‑1 (MQ=255)
                                                  ccGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACcaccac                                                                                                  >  2:689255/1‑149 (MQ=255)
                                                    gAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATACTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTa                                                                                                >  1:838395/1‑149 (MQ=255)
                                                         cTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCATCCACCACTACGCAt                                                                                           >  2:1129591/1‑149 (MQ=255)
                                                             ataAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCg                                                                                       >  2:1296236/1‑149 (MQ=255)
                                                              taAGTGCCACGGATCTCGACCGGGCTCTCCGGGGGCAAGGACGGTAATCCTTTGACATGGGCTGCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGTGCATCATCTTGCCAGGCATAGTGCTCAACCCCCACTACGCATGTCgg                                                                                      >  2:975316/1‑149 (MQ=255)
                                                                 gTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCCCCACTACGCATTTCGgcgc                                                                                   >  1:522592/1‑149 (MQ=255)
                                                                   gCCACGGATCTCGACGGTGATCTCCGGGGGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGGGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACAACGCATTTCGGCGCaa                                                                                 <  2:141557/149‑1 (MQ=255)
                                                                    ccACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCTCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAAt                                                                                >  2:1117330/1‑149 (MQ=255)
                                                                          aTCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATggcggc                                                                          <  2:239947/149‑1 (MQ=255)
                                                                               gACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCCCTACGCATTTCGGCGCAATGGCGGCGTTgg                                                                     >  2:170849/1‑148 (MQ=255)
                                                                                aCCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTc                                                                    <  2:1123332/149‑1 (MQ=255)
                                                                                  cGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATAAACAACCACTACGCATATCGGCGCAATGGCGGCGTTGTcgc                                                                  >  2:953863/1‑149 (MQ=255)
                                                                                    ggATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGa                                                                <  2:972107/149‑1 (MQ=255)
                                                                                     gATCTCCGGGTGCACGTCCGTTAATCCTTTGACCTGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCACTGGCGGCGTTGTCGCGAt                                                               <  2:568172/149‑1 (MQ=255)
                                                                                     gATCTCCGGGGGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCGGCCCCACGGCGGGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGAt                                                               <  2:539631/149‑1 (MQ=255)
                                                                                      aTCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATa                                                              <  1:1424490/149‑1 (MQ=255)
                                                                                      aTCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATa                                                              >  2:529189/1‑149 (MQ=255)
                                                                                      aTCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATa                                                              <  2:405351/149‑1 (MQ=255)
                                                                                              gTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGtt                                                      >  2:1509738/1‑149 (MQ=255)
                                                                                                gCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCTGCGCTATGGCGGCGTTGTCGCGATAGTCAGGTTCa                                                    >  1:194415/1‑149 (MQ=255)
                                                                                                gCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCa                                                    <  2:1308657/149‑1 (MQ=255)
                                                                                                gCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCa                                                    <  1:1108501/149‑1 (MQ=255)
                                                                                                 cAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAt                                                   <  2:466589/149‑1 (MQ=255)
                                                                                                 cAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCAAGTCATAGTGATCAACCACCACTACGAATTTCGGCGCATTGGAGGCGTTGTCGCGATACTCAGATTCAt                                                   >  1:539846/1‑149 (MQ=255)
                                                                                                  gaGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCCCGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAtt                                                  <  1:965473/148‑1 (MQ=255)
                                                                                                  aaGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAtt                                                  <  2:823112/149‑1 (MQ=255)
                                                                                                  aaGTACGTTAATCCTTTGACATCGGCTTCATAAATGATGGTGCTGCCCCACGGCGGGCGCGGCGGGGCATCATCTTCCCAGTCTTAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCAtt                                                  <  1:129498/149‑1 (MQ=255)
                                                                                                   aGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTa                                                 >  1:314849/1‑149 (MQ=255)
                                                                                                        gTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGAcc                                            >  1:294215/1‑149 (MQ=255)
                                                                                                        gTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCGTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGAcc                                            >  2:719451/1‑149 (MQ=255)
                                                                                                          tAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCgg                                          >  1:109757/1‑149 (MQ=255)
                                                                                                           aaTCCTTTGACATGGGCTGCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGGTCATTATGACCGGc                                         <  1:9797/149‑1 (MQ=255)
                                                                                                                    aCATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGcagcag                                >  2:26888/1‑149 (MQ=255)
                                                                                                                     cATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCagcagc                               >  1:540783/1‑149 (MQ=255)
                                                                                                                       tGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCgg                             <  1:1257891/149‑1 (MQ=255)
                                                                                                                       tGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGAATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACAGGCGTGCAGCAGCgg                             >  2:1349482/1‑149 (MQ=255)
                                                                                                                        gggCTTCATAAATGATGGAGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCggg                            <  1:26888/149‑1 (MQ=255)
                                                                                                                          gCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGtt                          >  2:859113/1‑149 (MQ=255)
                                                                                                                             tCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGTCATCATCTTCCCAGTCATAGTGATCACCCTCCACTACGCATTTCGGCGCATTGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATc                       <  1:170849/149‑1 (MQ=255)
                                                                                                                              cATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCAGTTCGGCGCAATGGCGGCGTTGTCGCGATTGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCt                      >  2:194990/1‑149 (MQ=255)
                                                                                                                                tAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCttt                    <  2:364959/149‑1 (MQ=255)
                                                                                                                                 aaaTGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGTCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTa                   <  2:869783/149‑1 (MQ=255)
                                                                                                                                 aaaTGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTa                   <  2:966423/149‑1 (MQ=255)
                                                                                                                                 aaaTGATGGTGAGTCCCCACGGCGTTCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTa                   <  2:367395/149‑1 (MQ=255)
                                                                                                                                  aaTGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTGGCGATAGTGAGGTTCATTATGTCCGGCGTGCTGCTGCGGGTTATCTTTaa                  >  1:367395/1‑149 (MQ=255)
                                                                                                                                   atgatgGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAAGCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTCTGACCGGCGCGCAGCAGCGGGTTATCTTTaaa                 <  1:295047/149‑1 (MQ=255)
                                                                                                                                   atgaggCTGCTGCCAGAAGGCGTGCGCGGCGGTGAAGCATCTTCCCAGTCATAGTTAGCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGGAGCTATAGTCAGGTGCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTaaa                 <  2:1094207/149‑1 (MQ=255)
                                                                                                                                      aagatgatgcCCCGCCGCGCACGCCGTGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTc              >  2:1013172/8‑149 (MQ=255)
                                                                                                                                       tgGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTcc             <  2:294215/149‑1 (MQ=255)
                                                                                                                                        gatgatgcCCCGCCGCGCACGCCGTGGGGCATCATCTTCCCAGTCATAGTTATCAAACACCACTACGCATTTCGGCGCAATGGCGGCGTTGTAGCGATTGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTccc            >  2:344568/6‑149 (MQ=255)
                                                                                                                                        gatgatgcCCCGCCGCGCACGCCGTGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGAGAGTCAGGGTCATTATGACCGGCGTGCAGCCGCGGGTGATCTTTAAACTccc            >  1:478309/6‑149 (MQ=255)
                                                                                                                                         gTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTcccc           >  1:1367423/1‑149 (MQ=255)
                                                                                                                                               ccccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAAt     >  1:269170/1‑149 (MQ=255)
                                                                                                                                               caccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAAt     >  1:954113/3‑149 (MQ=255)
                                                                                                                                                 ccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTCTCTTTAAACTCCCCGTCAAttt   >  1:1317801/1‑149 (MQ=255)
                                                                                                                                                 ccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAAttt   >  1:1036899/1‑149 (MQ=255)
                                                                                                                                                 ccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAAttt   >  1:953431/1‑149 (MQ=255)
                                                                                                                                                  cACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTg  <  1:771825/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTG  >  NC_000913/3570731‑3571025

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