Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,823,288 C→T G161D (GGC→GAC)  recA ← DNA recombination/repair protein RecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,823,2880CT98.3% 205.1 / ‑6.6 58G161D (GGC→GAC) recADNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/27);  minor base G (0/1);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCA  >  NC_000913/2823046‑2823531
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
cGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACTACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTTGTTTCCGGGTTACCGAACATCACACCATTTTGCAGACGGATCTGGTTGATGAAGATTAGCTGCTGGTTGGACTGCTTCAGGTTGCCGGCCAGCTTACGGATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCg                                                                                                                                                                                                                                             <  1:432239/251‑1 (MQ=255)
         ccacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGTTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGAtt                                                                                                                                                                                                                                    <  1:412859/251‑1 (MQ=255)
             gTTTTCGACCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCg                                                                                                                                                                                                                                >  1:339037/1‑251 (MQ=255)
                    cccTCTTTCCCCGCTCCGATACGACGAATTGCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTATGGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCgg                                                                                                                                                                                                                         <  2:324252/251‑1 (MQ=255)
                                 cgcCGATACGACGTATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCTTTACCCCCGGTATTGGTTTCCGGGTTACCGAACATCACCCCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCTTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTgccgcc                                                                                                                                                                                                            <  1:423820/251‑1 (MQ=255)
                                 cgcCGATACGACGGAAGTCGAGACGAACAGAGGCGTAGAATCTCAGCGCGCTACCACCGGTAGTGGTTTCCGGGTTGCCGAACATCACCCCAATTTTCATCAGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCAGCATAAGTGCCGCAAGGCCCATGGTAGAGGCGTCGATTTCGCCTCCGGTTTCCGCTTTCGGCGACAGTgccgcc                                                                                                                                                                                                            >  2:204510/1‑251 (MQ=255)
                                      aTACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACgg                                                                                                                                                                                                        >  1:430230/1‑250 (MQ=255)
                                           acgGATTTCGGGACGAACAGAGGCGTAGAATTTCACCGCGTTACCACCGGCAGTGGTTTCCGGGTTACCGATCATCACACCAATTTTCATACGGCTCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGATACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCa                                                                                                                                                                                                   >  2:123633/1‑250 (MQ=255)
                                               aTGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAacgac                                                                                                                                                                                              >  1:151964/1‑251 (MQ=255)
                                                      gaCGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACg                                                                                                                                                                                       <  1:369983/251‑1 (MQ=255)
                                                        cGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCACCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTc                                                                                                                                                                                     <  1:35335/251‑1 (MQ=255)
                                                          aaCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCt                                                                                                                                                                                    >  1:445652/1‑250 (MQ=255)
                                                           aCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTAc                                                                                                                                                                                  <  1:22525/251‑1 (MQ=255)
                                                            cAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACt                                                                                                                                                                                 <  2:384495/251‑1 (MQ=255)
                                                                 gCGTAGAATTTCAGCGCGTTACCCCCGGTAGTGTTTTCCGGGTTACCGAACATCACACCAATTTGCATCCGGATCGGGTTGATGAAGATCAGCAGCGTGTTGGACGGCTTCAGGGTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGAGAGAGTCGGCGATTTCGCCTTCGATTTCCGCGTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACGGCGcc                                                                                                                                                                            <  2:404496/251‑1 (MQ=255)
                                                                             aGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGcc                                                                                                                                                                  >  1:330583/1‑249 (MQ=255)
                                                                              gcgcGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAgg                                                                                                                                                               >  1:149037/1‑251 (MQ=255)
                                                                               cgcgTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAggg                                                                                                                                                              >  1:414540/1‑251 (MQ=255)
                                                                                    tACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCAGAACGCGCCAGGGCGTc                                                                                                                                                          >  2:40713/1‑250 (MQ=255)
                                                                                              aGTGGTTTCCGGTTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGGTTTCCa                                                                                                                                               <  1:127934/251‑1 (MQ=255)
                                                                                              aGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCa                                                                                                                                               <  1:353283/251‑1 (MQ=255)
                                                                                                          gTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGc                                                                                                                                   >  2:24646/1‑251 (MQ=255)
                                                                                                                gAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGCCTGCTTCAGGTTACCCGCCAGCTTACGCCTCGCCTGGCTCCTCATACGTGCCGCAAGGCCCATGTGCGAGTCGTCGCTTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGgtgt                                                                                                                             >  1:225575/1‑251 (MQ=255)
                                                                                                                  aCATCCCACCATTTTCCATACGGATCTGGTTGATGAAGATCAGTGGCTTGTTGGACTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCTCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTcc                                                                                                                           <  2:87291/251‑1 (MQ=255)
                                                                                                                  aCATCACACCAATTTTCATACGGATCTGGTTGATGAAGAGCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGCGTCGTCGCTTTCGCCTTCGATTTCCGCTTTCGACGTCAGTGCCGCCACGGAGTCAACGACGCTAACGTCTCCTGCGCCAGAACGCGCCCGGGCGTCCCAGATTTCCAGTTCCTGCGCGCCTGTGTcc                                                                                                                           >  2:308022/1‑251 (MQ=255)
                                                                                                                   cATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCg                                                                                                                          >  1:341537/1‑251 (MQ=255)
                                                                                                                     tCACACCAATTTTCATACGGATCTGGTTGATGACGATAAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGATTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGAGTCAGTGCCGCCACGGAGTCAAAGACGATAACGGCTACTGCGCCAGAACGCGCCAGGGCGCCACAGATTTCCAGGGCCTGCTCACCGGGGTCCGGc                                                                                                                        >  2:197544/1‑251 (MQ=255)
                                                                                                                                 tcagacgtaTCTGGTGGATGAAGATCAGCGGCGTGTTGGACTGCTTCAGGCTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGcac                                                                                                               <  2:12650/240‑1 (MQ=255)
                                                                                                                                          tCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCg                                                                                                   >  2:192587/1‑251 (MQ=255)
                                                                                                                                             gttggatgaAGTTGAGCAGCGGGTGGGCATGCTTCAGTGTCCCCGCCAGGTTACGCATCGCCTGGTGCATCATACGTGCCGCACGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTGCCGCGTTCGGCGTCAGGGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCGCAGATTTCCAGTGGCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTc                                                                                                    <  1:177535/243‑1 (MQ=255)
                                                                                                                                             ggTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGat                                                                                                 <  1:400738/250‑1 (MQ=255)
                                                                                                                                                 gaggaAGATCAGCAGCGTGTGGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAGTTCCAGGGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCg                                                                                             <  1:394732/247‑1 (MQ=255)
                                                                                                                                                    gaAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCATTCGATTTCCGCTTTCGGCGTAAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGACCGGCAGGGCGCACAGCAGGTTGTCGGTCTCGACGc                                                                                         >  2:34896/1‑251 (MQ=255)
                                                                                                                                                         tcagaggcGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAg                                                                                     <  2:341537/244‑1 (MQ=255)
                                                                                                                                                          cagcagCGTGTTGGACTGCTACATGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCCGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGACCGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGTTGTCCGGCTGGGCGCACAGCAGGTTGTCCATATCGAcacacagct                                                                                   >  2:49886/1‑243 (MQ=255)
                                                                                                                                                             cagcGTGTGGGACTGCTGCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGGGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTGCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACGGATTTCCAGTGCCTGCTCGCCTGTGTCCGGCTGGGAGCACAGCGGGTTGGCTATATCGACGCCCAGTTTAc                                                                                <  2:445652/251‑1 (MQ=255)
                                                                                                                                                                        aCTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGAtt                                                                     <  1:7051/251‑1 (MQ=255)
                                                                                                                                                                          tGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTc                                                                                                                                                 >  2:34142/1‑173 (MQ=255)
                                                                                                                                                                          tGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTc                                                                                                                                                 <  1:34142/173‑1 (MQ=255)
                                                                                                                                                                           gCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTGGATTgg                                                                   <  1:123633/250‑1 (MQ=255)
                                                                                                                                                                                aGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTcc                                                               >  1:397000/1‑249 (MQ=255)
                                                                                                                                                                                     aCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCACGTGAGAGTAGTCGCTTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGCGACCGGCTGGGAGCCCAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGcgcgg                                                        >  2:89574/1‑250 (MQ=255)
                                                                                                                                                                                             gCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCa                                                <  1:428742/251‑1 (MQ=255)
                                                                                                                                                                                               ttACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATc                                              >  1:99132/1‑251 (MQ=255)
                                                                                                                                                                                                 aCGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCg                                             >  1:223874/1‑250 (MQ=255)
                                                                                                                                                                                                    cATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGCGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATa                                          >  1:452895/1‑250 (MQ=255)
                                                                                                                                                                                                             gCTCACCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAgg                                <  2:339037/251‑1 (MQ=255)
                                                                                                                                                                                                                catcatTCGTGCCGTAAGGCCCATGGGAGAGTCGTCGATTTCGCCTGCGATGTCCGCGTTCGGCGTCAGTGCCGCCACGGAGGCAACGACGATAACGTCTACTTCGCCAGAACGCGCCAGGTCGTCACAGAGTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCTCAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCTATAAACGCGCAGGt                               <  2:193180/249‑1 (MQ=255)
                                                                                                                                                                                                                catcatACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTAc                                                                                                                                                                                  >  2:144843/1‑102 (MQ=255)
                                                                                                                                                                                                                catcatACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTAc                                                                                                                                                                                  <  1:144843/102‑1 (MQ=255)
                                                                                                                                                                                                                    atACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTAcc                         >  1:23207/1‑251 (MQ=255)
                                                                                                                                                                                                                     tACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCt                        >  1:168514/1‑251 (MQ=255)
                                                                                                                                                                                                                        gTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCa                     <  1:338989/251‑1 (MQ=255)
                                                                                                                                                                                                                        gTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCa                     <  1:192587/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATACCGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCCCAGGTTTTAACTTCAcg                   >  2:314604/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgc                  >  1:417283/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aTGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTgcggc       >  1:405291/1‑251 (MQ=255)
                                                                                                                                                                                                                                           agagTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCa  <  1:40713/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
CGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCA  >  NC_000913/2823046‑2823531

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