Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,398,679 C→T 5.6% P423L (CCG→CTG) ‡ mutL → DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,398,6790CT 5.6% 281.8 / 4.1 89P423L (CCG→CTG) ‡mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base C (41/42);  new base T (3/2);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

CCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGA  >  NC_000913/4398531‑4398828
                                                                                                                                                    |                                                                                                                                                     
ccGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCg                                                                                                                                                      >  1:1289596/1‑150 (MQ=255)
ccGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCg                                                                                                                                                      >  2:806623/1‑150 (MQ=255)
 cgcgCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGc                                                                                                                                                     <  2:409158/150‑1 (MQ=255)
  gcgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                                                                                                                                    >  2:839730/1‑150 (MQ=255)
   cgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAg                                                                                                                                                   <  2:446600/150‑1 (MQ=255)
   cgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAg                                                                                                                                                   <  1:745846/150‑1 (MQ=255)
    gcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAgg                                                                                                                                                  <  1:375988/150‑1 (MQ=255)
    gcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAgg                                                                                                                                                  >  2:366784/1‑150 (MQ=255)
      tACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGaa                                                                                                                                                >  2:1395956/1‑150 (MQ=255)
        cacTCCTGCGCCAGCATCAGGCAGGCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAAcc                                                                                                                                              <  2:856785/150‑1 (MQ=255)
          cTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTg                                                                                                                                            >  1:78178/1‑150 (MQ=255)
                 gcCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgc                                                                                                                                     >  2:1011974/1‑150 (MQ=255)
                   cAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcgg                                                                                                                                   <  2:1494016/150‑1 (MQ=255)
                    aGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcggc                                                                                                                                  <  1:665373/150‑1 (MQ=255)
                        tCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGCAAGTGCAGGTTCCTGCGGTTccggcgctt                                                                                                                              <  1:305564/150‑4 (MQ=255)
                        gggggCAGCCGGACGACTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAAc                                                                                                                              >  2:413254/4‑150 (MQ=255)
                            gCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTc                                                                                                                          <  1:818548/150‑1 (MQ=255)
                             cAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCa                                                                                                                         >  1:1251083/1‑150 (MQ=255)
                                 cgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGt                                                                                                                     >  1:682407/1‑150 (MQ=255)
                                 cgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGt                                                                                                                     >  2:1417943/1‑150 (MQ=255)
                                 cgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTNNNNAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGt                                                                                                                     >  2:241268/1‑150 (MQ=255)
                                   tcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGttt                                                                                                                   <  2:682407/150‑1 (MQ=255)
                                     cGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTg                                                                                                                 <  2:1152827/150‑1 (MQ=255)
                                      ggCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTgg                                                                                                                <  1:806623/150‑1 (MQ=255)
                                        cTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTc                                                                                                              <  1:1011974/150‑1 (MQ=255)
                                          gCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCgg                                                                                                            >  2:54737/1‑150 (MQ=255)
                                           cccccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCggg                                                                                                           <  2:240748/150‑1 (MQ=255)
                                            ccccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGt                                                                                                          >  2:494624/1‑150 (MQ=255)
                                             cccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                                                                                         <  2:653493/150‑1 (MQ=255)
                                              ccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTAc                                                                                                        >  1:1399577/1‑150 (MQ=255)
                                                  gCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAAGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGAc                                                                                                    >  1:1622833/1‑150 (MQ=255)
                                                     gAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCggg                                                                                                           >  1:469093/1‑140 (MQ=255)
                                                     gAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCggg                                                                                                           <  2:469093/140‑1 (MQ=255)
                                                         gcgcAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTc                                                                                             <  2:1289596/150‑1 (MQ=255)
                                                             aGCCAGGCTACCAGAAACAGCAAGGGGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCAtt                                                                                         <  2:482241/150‑1 (MQ=255)
                                                              gCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCCCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTc                                                                                        <  1:902595/150‑1 (MQ=255)
                                                              accAGCCTACAAGAAAAAGCAAGCTGAACTGTATCGCCACCTTTTGAAAACGCCCGCGCAGATGCAAAAATTAAAAGCGCAGAAACAGAAGGAACCTGCACTTGCGGCGAACAGTGAGAGTTTTGGTCGGGTACTGACTATCGTCCATTc                                                                                        >  2:199743/2‑150 (MQ=255)
                                                               ccAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTcc                                                                                       <  1:324865/150‑1 (MQ=255)
                                                               ccAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTcc                                                                                       >  2:1178500/1‑150 (MQ=255)
                                                                  ggCTACCAGAAACAGCAAGGTGAAGTGTATAGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATGAAAAGCGCCGGAACCGCAGGCACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGCCCATTCCGAc                                                                                    >  2:832859/1‑150 (MQ=255)
                                                                      aCCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtgtg                                                                                <  1:1128607/150‑1 (MQ=255)
                                                                       ccAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGc                                                                               <  1:413254/150‑1 (MQ=255)
                                                                       ccAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGc                                                                               >  2:257582/1‑150 (MQ=255)
                                                                       ccAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGc                                                                               >  1:750258/1‑150 (MQ=255)
                                                                        cAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCg                                                                              >  2:667610/1‑150 (MQ=255)
                                                                          gAAACAGCAAGGTGAAGTGTATCGCCAGCTCTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                                                                                                                                    >  1:1144205/1‑78 (MQ=255)
                                                                          gAAACAGCAAGGTGAAGTGTATCGCCAGCTCTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                                                                                                                                    <  2:1144205/78‑1 (MQ=255)
                                                                           aaaCAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg                                                                           <  1:732002/150‑1 (MQ=255)
                                                                                 cAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGAAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAg                                                                     >  2:192464/1‑150 (MQ=255)
                                                                                    ggTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc                                                                  <  1:528038/150‑1 (MQ=255)
                                                                                       gAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGAc                                                               <  1:70828/150‑1 (MQ=255)
                                                                                            gtATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCaa                                                          <  1:1395956/150‑1 (MQ=255)
                                                                                             tATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAAc                                                         <  2:1402014/150‑1 (MQ=255)
                                                                                             tATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAAc                                                         <  2:1524596/150‑1 (MQ=255)
                                                                                                 gCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACAttt                                                     >  1:1381115/1‑150 (MQ=255)
                                                                                                 gCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACAttt                                                     >  1:1314069/1‑150 (MQ=255)
                                                                                                 gCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTAACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACAttt                                                     <  1:494624/150‑1 (MQ=255)
                                                                                                  ccAGCTTTTGCAACCGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTc                                                    >  2:709153/1‑150 (MQ=255)
                                                                                                      cTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGTGTTGCTGGAGCGCGACGGCAACATTTCACtt                                                >  1:780703/1‑150 (MQ=255)
                                                                                                      cTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACtt                                                >  2:380726/1‑150 (MQ=255)
                                                                                                           gCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATc                                           >  2:1131563/1‑150 (MQ=255)
                                                                                                            cAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATcc                                          >  1:1489607/1‑150 (MQ=255)
                                                                                                            cAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATcc                                          >  1:911918/1‑150 (MQ=255)
                                                                                                              aaCGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCtt                                        >  1:229922/1‑150 (MQ=255)
                                                                                                               aCGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTg                                       <  1:337247/150‑1 (MQ=255)
                                                                                                               aCGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTg                                       <  1:299892/150‑1 (MQ=255)
                                                                                                                cGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGc                                      <  1:667610/150‑1 (MQ=255)
                                                                                                                    cgcgCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGt                                  >  1:870152/1‑150 (MQ=255)
                                                                                                                       gcCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGATGGCAACATTTCACTTTTATCCTTGCCAGTGGc                               <  1:839730/150‑1 (MQ=255)
                                                                                                                       gcCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGc                               <  2:1251700/150‑1 (MQ=255)
                                                                                                                        ccGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCa                              <  2:1374650/150‑1 (MQ=255)
                                                                                                                        ccGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCa                              <  2:270736/150‑1 (MQ=255)
                                                                                                                          gATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGa                            <  2:532766/150‑1 (MQ=255)
                                                                                                                           aTGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGATCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGaa                           <  2:765226/150‑1 (MQ=255)
                                                                                                                             gCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACg                         >  1:1287542/1‑150 (MQ=255)
                                                                                                                             gCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACg                         <  1:1306515/150‑1 (MQ=255)
                                                                                                                              cAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGt                        >  2:360804/1‑150 (MQ=255)
                                                                                                                                aaaaTTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTg                      >  2:205993/1‑150 (MQ=255)
                                                                                                                                 aaaTTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTgg                     <  2:826109/150‑1 (MQ=255)
                                                                                                                                   aTTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCt                   >  2:361001/1‑150 (MQ=255)
                                                                                                                                     actgttcgccgcAAGTGCAGGTTCCTGCGGTTCCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGc                 >  1:1156864/7‑150 (MQ=38)
                                                                                                                                       tgttcgccgcAAGTGCAGGTTCCTGCGGTTCCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGt               >  1:1627871/5‑150 (MQ=38)
                                                                                                                                        aaGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTc              >  2:1042257/1‑150 (MQ=255)
                                                                                                                                        aaGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTc              >  1:203762/1‑150 (MQ=255)
                                                                                                                                        aaGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTc              >  1:209649/1‑150 (MQ=255)
                                                                                                                                         ttcgccgcAAGTGCAGGTTCCTGCGGTTCCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCa             <  2:1156864/148‑1 (MQ=255)
                                                                                                                                         ttcgccgcAAGTGCAGGTTCCTGCGGTTCCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCa             >  2:28417/3‑150 (MQ=255)
                                                                                                                                           cgcCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAgg           <  1:954607/150‑1 (MQ=255)
                                                                                                                                            gcCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGc          <  1:29867/150‑1 (MQ=255)
                                                                                                                                                    cGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGa  >  2:632532/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGA  >  NC_000913/4398531‑4398828

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