Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,397,626 T→C 6.3% L72P (CTG→CCG) ‡ mutL → DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,397,6260TC 6.3% 228.6 / 5.0 79L72P (CTG→CCG) ‡mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base T (43/31);  new base C (2/3);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

CCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGC  >  NC_000913/4397478‑4397770
                                                                                                                                                     |                                                                                                                                                
ccTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                                                                                                                                                 <  1:1270354/150‑1 (MQ=255)
 cTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctgg                                                                                                                                                >  2:807226/1‑150 (MQ=255)
  tGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctggt                                                                                                                                               <  1:664861/150‑2 (MQ=255)
     gTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGTGCTGGCGCTGGCGCTGGCTCg                                                                                                                                            >  1:596657/1‑150 (MQ=255)
     gTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCg                                                                                                                                            >  1:1470847/1‑150 (MQ=255)
        ggTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa                                                                                                                                         >  2:997192/1‑150 (MQ=255)
         gTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCAt                                                                                                                                        <  2:905729/150‑1 (MQ=255)
              cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac                                                                                                                                   >  2:139102/1‑150 (MQ=255)
              cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac                                                                                                                                   >  1:526825/1‑150 (MQ=255)
              cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac                                                                                                                                   <  1:328971/150‑1 (MQ=255)
              cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCCAGCGCCAGCGCCAGCTCATcttttttg                                                                                                                                   <  1:512524/150‑8 (MQ=255)
              cAAAGAACTAGTGGAAAACAGACTCGATGCAGGTGAGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACCGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac                                                                                                                                   >  2:75609/1‑150 (MQ=255)
                aaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccacca                                                                                                                                 <  2:792964/150‑1 (MQ=255)
                      tAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTaaaa                                                                                                                           >  1:726638/1‑150 (MQ=255)
                       aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCCAGCGCCAGCGCCAGCTCATcttttttg                                                                                                                                   >  2:574823/1‑134 (MQ=255)
                       aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCCAGCGCCAGCGCCAGCTCATcttttttg                                                                                                                                   <  1:574823/141‑8 (MQ=255)
                        gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATc                                                                                                                         >  1:1089723/1‑150 (MQ=255)
                        gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACAGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                                                                                                                                       <  1:176979/136‑1 (MQ=255)
                        gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACAGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg                                                                                                                                       >  2:176979/1‑136 (MQ=255)
                          ggAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGc                                                                                                                       >  2:1325108/1‑150 (MQ=255)
                            aaaaCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCct                                                                                                                     <  2:765002/150‑1 (MQ=255)
                                   ccTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAAATTAACCGCATTAGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTggc                                                                                                              >  2:987204/1‑149 (MQ=255)
                                           cAGGTGCTACGCGTATCTATATTGATATCTAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGGTGCGGTATCAAATAAGATGAGCTTGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGTTCTCg                                                                                                      <  2:150728/150‑1 (MQ=255)
                                           cAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCg                                                                                                      >  1:1506101/1‑150 (MQ=255)
                                            aGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGa                                                                                                     >  1:1053073/1‑150 (MQ=255)
                                              gTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAg                                                                                                   <  1:997192/150‑1 (MQ=255)
                                                gCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGcc                                                                                                 >  1:140914/1‑150 (MQ=255)
                                                     cgcgttCGATATCAATATCGATACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTa                                                                                             >  1:1152335/6‑150 (MQ=255)
                                                     cgcgttCGATATCAATATCGATACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATAAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTa                                                                                             >  2:1141979/6‑150 (MQ=255)
                                                      cgTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTAta                                                                                           >  2:174015/1‑149 (MQ=255)
                                                      cgTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCCAGCGCCAGCGCCAGCTcat                                                                                                                                           >  2:690189/1‑100 (MQ=255)
                                                      cgTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCCAGCGCCAGCGCCAGCTcat                                                                                                                                           <  1:690189/102‑3 (MQ=255)
                                                      cgTATCGATATTGATATCCAACGCGGTGGGGCGAAACTTATCCCCATTCGTGATAACGCCTGCGCTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTAta                                                                                           >  2:168677/1‑149 (MQ=255)
                                                       gTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCa                                                                                          <  1:349131/150‑1 (MQ=255)
                                                       gTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCa                                                                                          <  1:1138860/150‑1 (MQ=255)
                                                         aTCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGc                                                                                        >  2:1402786/1‑150 (MQ=255)
                                                         aTCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCGGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGc                                                                                        <  2:1352980/150‑1 (MQ=255)
                                                               atTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGc                                                                                  >  2:1424854/1‑150 (MQ=255)
                                                                ttGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCt                                                                                 >  2:1342146/1‑150 (MQ=255)
                                                                   atatCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCttta                                                                              >  1:429127/1‑149 (MQ=255)
                                                                    tatCGAACGCGGTGGGGAGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGTCTTTcg                                                                             <  1:75609/150‑1 (MQ=255)
                                                                      tCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgcg                                                                           <  2:177647/150‑1 (MQ=255)
                                                                      tCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgcg                                                                           >  2:1235650/1‑150 (MQ=255)
                                                                       cGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCgg                                                                          >  1:636359/1‑150 (MQ=255)
                                                                       cGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCgg                                                                          <  2:895169/150‑1 (MQ=255)
                                                                         aaCGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGTCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTg                                                                        >  1:1444000/1‑150 (MQ=255)
                                                                         aaCGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTg                                                                        <  1:817482/150‑1 (MQ=255)
                                                                         aaCGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCGGGCGCTGGCGCTGGTTCGTGATGCCTCCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTg                                                                        <  2:664634/150‑1 (MQ=255)
                                                                             gcgGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAgg                                                                     <  2:1470847/150‑1 (MQ=255)
                                                                              cgGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGgc                                                                    >  1:1323242/1‑150 (MQ=255)
                                                                                   gggCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTgg                                                               <  2:62835/150‑1 (MQ=255)
                                                                                    ggCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGc                                                              >  2:325166/1‑150 (MQ=255)
                                                                                    ggCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATACCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGc                                                              >  1:1510602/1‑150 (MQ=255)
                                                                                     gCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCg                                                             <  1:325166/150‑1 (MQ=255)
                                                                                       gAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAg                                                           <  2:526825/150‑1 (MQ=255)
                                                                                       gAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAg                                                           <  2:1494892/150‑1 (MQ=255)
                                                                                        aaaCTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGt                                                          <  2:191287/150‑1 (MQ=255)
                                                                                        aaaCTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGt                                                          <  2:185763/150‑1 (MQ=255)
                                                                                           cTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATc                                                       <  1:807226/150‑1 (MQ=255)
                                                                                             ataccGCAGCCGTTATCACGAATGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAg                                                     >  2:1414809/4‑150 (MQ=255)
                                                                                                ccGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTc                                                  <  2:1269755/150‑1 (MQ=255)
                                                                                                    aTTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGtt                                              >  1:112772/1‑150 (MQ=255)
                                                                                                       cGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                                           >  1:1407335/1‑150 (MQ=255)
                                                                                                       cGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                                           >  1:239979/1‑150 (MQ=255)
                                                                                                       cGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                                           >  1:1433854/1‑150 (MQ=255)
                                                                                                       cGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                                           >  1:399926/1‑150 (MQ=255)
                                                                                                                 gCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGc                                 <  2:1506101/150‑1 (MQ=255)
                                                                                                                  cTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCt                                >  2:171088/1‑150 (MQ=255)
                                                                                                                  cTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCt                                >  1:901646/1‑150 (MQ=255)
                                                                                                                    gCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCa                              >  1:1377791/1‑150 (MQ=255)
                                                                                                                      ggTATCAAAAAAGATGAGCTGGAGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCACGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACt                            <  1:1157015/150‑1 (MQ=255)
                                                                                                                          tCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGCTCTCGACCCCAGTATCCGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCAc                        >  1:587808/1‑150 (MQ=255)
                                                                                                                          tCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCAc                        <  2:1152335/150‑1 (MQ=255)
                                                                                                                           cAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGGCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGACTAACAGATCGGTTTCCCGCCTGACGCTCACTTCACg                       >  1:1122691/1‑150 (MQ=255)
                                                                                                                            aaaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGc                      >  1:1453431/1‑150 (MQ=255)
                                                                                                                            aaaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGc                      <  1:1325108/150‑1 (MQ=255)
                                                                                                                               aaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGGCGCTCACTTGACGCAcc                   >  2:1379135/1‑150 (MQ=255)
                                                                                                                                aaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCGGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCg                  <  2:1095665/150‑1 (MQ=255)
                                                                                                                                 aGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGc                 <  2:1089723/150‑1 (MQ=255)
                                                                                                                                   aTGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAg               <  2:1323242/150‑1 (MQ=255)
                                                                                                                                   aTGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAg               >  1:1382574/1‑150 (MQ=255)
                                                                                                                                    tGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGa              >  2:642131/1‑150 (MQ=255)
                                                                                                                                       gctggcgctggcgctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAca           >  1:1414342/1‑150 (MQ=255)
                                                                                                                                       gctggcgctggcgctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAca           >  1:1350431/1‑150 (MQ=255)
                                                                                                                                                ggcgctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGTTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACATCAGTAAGc  <  2:511659/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGC  >  NC_000913/4397478‑4397770

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