Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,823,589 G→A R61C (CGT→TGT)  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,5890GA97.7% 148.1 / NA 43R61C (CGT→TGT) recADNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/22);  minor base C (1/0);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATA  >  NC_000913/2823345‑2823830
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
aTAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATCCAGCCCAt                                                                                                                                                                                                                                             >  1:285727/1‑251 (MQ=255)
               ccAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGc                                                                                                                                                                                                                              >  1:297227/1‑251 (MQ=255)
               ccAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCAACGCAGGATGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCAAAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACATGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTCGATTTCGACGATACCGCCTATCGGCAGACCACCTGc                                                                                                                                                                                                                              >  2:198814/1‑251 (MQ=255)
                      gcgcCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAgc                                                                                                                                                                                                                       >  1:83067/1‑251 (MQ=255)
                          cAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGcgcg                                                                                                                                                                                                                    >  2:301717/1‑250 (MQ=255)
                             ggCGTCAACGATTTCCAGAGCCTGATCCCCGGTGTCCGGCTGGGCGCACAGCAGGTTGTCGCTATCGCAGCCCAGTTTACGTGCGTAGCTTGGGTCCAGCGCGTGTTCAGCATCGATAAACTCACAGGTTTTACCTTCACGCTGCGCTGCCGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATGTCGACGATACAGCCCCTCTGCAGACCAACTGCCCCAAGCGCGatat                                                                                                                                                                                                                >  2:219617/1‑251 (MQ=255)
                                gTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGatatac                                                                                                                                                                                                              >  1:125043/1‑248 (MQ=255)
                                                    gCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTagag                                                                                                                                                                                         <  1:275602/251‑1 (MQ=255)
                                                       cGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGat                                                                                                                                                                                                                  <  2:244456/223‑1 (MQ=255)
                                                       cGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGat                                                                                                                                                                                                                  >  1:244456/1‑223 (MQ=255)
                                                         ccGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGt                                                                                                                                                                                    <  1:253931/251‑1 (MQ=255)
                                                              gtCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTcc                                                                                                                                                                                <  1:43658/250‑1 (MQ=255)
                                                                                        cGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAgg                                                                                                                                                     <  1:182217/251‑1 (MQ=255)
                                                                                           tatCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGgcgc                                                                                                                                                  >  1:165871/1‑251 (MQ=255)
                                                                                              cGACGCCCAGTGTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCatg                                                                                                                                               <  2:316887/251‑1 (MQ=255)
                                                                                                     cAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGAGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCACGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCCTGATGGAGc                                                                                                                                        >  1:24737/1‑251 (MQ=255)
                                                                                                                                  cgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCTCCTGCCCTCAGCGCGATATCTAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGTAACgttc                                                                                                                                                               <  1:257805/199‑4 (MQ=255)
                                                                                                                                  cgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTc                                                                                                                                                               <  1:257807/199‑1 (MQ=255)
                                                                                                                                  cgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTc                                                                                                                                                               >  2:257807/1‑199 (MQ=255)
                                                                                                                                  cgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCGGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGccc                                                                                                           <  1:76981/251‑1 (MQ=255)
                                                                                                                                        cAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCAGCACAGTAAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGACCCGGAAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCACTCTGGCCCAGtgct                                                                                                     >  2:229909/1‑251 (MQ=255)
                                                                                                                                           cATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCCGCGTGGTTTTACCGGAAGATTACGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATAGTTTCTCAATCTGGCCCAGTGCTgcc                                                                                                  >  1:20257/1‑251 (MQ=255)
                                                                                                                                                  aaaCGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTTTTTCTCAATCTGGCCCAGTGCTGCCGCCAAc                                                                                            <  1:364002/250‑1 (MQ=255)
                                                                                                                                                     cGCACAGTTTTGACCTTCACGCTGCGCTGCGGCGACCCCCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTGTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCtt                                                                                        <  2:206383/251‑1 (MQ=255)
                                                                                                                                                       cacaGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCCGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCCTCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTc                                                                                      >  2:84419/1‑251 (MQ=255)
                                                                                                                                                                            gctctgcgGCGATCACCGGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCAACTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAAAGCTTTCTGTTTGTGTTCGTCGATAg                                                                   <  2:274092/246‑1 (MQ=255)
                                                                                                                                                                             cgctgcgGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGCCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTGCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCGTTTTGTTTTCGTCGATAGCCa                                                                <  2:364874/251‑1 (MQ=255)
                                                                                                                                                                                         cACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTg                                                    <  2:195607/251‑1 (MQ=255)
                                                                                                                                                                                          aCCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGt                                                   >  1:41064/1‑251 (MQ=255)
                                                                                                                                                                                           ccTGCAGCGTCAGCGTGGTTTCACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTc                                                  >  1:96102/1‑251 (MQ=255)
                                                                                                                                                                                            cTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTc                                                  >  2:278521/1‑250 (MQ=255)
                                                                                                                                                                                              tcagcgtcagcgtGGTTTTACCGGAATATTCCGGTCCGTAGAATTCGACGATACAGCCCATCGGCAGCCCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCa                                                 <  2:195493/248‑1 (MQ=255)
                                                                                                                                                                                              gCAGCGTCAGCGTGTTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGAGACAGCCCAACGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCGGGCCCAGTGCTGCCGCCAACGCTTTCGGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCa                                                 <  2:24737/249‑1 (MQ=255)
                                                                                                                                                                                                 gcgtcagcgtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCg                                            <  1:195386/251‑1 (MQ=255)
                                                                                                                                                                                                 gcgtcagcgtGGTGTTACAGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGTCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGACGATAGCCATTTTTACTCCTGTCATGCCg                                            <  2:115882/251‑1 (MQ=255)
                                                                                                                                                                                                  cgtcagcgtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCgg                                           >  1:188687/1‑251 (MQ=255)
                                                                                                                                                                                                  cgtcagcgtGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCgg                                           <  1:301717/251‑1 (MQ=255)
                                                                                                                                                                                                           ggTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCg                                  <  1:278521/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCAATCGGCAGACCAACTGCCCCAAGAGCGATATCCAGTGAAAGCGAACAGGTAGAGATGGTTTCCACCTCCATGGAACGGTCTTCACCCAGGCGCCTGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCACCGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGCTAGTc                           >  1:245959/1‑251 (MQ=255)
                                                                                                                                                                                                                   cGGAAGATTCCGGTCCGTAGATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTc                           >  2:352169/1‑250 (MQ=255)
                                                                                                                                                                                                                       aGATTCCGGTCCGTAGATTTAGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCACATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGCTATCCATTTTTACTCCTGTCCTGCCGGGTAATACCGGATAGCCAatat                      >  2:129118/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gATTTCGACGATACAGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCaa       <  2:220843/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cgacgaTACAGCAAATCGGAAGACCACCTTCCCCAAGCGCTATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGCGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCAATAGCCATTTTTACTCCTGTCATGCCGGGTAATTCCGGATAGTCAATATGTTCTGTTGAAGCAATtata  <  2:83067/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
ATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATA  >  NC_000913/2823345‑2823830

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