Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,548,181 A→T 5.8% S630R (AGT→AGA malQ ← 4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,548,1810AT 5.8% 369.6 / 7.7 102S630R (AGT→AGAmalQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base A (50/45);  new base T (3/3);  total (54/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

TTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGCTTCGCCAGTTCGCGATCC  >  NC_000913/3548041‑3548328
                                                                                                                                            |                                                                                                                                                   
ttCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTa                                                                                                                                               >  3:97659/1‑147 (MQ=255)
 tCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACt                                                                                                                                             >  4:667189/1‑148 (MQ=255)
            tcGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTg                                                                                                                                            <  3:101395/138‑1 (MQ=255)
            tcGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGt                                                                                                                                  >  1:1071020/1‑148 (MQ=255)
            tcGGCAAACAGCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTAGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTg                                                                                                                                            >  4:101395/1‑138 (MQ=255)
                cAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAgcg                                                                                                                              >  3:28734/1‑148 (MQ=255)
                   aCATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTg                                                                                                                           >  3:279227/1‑148 (MQ=255)
                      tCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAg                                                                                                                                >  4:970414/1‑140 (MQ=255)
                      tCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAg                                                                                                                                <  3:970414/140‑1 (MQ=255)
                       cGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGTAATGTAgcgc                                                                                                                             >  3:1034869/1‑142 (MQ=255)
                       cGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGTAATGTAgcgc                                                                                                                             <  4:1034869/142‑1 (MQ=255)
                       cGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCaaa                                                                                                                       <  2:1071020/148‑1 (MQ=255)
                       cGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCaaa                                                                                                                       >  1:13415/1‑148 (MQ=255)
                          cTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAg                                                                                                                                <  4:444995/136‑1 (MQ=255)
                          cTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAg                                                                                                                                >  3:444995/1‑136 (MQ=255)
                            cAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACg                                                                                                                  >  3:231845/1‑148 (MQ=255)
                            cAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACg                                                                                                                  <  2:13415/148‑1 (MQ=255)
                              aGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCa                                                                                                                                      >  2:938014/1‑126 (MQ=255)
                              aGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCa                                                                                                                                      <  1:938014/126‑1 (MQ=255)
                                cGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATACAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCCCGGTTc                                                                                                              <  4:28734/148‑1 (MQ=255)
                                  ttGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAg                                                                                                            <  4:231845/148‑1 (MQ=255)
                                   tGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAAcc                                                                                                                     >  3:545065/1‑138 (MQ=255)
                                   tGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAAcc                                                                                                                     <  4:545065/138‑1 (MQ=255)
                                   tGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGc                                                                                                           <  4:279227/148‑1 (MQ=255)
                                       gAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGTGCTCTGTTagg                                                                                                                                             <  3:899129/110‑3 (MQ=255)
                                       gAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGTGCTCTGTTagg                                                                                                                                             >  4:899129/1‑108 (MQ=255)
                                       gAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGt                                                                                                                                                 <  1:1000874/106‑1 (MQ=255)
                                       gAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGt                                                                                                                                                 >  2:1000874/1‑106 (MQ=255)
                                              tGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCAt                                                                                                >  2:89623/1‑148 (MQ=255)
                                                  tCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTc                                                                                                  >  3:760562/1‑142 (MQ=255)
                                                  tCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTc                                                                                                  <  4:760562/142‑1 (MQ=255)
                                                      cAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATTCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGt                                                                                                                                                 >  3:1026498/1‑91 (MQ=255)
                                                      cAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACt                                                                                                                                                   >  1:87203/1‑89 (MQ=255)
                                                      cAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACt                                                                                                                                                   <  2:87203/89‑1 (MQ=255)
                                                      cAGTTATTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGt                                                                                                                                                 <  4:1026498/91‑1 (MQ=255)
                                                         tttttATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGt                                                                                                                <  2:24846/121‑1 (MQ=255)
                                                         tttttATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGt                                                                                                                >  1:24846/1‑121 (MQ=255)
                                                         tttttATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACg                                                                                     <  1:89623/148‑1 (MQ=255)
                                                          ttttATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGa                                                                                    >  3:573650/1‑148 (MQ=255)
                                                                aCTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTg                                                                                                                                            <  4:1094623/86‑1 (MQ=255)
                                                                aCTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTg                                                                                                                                            >  3:1094623/1‑86 (MQ=255)
                                                                   ggTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGAt                                                                                   >  3:158301/1‑140 (MQ=255)
                                                                   ggTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGAt                                                                                   <  4:158301/140‑1 (MQ=255)
                                                                   ggTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGCCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATg                                                                           >  2:989630/1‑148 (MQ=255)
                                                                       aCTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCg                                                                       <  1:989630/148‑1 (MQ=255)
                                                                       aCTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCg                                                                       >  3:368363/1‑148 (MQ=255)
                                                                                 ggAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTg                                                                                                                                                  >  1:321585/1‑63 (MQ=255)
                                                                                 ggAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTg                                                                                                                                                  >  1:755737/1‑63 (MQ=255)
                                                                                 ggAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTg                                                                                                                                                  <  2:755737/63‑1 (MQ=255)
                                                                                 ggAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTg                                                                                                                                                  <  2:321585/63‑1 (MQ=255)
                                                                                        tCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAg                                                                                                            <  1:72787/94‑1 (MQ=255)
                                                                                        tCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAg                                                                                                            >  2:72787/1‑94 (MQ=255)
                                                                                        tCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTa                                                                             <  3:354009/125‑1 (MQ=255)
                                                                                        tCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTa                                                                             >  4:354009/1‑125 (MQ=255)
                                                                                           ccGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGc                                                                                 <  1:483028/118‑1 (MQ=255)
                                                                                           ccGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGc                                                                                 >  2:483028/1‑118 (MQ=255)
                                                                                            cGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGcc                                                                         >  1:823929/1‑125 (MQ=255)
                                                                                            cGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGcc                                                                         <  2:823929/125‑1 (MQ=255)
                                                                                                     cATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGc                                         <  4:573650/148‑1 (MQ=255)
                                                                                                        atCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGTGCTCTGTTa                                                                                                                                               >  3:319923/1‑43 (MQ=255)
                                                                                                        atCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGTGCTCTGTTa                                                                                                                                               <  4:319923/43‑1 (MQ=37)
                                                                                                             gccagTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATc                                 >  3:129229/1‑148 (MQ=255)
                                                                                                              ccagTCTTCCGGCTGTAGTCCTAACAGGGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACGACGGTTCAGCGTCGGCGTGATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAg                                                         >  2:808543/1‑123 (MQ=255)
                                                                                                              ccagTCTTCCGGCTGTAGTCCTAACAGAGCCCTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATcc                                >  3:51007/1‑148 (MQ=255)
                                                                                                              ccagTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAg                                                         <  1:808543/123‑1 (MQ=255)
                                                                                                               cagTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCca                               >  2:544200/1‑148 (MQ=255)
                                                                                                                  tCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCcagca                            >  3:306844/1‑148 (MQ=255)
                                                                                                                   cTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCcagcag                           <  1:544200/148‑1 (MQ=255)
                                                                                                                    ttCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCagcagc                          <  4:306844/148‑1 (MQ=255)
                                                                                                                    ttCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCagcagc                          >  4:232614/1‑148 (MQ=255)
                                                                                                                    ttCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCagcagc                          <  4:129229/148‑1 (MQ=255)
                                                                                                                     tCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGa                                                                                             >  2:453974/1‑80 (MQ=255)
                                                                                                                     tCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGa                                                                                             <  1:453974/80‑1 (MQ=255)
                                                                                                                     tCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACg                                                                  <  4:1004243/107‑1 (MQ=255)
                                                                                                                     tCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACg                                                                  >  3:1004243/1‑107 (MQ=255)
                                                                                                                     tCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATc                                 >  2:575197/1‑140 (MQ=255)
                                                                                                                     tCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATc                                 <  1:575197/140‑1 (MQ=255)
                                                                                                                     tCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGcc                         >  1:1111963/1‑148 (MQ=255)
                                                                                                                       cGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACAt                                                                                          >  2:867899/1‑81 (MQ=255)
                                                                                                                       cGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACAt                                                                                          <  1:867899/81‑1 (MQ=255)
                                                                                                                         gCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCagcagc                          >  2:427813/1‑143 (MQ=255)
                                                                                                                         gCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCagcagc                          <  1:427813/143‑1 (MQ=255)
                                                                                                                          cTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGc                    >  4:808312/1‑148 (MQ=255)
                                                                                                                          cTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGc                    >  3:884015/1‑148 (MQ=255)
                                                                                                                                 ccTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCaaa                                                                                                                       <  2:520875/42‑1 (MQ=255)
                                                                                                                                 ccTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCaa                                                                                       <  3:856030/74‑1 (MQ=255)
                                                                                                                                 ccTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCaa                                                                                       >  4:856030/1‑74 (MQ=255)
                                                                                                                                 ccTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGCTTCGCCa             >  2:1115476/1‑148 (MQ=255)
                                                                                                                                 ccTAACAGAGCACTGTTACTGTCGGCAATGAAGCGCTGCaaa                                                                                                                       >  1:520875/1‑42 (MQ=255)
                                                                                                                                  gtAACAGTGCTCTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTc                                                                                                              <  1:492953/49‑1 (MQ=38)
                                                                                                                                   tAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGCTTCGCCAGt           <  4:884015/148‑1 (MQ=255)
                                                                                                                                    aaCAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCat                                                 >  3:834205/1‑109 (MQ=255)
                                                                                                                                    aaCAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCat                                                 <  4:834205/109‑1 (MQ=255)
                                                                                                                                    aaCAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGCTTCGCCAGtt          >  3:463091/1‑148 (MQ=255)
                                                                                                                                     aCAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGc                                                                                 <  2:1126166/76‑1 (MQ=255)
                                                                                                                                     aCAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGc                                                                                 >  1:1126166/1‑76 (MQ=255)
                                                                                                                                     aCAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGc                                                                                 <  4:285882/76‑1 (MQ=255)
                                                                                                                                     aCAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGc                                                                                 >  3:285882/1‑76 (MQ=255)
                                                                                                                                        gagCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGCTTCGCCAGTTcgcg      <  1:1115476/148‑1 (MQ=255)
                                                                                                                                            aCTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGCTTCGCCAGTTCGCGATcc  <  3:232614/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCAGGAATATTCACCGGTTCGGCCATATCCAGCCAGTCTTCCGGCTGTAGTCCTAACAGAGCACTGTTACTGTCGGCAATGTAGCGCTGCAAACCACGGTTCAGCGTCGGCGTCATCGACATCAACGATGCCTTATGCCCGGCACGTTTCGGCAGACAACCATATTTATGCAGTGCATCCAGCAGCCCTTGCTTCGCCAGTTCGCGATCC  >  NC_000913/3548041‑3548328

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