Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,891 C→T 7.7% P143P (CCG→CCA glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,8910CT 7.7% 246.0 / 8.1 93P143P (CCG→CCAglgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base C (38/46);  new base T (4/3);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGG  >  NC_000913/3570743‑3571037
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                             <  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: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                                                                          <  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:1036899/1‑149 (MQ=255)
                                                                                                                                     ccACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAAttt               >  1:953431/1‑149 (MQ=255)
                                                                                                                                      cACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTg              <  1:771825/149‑1 (MQ=255)
                                                                                                                                             tGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTAGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGAcgcgcc       >  1:322109/1‑148 (MQ=255)
                                                                                                                                                 cggcggGGCATCATCTTCCCAGGCATCGTGAGCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTGCATTATGCCCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAg   <  2:994871/149‑1 (MQ=255)
                                                                                                                                                  ggcggGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGCCCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAgg  <  1:1085340/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGG  >  NC_000913/3570743‑3571037

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: