Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,398,685 A→T 5.6% E425V (GAA→GTA) ‡ 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,6850AT 5.6% 272.9 / 3.1 85E425V (GAA→GTA) ‡mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base T (3/2);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGACGCC  >  NC_000913/4398537‑4398832
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                               <  2:469093/140‑1 (MQ=255)
                                               gAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCggg                                                                                                               >  1:469093/1‑140 (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                                                                                           >  2:1178500/1‑150 (MQ=255)
                                                         ccAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTcc                                                                                           <  1:324865/150‑1 (MQ=255)
                                                            ggCTACCAGAAACAGCAAGGTGAAGTGTATAGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATGAAAAGCGCCGGAACCGCAGGCACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGCCCATTCCGAc                                                                                        >  2:832859/1‑150 (MQ=255)
                                                                aCCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtgtg                                                                                    <  1:1128607/150‑1 (MQ=255)
                                                                 ccAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGc                                                                                   >  2:257582/1‑150 (MQ=255)
                                                                 ccAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGc                                                                                   >  1:750258/1‑150 (MQ=255)
                                                                 ccAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGc                                                                                   <  1:413254/150‑1 (MQ=255)
                                                                  cAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCg                                                                                  >  2:667610/1‑150 (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:299892/150‑1 (MQ=255)
                                                                                                         aCGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTg                                           <  1:337247/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:209649/1‑150 (MQ=255)
                                                                                                                                  aaGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTc                  >  1:203762/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)
                                                                                                                                               gCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGAc     >  2:1458899/1‑150 (MQ=255)
                                                                                                                                               gCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGAc     >  1:697956/1‑150 (MQ=255)
                                                                                                                                               gCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGAc     >  2:811766/1‑150 (MQ=255)
                                                                                                                                                cAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGACg    >  2:579553/1‑150 (MQ=255)
                                                                                                                                                 aGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGACGc   <  2:1399577/150‑1 (MQ=255)
                                                                                                                                                  ggAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGACGcc  <  2:870152/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCAACATTTCACTTTTATCCTTGCCAGTGGCAGAACGTTGGCTGCGTCAGGCACAATTGACGCC  >  NC_000913/4398537‑4398832

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: