Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,569,437 T→C 5.1% H628R (CAC→CGC)  glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,569,4370TC 5.1% 492.5 / 10.4 138H628R (CAC→CGC) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base T (66/65);  new base C (4/3);  total (70/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 142 total aligned reads displayed.

CAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGATTTAG  >  NC_000913/3569293‑3569582
                                                                                                                                                |                                                                                                                                                 
cAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGcc                                                                                                                                                >  2:871443/1‑148 (MQ=255)
     ttAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTcc                                                                                                                                           <  1:871443/148‑1 (MQ=255)
           gATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGt                                                                                                                                                   >  4:265251/1‑134 (MQ=255)
           gATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGt                                                                                                                                                   <  3:265251/134‑1 (MQ=255)
                       ctAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCcagcag                                                                                                                                    <  2:283363/137‑1 (MQ=255)
                       ctAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCcagcag                                                                                                                                    >  1:487621/1‑137 (MQ=255)
                       ctAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCcagcag                                                                                                                                    <  2:487621/137‑1 (MQ=255)
                       ctAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCcagcag                                                                                                                                    >  1:283363/1‑137 (MQ=255)
                                  aTGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTc                                                                                                              >  4:198775/1‑148 (MQ=255)
                                   tGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCa                                                                                                             >  3:851252/1‑148 (MQ=255)
                                      cTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTATCt                                                                                                               >  4:116176/1‑143 (MQ=255)
                                           tCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGc                                                                                                     >  2:827022/1‑148 (MQ=255)
                                            ccTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGt                                                                                                        >  4:1034377/1‑144 (MQ=255)
                                            ccTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCg                                                                                                    >  1:14741/1‑148 (MQ=255)
                                             cTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGt                                                                                                   >  4:308815/1‑148 (MQ=255)
                                             cTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGt                                                                                                   <  1:827022/148‑1 (MQ=255)
                                              ttttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTcc                                                                                                                                           >  3:718176/1‑107 (MQ=255)
                                              ttttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtct                                                                                                                        >  3:628649/1‑126 (MQ=255)
                                              ttttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtct                                                                                                                        <  4:628649/126‑1 (MQ=255)
                                              ttttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCAGTAATCACTGGGTTATCCTCTCCAGCGAATGGAGGAGTGGCACGCCATTcc                                                                                                                                           <  4:718176/107‑1 (MQ=255)
                                               tttttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTa                                                                                                 >  1:730150/1‑148 (MQ=255)
                                               tttttATCATCTCTGGAACACACACAATCCGGGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAACCGTGGCGTTa                                                                                                 >  4:975977/1‑148 (MQ=255)
                                                        tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtct                                                                                                                        >  2:149209/1‑116 (MQ=255)
                                                        tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtct                                                                                                                        >  2:455788/1‑116 (MQ=255)
                                                        tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtct                                                                                                                        <  1:455788/116‑1 (MQ=255)
                                                        tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtct                                                                                                                        <  1:149209/116‑1 (MQ=255)
                                                        tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTaa                                                                                                <  4:1134788/140‑1 (MQ=255)
                                                        tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTaa                                                                                                >  3:1134788/1‑140 (MQ=255)
                                                        tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAGCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATc                                                                                        >  4:504201/1‑148 (MQ=255)
                                                         ctctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTc                                                                                                              >  3:995303/1‑125 (MQ=255)
                                                           ctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTc                                                                                                                                            >  4:1005539/1‑93 (MQ=255)
                                                           ctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTc                                                                                                                                            <  3:1005539/93‑1 (MQ=255)
                                                           ctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCg                                                                                           >  4:567331/1‑142 (MQ=255)
                                                           ctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCg                                                                                           <  3:567331/142‑1 (MQ=255)
                                                             ggAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaaaa                                                                                    <  2:635715/147‑1 (MQ=255)
                                                             ggAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaaaa                                                                                    >  1:635715/1‑147 (MQ=255)
                                                              gAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTGATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTa                                                                                                 <  4:987666/133‑1 (MQ=255)
                                                               aaCACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGc                                                                                 >  2:379622/1‑148 (MQ=255)
                                                                 cacacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTcccc                                                                                                                                         <  4:786415/90‑1 (MQ=255)
                                                                 cacacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTcccc                                                                                                                                         <  4:946070/90‑1 (MQ=255)
                                                                 cacacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGa                                                                               <  2:14741/148‑1 (MQ=255)
                                                                  acacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTa                                                                                                 <  2:425822/129‑1 (MQ=255)
                                                                  acacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTa                                                                                                 >  1:425822/1‑129 (MQ=255)
                                                                     cacaATCCGTTTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAAtt                                                                                              <  2:1119811/129‑1 (MQ=255)
                                                                     cacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCa                                                                                                             <  1:831398/114‑1 (MQ=255)
                                                                     cacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCa                                                                                                             <  3:475362/114‑1 (MQ=255)
                                                                     cacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAAtt                                                                                              >  1:1119811/1‑129 (MQ=255)
                                                                      acaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaaaaa                                                                                   >  1:632760/1‑139 (MQ=255)
                                                                      acaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaaaaa                                                                                   <  2:632760/139‑1 (MQ=255)
                                                                      acaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaaa                                                                                     <  1:1117253/137‑1 (MQ=255)
                                                                      acaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaaa                                                                                     >  2:1117253/1‑137 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCGTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCgag                                                                         >  3:946449/1‑148 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTcc                                                                                                                                           >  1:1023577/1‑82 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGt                                                                                                                                                   <  3:212355/74‑1 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACaa                                                                                                                           <  4:374323/98‑1 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACaa                                                                                                                           >  3:374323/1‑98 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTg                                                                                                                     >  2:628530/1‑104 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTg                                                                                                                     <  1:628530/104‑1 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGt                                                                                                   <  3:690612/122‑1 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGt                                                                                                   >  4:690612/1‑122 (MQ=255)
                                                                       caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGAt                                                                              >  3:71057/1‑143 (MQ=255)
                                                                               gtgCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTACAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTTCGATCAAAAAGCGATCCGAGAGCAGAAt                                                                 <  3:975977/148‑1 (MQ=255)
                                                                                    ggTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAAt                                                                                               <  2:926011/113‑1 (MQ=255)
                                                                                    ggTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAAt                                                                                               >  1:926011/1‑113 (MQ=255)
                                                                                       cccTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtctc                                                                                                                       >  4:240925/1‑86 (MQ=255)
                                                                                       cccTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtctc                                                                                                                       <  3:240925/86‑1 (MQ=255)
                                                                                       cccTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCg                                                                                           >  3:57862/1‑114 (MQ=255)
                                                                                         cTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAGCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGt                                                       <  3:504201/148‑1 (MQ=255)
                                                                                                     cGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTAcc                                                                                                                <  4:511437/79‑1 (MQ=255)
                                                                                                           cACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTg                                                                                                                     >  4:641792/1‑68 (MQ=255)
                                                                                                            aCTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCa                                                                                                             <  3:887564/75‑1 (MQ=255)
                                                                                                            aCTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCa                                                                                                             >  4:887564/1‑75 (MQ=255)
                                                                                                            aCTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAAt                                                                                               <  1:500687/89‑1 (MQ=255)
                                                                                                            aCTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAAt                                                                                               >  2:500687/1‑89 (MQ=255)
                                                                                                             cTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTc                                   <  2:730150/148‑1 (MQ=255)
                                                                                                             cTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTc                                   <  4:416705/148‑1 (MQ=255)
                                                                                                                   tATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAg                                                                                                           >  2:859502/1‑70 (MQ=255)
                                                                                                                      ccTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGtt                                                                                                                   >  1:669728/1‑59 (MQ=255)
                                                                                                                      ccTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGtt                                                                                                                   <  2:669728/59‑1 (MQ=255)
                                                                                                                         ctcCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGa                                                                               <  4:184813/92‑1 (MQ=255)
                                                                                                                         ctcCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGa                                                                               >  3:184813/1‑92 (MQ=255)
                                                                                                                            cAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCg                                              <  1:331223/122‑1 (MQ=255)
                                                                                                                             aGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGc                                                                      >  3:11873/1‑97 (MQ=255)
                                                                                                                              gCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGtt                                            <  1:988286/122‑1 (MQ=255)
                                                                                                                              gCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGtt                                            >  2:988286/1‑122 (MQ=255)
                                                                                                                              gCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGcc                                       >  1:824474/1‑127 (MQ=255)
                                                                                                                              gCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGcc                                       <  2:824474/127‑1 (MQ=255)
                                                                                                                                     ggggAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCata           >  3:1082584/1‑148 (MQ=255)
                                                                                                                                       ggAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCg                                                                                <  1:254020/77‑1 (MQ=255)
                                                                                                                                       ggAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCg                                                                                >  2:254020/1‑77 (MQ=255)
                                                                                                                                        gAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCata           >  3:647044/1‑145 (MQ=255)
                                                                                                                                         aaTGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAg                     >  4:691869/1‑134 (MQ=255)
                                                                                                                                          aTGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCa                                                                                                             <  3:1123083/45‑1 (MQ=255)
                                                                                                                                           tGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGt                             >  2:869260/1‑124 (MQ=255)
                                                                                                                                           tGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGt                             <  1:869260/124‑1 (MQ=255)
                                                                                                                                           tGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGt                   <  3:326730/134‑1 (MQ=255)
                                                                                                                                            ggCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGAttt    >  3:739910/1‑148 (MQ=255)
                                                                                                                                              cGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGt                   >  2:910294/1‑131 (MQ=255)
                                                                                                                                              cGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGt                   <  1:910294/131‑1 (MQ=255)
                                                                                                                                              cGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGATTTAg  <  4:739910/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGATTTAG  >  NC_000913/3569293‑3569582

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: