Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,880 G→C 7.2% P147R (CCG→CGG)  glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,8800GC 7.2% 269.6 / 12.3 97P147R (CCG→CGG) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base G (42/47);  new base C (4/3);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

CAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTG  >  NC_000913/3570732‑3571025
                                                                                                                                                    |                                                                                                                                                 
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: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                                                              <  1:1424490/149‑1 (MQ=255)
                                                                                     aTCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATa                                                              >  2:529189/1‑149 (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)
                                                                                                                                                    |                                                                                                                                                 
CAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTG  >  NC_000913/3570732‑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: