Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,961,248 A→T 5.7% L394Q (CTG→CAG)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,961,2480AT 5.7% 303.4 / 7.2 87L394Q (CTG→CAG) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base T (2/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTTTCCA  >  NC_000913/2961109‑2961395
                                                                                                                                           |                                                                                                                                                   
ggAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCTGCTGGCGAt                                                                                                                                            <  1:270403/149‑1 (MQ=255)
ggAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGAt                                                                                                                                            >  1:78822/1‑149 (MQ=255)
ggAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGAt                                                                                                                                            <  2:110157/149‑1 (MQ=255)
ggAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGAt                                                                                                                                            >  1:40736/1‑149 (MQ=255)
   aaTGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCgt                                                                                                                                         <  2:552442/149‑1 (MQ=255)
         atAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTaaaa                                                                                                                                   <  1:645357/148‑1 (MQ=255)
          tAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaac                                                                                                                                  <  1:363935/149‑1 (MQ=255)
                   taACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaacttc                                                                                                                         >  1:187183/1‑149 (MQ=255)
                   taACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaacttc                                                                                                                         >  3:457144/1‑149 (MQ=255)
                    aaCGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAacttcaacttca                                                                                                                        >  2:523114/1‑149 (MQ=255)
                    aaCGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAacttcaacttca                                                                                                                        >  1:580507/1‑149 (MQ=255)
                     aCGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCAc                                                                                                                       >  2:657907/1‑149 (MQ=255)
                       gATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGc                                                                                                                     >  2:403913/1‑149 (MQ=255)
                         tCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTg                                                                                                                   >  2:301161/1‑149 (MQ=255)
                         cccGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCCGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTg                                                                                                                   >  2:594566/2‑149 (MQ=255)
                            gCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCgg                                                                                                                <  1:403913/149‑1 (MQ=255)
                               ggTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCt                                                                                                             <  1:523114/149‑1 (MQ=255)
                                 tGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTAt                                                                                                           <  1:657907/149‑1 (MQ=255)
                                   cgcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATgg                                                                                                         <  4:372481/149‑1 (MQ=255)
                                    gcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGc                                                                                                        <  2:441257/149‑1 (MQ=255)
                                      tACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCaa                                                                                                      >  1:87124/1‑149 (MQ=255)
                                        ccAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCcagcag                                                                                                                                                 >  1:289374/1‑104 (MQ=255)
                                        ccAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAAc                                                                                                    <  4:17294/149‑1 (MQ=255)
                                        ccAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAAc                                                                                                    <  2:368515/149‑1 (MQ=255)
                                        ccAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCAGCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCcagcag                                                                                                                                                 <  2:289374/104‑1 (MQ=255)
                                         cAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACg                                                                                                   >  2:56259/1‑149 (MQ=255)
                                           aaCACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTg                                                                                                 <  2:489286/149‑1 (MQ=255)
                                             cacaGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGa                                                                                               <  3:680033/149‑1 (MQ=255)
                                              acaGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGaa                                                                                              <  4:257483/149‑1 (MQ=255)
                                                  ccTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCTGCTGGCg                                                                                                                                              <  1:333728/97‑1 (MQ=255)
                                                  ccTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCTGCTGGCg                                                                                                                                              >  2:333728/1‑97 (MQ=255)
                                                  ccTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                          >  3:210172/1‑149 (MQ=255)
                                                     gAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATa                                                                                       >  4:223869/1‑149 (MQ=255)
                                                         aaaCGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTTCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGc                                                                                   >  2:549965/1‑149 (MQ=255)
                                                          aaCGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCt                                                                                  >  4:100365/1‑149 (MQ=255)
                                                           aCGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTa                                                                                 >  1:460071/1‑149 (MQ=255)
                                                             ggACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATc                                                                               <  1:442367/149‑1 (MQ=255)
                                                               aCTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCa                                                                             >  3:206603/1‑149 (MQ=255)
                                                                 tGTAGCTGTCGATATCAGCCACCATCATGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGt                                                                           <  3:399347/149‑1 (MQ=255)
                                                                  gTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAGCTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTg                                                                          <  3:345616/149‑1 (MQ=255)
                                                                  gTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTg                                                                          <  3:570363/149‑1 (MQ=255)
                                                                  gTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTATCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTg                                                                          <  2:700629/149‑1 (MQ=255)
                                                                     gCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGAt                                                                       >  1:460206/1‑149 (MQ=255)
                                                                      cTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATc                                                                      >  3:25237/1‑149 (MQ=255)
                                                                       tGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCa                                                                     >  2:13726/1‑149 (MQ=255)
                                                                         tCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAg                                                                   >  2:640839/1‑149 (MQ=255)
                                                                            atgtCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCgg                                                                >  4:287661/4‑149 (MQ=255)
                                                                                 aGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGtt                                                           <  3:374642/149‑1 (MQ=255)
                                                                                  gCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGttt                                                          <  1:467218/149‑1 (MQ=255)
                                                                                  gCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGttt                                                          <  2:65959/149‑1 (MQ=255)
                                                                                  gCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGttt                                                          >  2:353424/1‑149 (MQ=255)
                                                                                   ccaccaTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTg                                                         >  2:167287/1‑149 (MQ=255)
                                                                                   ccaccaTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTg                                                         >  1:41986/1‑149 (MQ=255)
                                                                                    caccaTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGt                                                        <  1:640839/149‑1 (MQ=255)
                                                                                      ccaTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTa                                                      <  3:392862/149‑1 (MQ=255)
                                                                                       caTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTAt                                                     <  1:13726/149‑1 (MQ=255)
                                                                                       caTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTAt                                                     <  1:664085/149‑1 (MQ=255)
                                                                                        aTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATc                                                    <  2:87124/149‑1 (MQ=255)
                                                                                        aTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATc                                                    <  4:210172/149‑1 (MQ=255)
                                                                                         tCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaac                                                                                                                            >  3:315123/1‑76 (MQ=255)
                                                                                         tCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaac                                                                                                                            <  4:315123/76‑1 (MQ=255)
                                                                                         tCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTgg                                                                         <  2:374359/127‑1 (MQ=255)
                                                                                         tCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTgg                                                                         >  1:374359/1‑127 (MQ=255)
                                                                                         tCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCg                                                   >  1:408439/1‑149 (MQ=255)
                                                                                          cACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGc                                                  <  1:206310/149‑1 (MQ=255)
                                                                                             gatgatGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTAc                                               >  3:247970/1‑149 (MQ=255)
                                                                                             gatgatGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTAc                                               >  3:357680/1‑149 (MQ=255)
                                                                                                gatgTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATc                                                    >  3:17933/1‑141 (MQ=255)
                                                                                                gatgTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATc                                                    <  4:17933/141‑1 (MQ=255)
                                                                                                       cgcgGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAActc                                     >  3:358714/1‑149 (MQ=255)
                                                                                                         cgGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCtt                                   <  2:580507/149‑1 (MQ=255)
                                                                                                          ggAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTc                                  <  1:534120/149‑1 (MQ=255)
                                                                                                             gcaaGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGCGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGAt                               >  2:480755/3‑149 (MQ=255)
                                                                                                             gTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGAt                               >  1:554010/1‑149 (MQ=255)
                                                                                                             gTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGAt                               <  2:187183/149‑1 (MQ=255)
                                                                                                              tAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATg                              >  3:625242/1‑149 (MQ=255)
                                                                                                                aGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGtt                            <  1:551247/149‑1 (MQ=255)
                                                                                                                 gtgtCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTc                           >  4:256006/1‑149 (MQ=255)
                                                                                                                   gtCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTcac                         <  1:120280/149‑1 (MQ=255)
                                                                                                                           ttCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCaa                 >  3:620909/1‑149 (MQ=255)
                                                                                                                           ttCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCaa                 >  1:535974/1‑149 (MQ=255)
                                                                                                                               tccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGgc             >  3:360389/1‑149 (MQ=255)
                                                                                                                                 cAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGtt                            <  2:241872/132‑1 (MQ=255)
                                                                                                                                 cAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGtt                            >  1:241872/1‑132 (MQ=255)
                                                                                                                                    gaTCGCCTGCTGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGttt                                                          <  2:490417/98‑1 (MQ=255)
                                                                                                                                    gaTCGCCTGCTGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGttt                                                          >  1:490417/2‑99 (MQ=255)
                                                                                                                                        gCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttc    <  3:223869/149‑1 (MQ=255)
                                                                                                                                        gCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttc    >  2:453784/1‑149 (MQ=255)
                                                                                                                                         ccAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttcc   >  1:39829/1‑149 (MQ=255)
                                                                                                                                          cagcagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttcca  <  1:301161/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTTTCCA  >  NC_000913/2961109‑2961395

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