Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,397,727:1 +C 100% coding (316/1848 nt) mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,397,7271.C98.3% 2304.7 / 23.1 630coding (316/1848 nt)mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base . (7/4);  new base C (309/309);  total (317/313)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 671 total aligned reads displayed.

TCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTT‑CCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAAC  >  NC_000913/4397579‑4397870
                                                                                                                                                      |                                                                                                                                               
  gTGATAACGGCTGCGGTATCAAAAGAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTT‑cc                                                                                                                                               <  3:1546809/149‑1 (MQ=255)
  gTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTT‑cc                                                                                                                                               <  2:5966520/149‑1 (MQ=255)
    gATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTT‑cccc                                                                                                                                             <  1:2913207/149‑2 (MQ=255)
                 gTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCAc                                                                                                                                 >  4:5737416/1‑149 (MQ=255)
                 gTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCAc                                                                                                                                 <  4:3295362/149‑1 (MQ=255)
                    tCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCAATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTc                                                                                                                              >  4:142204/1‑149 (MQ=255)
                      aaaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACtt                                                                                                                               >  1:141558/1‑146 (MQ=255)
                      aaaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCAc                                                                                                                            <  3:6874964/149‑1 (MQ=255)
                       aaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACg                                                                                                                           >  3:6131527/1‑149 (MQ=255)
                         aaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCa                                                                                                                         <  4:405204/149‑1 (MQ=255)
                            gATGAGCTGGCGCTGGCGCTGGCTCAGCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCCAAGCCATTATTAGCCTCGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGATCGGCTTC‑CCCGCCTGACGCTCACTTCACGCACCg                                                                                                                      <  4:4417706/149‑1 (MQ=255)
                                aGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGa                                                                                                                  >  2:195650/1‑149 (MQ=255)
                                      cgctggcgctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAAcagca                                                                                                            <  3:1174465/149‑1 (MQ=255)
                                          ggcgctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGaa                                                                                                        <  1:1204284/149‑1 (MQ=255)
                                          ggcgctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGaa                                                                                                        <  1:3240572/149‑1 (MQ=255)
                                            cgctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGc                                                                                                      >  2:6437230/1‑149 (MQ=255)
                                              ctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGaa                                                                                                        <  3:3533208/145‑1 (MQ=255)
                                              ctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGaa                                                                                                        >  4:3533208/1‑145 (MQ=255)
                                              ctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCt                                                                                                    <  1:2421307/149‑1 (MQ=255)
                                               tggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTg                                                                                                   >  1:1662943/1‑149 (MQ=255)
                                                ggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTgg                                                                                                  <  2:4895213/149‑1 (MQ=255)
                                                ggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTgg                                                                                                  >  3:930596/1‑149 (MQ=255)
                                                   tCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAg                                                                                               >  1:872671/1‑149 (MQ=255)
                                                    cGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAgg                                                                                              <  4:5322647/149‑1 (MQ=255)
                                                     gTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGc                                                                                             <  4:6112509/149‑1 (MQ=255)
                                                       cATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCt                                                                                           <  1:5727345/149‑1 (MQ=255)
                                                        aTGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTa                                                                                          >  1:4613893/1‑149 (MQ=255)
                                                             accaGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCg                                                                                     >  2:2315684/1‑149 (MQ=255)
                                                               caGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGaa                                                                                   >  1:4167254/1‑149 (MQ=255)
                                                               caGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCGCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGaa                                                                                   <  4:6522763/149‑1 (MQ=255)
                                                                aGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCAc                                                                                                                            <  3:2463792/107‑1 (MQ=255)
                                                                 gTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAgg                                                                                 <  3:5061192/149‑1 (MQ=255)
                                                                  tAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAggg                                                                                >  3:6468944/1‑149 (MQ=255)
                                                                     aaTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGgcgc                                                                             <  1:3602230/149‑1 (MQ=255)
                                                                     aaTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGgcgc                                                                             <  4:1743021/149‑1 (MQ=255)
                                                                     aaTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGgcgc                                                                             >  3:4359089/1‑149 (MQ=255)
                                                                     aaTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGgcgc                                                                             <  1:4923027/149‑1 (MQ=255)
                                                                       tCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGa                                                                           <  1:2691256/149‑1 (MQ=255)
                                                                       tCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGa                                                                           >  4:5184572/1‑149 (MQ=255)
                                                                         gCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGata                                                                         <  2:4425261/149‑1 (MQ=255)
                                                                          ccTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGatat                                                                        >  4:2551586/1‑149 (MQ=255)
                                                                          ccTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGatat                                                                        <  3:2804348/149‑1 (MQ=255)
                                                                          ccTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGatat                                                                        >  1:4334903/1‑149 (MQ=255)
                                                                            tctctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGa                                                                      >  2:7241261/1‑149 (MQ=255)
                                                                             ctctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCGGGAAGCCTGGCAGGCCTATGCCGAAGGGCGTGATGTGaa                                                                     >  4:5997164/1‑149 (MQ=255)
                                                                              tctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAAc                                                                    <  4:2069126/149‑1 (MQ=255)
                                                                                gaggcGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGt                                                                  <  1:2060506/145‑1 (MQ=255)
                                                                                 ggACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCGTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTg                                                                 <  3:4628795/149‑1 (MQ=255)
                                                                                  gACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGa                                                                <  2:5434779/149‑1 (MQ=255)
                                                                                  gACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGa                                                                >  1:5441126/1‑149 (MQ=255)
                                                                                     gATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGaa                                                                                   >  3:32823/1‑127 (MQ=255)
                                                                                     gATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACgg                                                             >  3:1474283/1‑149 (MQ=255)
                                                                                      aTCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGt                                                            <  4:3405002/149‑1 (MQ=255)
                                                                                      aTCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGt                                                            >  2:3143582/1‑149 (MQ=255)
                                                                                           gAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAAc                                                       >  3:5835644/1‑149 (MQ=255)
                                                                                              gCCATTATCAGCCTGGGCTTTCGTGGTGAGGCGCTGGCGTGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTGACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACcgg                                                    >  2:772531/1‑149 (MQ=255)
                                                                                                  ttATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCggcggc                                                >  3:5403853/1‑149 (MQ=255)
                                                                                                   tATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGgcggcg                                               >  2:1536813/1‑149 (MQ=255)
                                                                                                     tCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCa                                             >  3:415895/1‑149 (MQ=255)
                                                                                                       aGCCTGGGCATTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACCTGACGGTAAAACCGGCGGCGCATc                                           <  3:6011467/149‑1 (MQ=255)
                                                                                                        gCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATcc                                          <  2:2633913/149‑1 (MQ=255)
                                                                                                         ccTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCt                                         >  4:1932128/1‑149 (MQ=255)
                                                                                                         ccTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCt                                         <  2:4657340/149‑1 (MQ=255)
                                                                                                           tGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCtgt                                       <  1:594573/149‑1 (MQ=255)
                                                                                                             ggCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTgg                                     <  1:6437230/149‑1 (MQ=255)
                                                                                                              gCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTggg                                    <  3:3930558/149‑1 (MQ=255)
                                                                                                                tttCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGga                                  >  1:5199510/1‑149 (MQ=255)
                                                                                                                tttCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGga                                  >  1:3935357/1‑149 (MQ=255)
                                                                                                                tttCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGga                                  <  2:3579163/149‑1 (MQ=255)
                                                                                                                        tGAGGCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTg                          <  4:88593/149‑1 (MQ=255)
                                                                                                                           ggCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGATGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAg                       <  4:4831740/149‑1 (MQ=255)
                                                                                                                           ggCGCTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAg                       <  3:1853623/149‑1 (MQ=255)
                                                                                                                            gcgcTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAgg                      >  2:877697/1‑149 (MQ=255)
                                                                                                                            gcgcTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAgg                      >  4:2823886/1‑149 (MQ=255)
                                                                                                                              gcTGGCGAGTATCAGTTAGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTg                    <  2:6742099/149‑1 (MQ=255)
                                                                                                                               cTGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGc                   <  4:2197423/149‑1 (MQ=255)
                                                                                                                                tGGCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCt                  <  3:2170241/149‑1 (MQ=255)
                                                                                                                                  gCGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTgg                >  1:1031043/1‑149 (MQ=255)
                                                                                                                                   cGAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGa               >  4:6982823/1‑149 (MQ=255)
                                                                                                                                    gAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAg                       >  1:1823488/1‑140 (MQ=255)
                                                                                                                                    gAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGAt              <  3:662631/149‑1 (MQ=255)
                                                                                                                                    gAGTATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGAt              <  2:1302607/149‑1 (MQ=255)
                                                                                                                                       tATCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAAc                                                                                                                >  1:6695313/1‑48 (MQ=255)
                                                                                                                                        aTCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGt          <  1:719690/149‑1 (MQ=255)
                                                                                                                                         tCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGtt         >  1:3406354/1‑149 (MQ=255)
                                                                                                                                         tCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGtt         <  3:567204/149‑1 (MQ=255)
                                                                                                                                         tCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGtt         <  3:5197241/149‑1 (MQ=255)
                                                                                                                                         tCAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGtt         <  2:5217367/149‑1 (MQ=255)
                                                                                                                                          cAGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTc        >  2:331683/1‑149 (MQ=255)
                                                                                                                                           aGTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCt       >  2:5653520/1‑149 (MQ=255)
                                                                                                                                            gTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTa      >  4:1249301/1‑149 (MQ=255)
                                                                                                                                            gTTCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTa      >  3:3435131/1‑149 (MQ=255)
                                                                                                                                             ttCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAg                       >  1:6901364/1‑131 (MQ=255)
                                                                                                                                             ttCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTac     >  4:3859205/1‑149 (MQ=255)
                                                                                                                                             ttCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTac     >  4:395766/1‑149 (MQ=255)
                                                                                                                                             ttCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTac     >  2:2855301/1‑149 (MQ=255)
                                                                                                                                             ttCGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTac     >  4:6627438/1‑149 (MQ=255)
                                                                                                                                               cGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCt                                                                                           >  2:40951/1‑61 (MQ=255)
                                                                                                                                               cGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGc                                                                                       <  2:2459862/65‑1 (MQ=255)
                                                                                                                                               cGG‑TTTCCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGc                                                                                       >  1:2459862/1‑65 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
TCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGG‑TTT‑CCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAACGTGACGGTAAAACCGGCGGCGCATCCTGTGGGGACGACGCTGGAGGTGCTGGATCTGTTCTACAAC  >  NC_000913/4397579‑4397870

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