Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,397,615 G→A 6.4% L68L (CTG→CTA) ‡ 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,6150GA 6.4% 242.0 / 4.3 79L68L (CTG→CTA) ‡mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base G (43/30);  new base A (2/3);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

TGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACA  >  NC_000913/4397467‑4397761
                                                                                                                                                     |                                                                                                                                                  
tggtCGAGCGACCTGCGTCGGTAGTCAAATAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGATCGCGGTGGGGCGAATCTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgg                                                                                                                                                   <  2:497848/150‑1 (MQ=255)
 ggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggc                                                                                                                                                  <  1:720599/150‑1 (MQ=255)
   tCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgc                                                                                                                                                <  2:999111/150‑1 (MQ=255)
        cGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcg                                                                                                                                           <  1:146134/150‑1 (MQ=255)
           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:1235650/1‑150 (MQ=255)
                                                                                 tCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTcgcg                                                                  <  2:177647/150‑1 (MQ=255)
                                                                                  cGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCgg                                                                 <  2:895169/150‑1 (MQ=255)
                                                                                  cGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCgg                                                                 >  1:636359/1‑150 (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:239979/1‑150 (MQ=255)
                                                                                                                  cGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                                  >  1:1433854/1‑150 (MQ=255)
                                                                                                                  cGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                                  >  1:399926/1‑150 (MQ=255)
                                                                                                                  cGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTcc                                  >  1:1407335/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:1325108/150‑1 (MQ=255)
                                                                                                                                       aaaaaaGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGc             >  1:1453431/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                  
TGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACA  >  NC_000913/4397467‑4397761

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: