Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,885 G→A 7.4% R145R (CGC→CGT glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,8850GA 7.4% 266.4 / 9.4 95R145R (CGC→CGTglgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base G (42/45);  new base A (4/3);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

CGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCG  >  NC_000913/3570737‑3571032
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                     >  2:637326/1‑148 (MQ=255)
                tttCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGATGCCCCGCCGCGCACGCCGTGGGgcagcac                                                                                                                                     <  1:637326/149‑2 (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                                                                                                                        <  1:1473166/149‑1 (MQ=255)
                             atcatcACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATa                                                                                                                        <  2:186373/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:689255/1‑149 (MQ=255)
                                            ccGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACcaccac                                                                                                         <  2:1135371/149‑1 (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                                                                     <  2:405351/149‑1 (MQ=255)
                                                                                aTCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATa                                                                     >  2:529189/1‑149 (MQ=255)
                                                                                aTCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATa                                                                     <  1:1424490/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:953431/1‑149 (MQ=255)
                                                                                                                                           ccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAAttt          >  1:1036899/1‑149 (MQ=255)
                                                                                                                                            cACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTg         <  1:771825/149‑1 (MQ=255)
                                                                                                                                                   tGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTAGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGAcgcgcc  >  1:322109/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCG  >  NC_000913/3570737‑3571032

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: