Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,308 C→T R81R (CGG→CGA)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,3080CT98.0% 171.7 / NA 52R81R (CGG→CGA) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/22);  total (30/22)
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

CATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATG  >  NC_000913/66062‑66551
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTGCCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCGAGATAGTGGCGTGGCGCGATTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGt                                                                                                                                                                                                                                                 <  1:40167/251‑1 (MQ=255)
            aTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgtat                                                                                                                                                                                                                                     >  1:290911/1‑249 (MQ=255)
              cGATACCCTGTTTTTCGAAGGTTTCTACGATGCCGTTCCCGGTTTCCCACTCAAATTCCCCTTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGTAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgtaacg                                                                                                                                                                                                                                   <  2:244150/251‑5 (MQ=255)
              cGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgtatc                                                                                                                                                                                                                                    <  1:330559/250‑4 (MQ=255)
                    ccTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATCGAGCAGGCGGGGAGTTGGCGTATCggagggg                                                                                                                                                                                                                             >  2:333056/1‑249 (MQ=255)
                    ccTGTGTTTCGAAGGTTTCTACGATGCCGTTCCAGGTTTCCCACTCATAGTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGGGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGgaggag                                                                                                                                                                                                                             <  1:86522/251‑1 (MQ=255)
                        tttttCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCt                                                                                                                                                                                                                         >  1:51646/1‑251 (MQ=255)
                           ttcgaaGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCttt                                                                                                                                                                                                                       <  2:91402/246‑1 (MQ=255)
                                aGGTTTCTAAGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGGGCGAGTGTGTGTGCACAATGCCGCCCATAGACGGGAATTCCTGATCGCGCAGCCGGTGGATTGGCGTATAGGAGGAGGGCTTTTTCgc                                                                                                                                                                                                                  >  2:331684/1‑249 (MQ=255)
                                 ggTTGCTACGATGACGTTAACGGTTTCCCACTCATATTCACCATTGATTTCTGCGTCGGGCATTGTACGGTTGCATGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGTAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCTGCCAATGGACGGGAATGCCTTATAGAGCAGTCGGTGAGTTGTCGTATCGAAGTAGGGCTTTTTCgcac                                                                                                                                                                                                                <  1:321440/251‑4 (MQ=255)
                                            tGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACt                                                                                                                                                                                                     >  1:157564/1‑251 (MQ=255)
                                              aCGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTc                                                                                                                                                                                                   <  2:123996/251‑1 (MQ=255)
                                                                    atTCCCCGTGGATTTCTGCGTCGTTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACg                                                                                                                                                                             <  2:86959/251‑1 (MQ=255)
                                                                              gATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGt                                                                                                                                                                   <  1:224861/251‑1 (MQ=255)
                                                                               aTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTc                                                                                                                                                                  >  1:193167/1‑251 (MQ=255)
                                                                                             cATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAAACACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTAATGatgc                                                                                                                                                    >  2:17398/1‑248 (MQ=255)
                                                                                                tttGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGt                                                                                                                                                 >  1:70166/1‑251 (MQ=255)
                                                                                                      ggTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCg                                                                                                                                            >  1:306847/1‑250 (MQ=255)
                                                                                                         gCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACg                                                                                                                                         <  1:168402/250‑1 (MQ=255)
                                                                                                               aaTGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCCATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGa                                                                                                                                    >  1:161705/1‑249 (MQ=255)
                                                                                                                      ccGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGAt                                                                                                                           <  1:9213/251‑1 (MQ=255)
                                                                                                                       cGTAGAACTAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGCACATATAGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATc                                                                                                                          >  2:222615/1‑251 (MQ=255)
                                                                                                                                 gTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGGCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAcgc                                                                                                                >  2:150347/1‑251 (MQ=255)
                                                                                                                                           gtggtgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTcgcg                                                                                                      >  1:288926/1‑251 (MQ=255)
                                                                                                                                              gtgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATc                                                                                                   >  2:109231/1‑251 (MQ=255)
                                                                                                                                               tgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCCCTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCCCGCGATcc                                                                                                  >  2:374890/1‑250 (MQ=255)
                                                                                                                                               tgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGGTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTcgcgctcc                                                                                                  >  1:81305/1‑247 (MQ=255)
                                                                                                                                                   ggTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACgg                                                                                              >  1:118829/1‑251 (MQ=255)
                                                                                                                                                           gAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGAACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGGTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACg                                                                                      >  1:162281/1‑251 (MQ=255)
                                                                                                                                                                   tggcccGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGTTGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGcccc                                                                               >  1:215944/4‑250 (MQ=255)
                                                                                                                                                                      cccGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCAtgt                                                                           <  1:119522/251‑1 (MQ=255)
                                                                                                                                                                          ccTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGc                                                                       <  2:290911/251‑1 (MQ=255)
                                                                                                                                                                            tGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGt                                                                     >  1:37470/1‑251 (MQ=255)
                                                                                                                                                                                 ccAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGAccc                                                                >  2:234011/1‑250 (MQ=255)
                                                                                                                                                                                   aGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAgg                                                              <  1:17398/251‑1 (MQ=255)
                                                                                                                                                                                     atagtggcgtggcGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGtt                                                            >  1:354415/4‑251 (MQ=255)
                                                                                                                                                                                               ggcgCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGc                                                  <  2:113131/251‑1 (MQ=255)
                                                                                                                                                                                                          gtgtGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAgc                                               >  1:318760/4‑243 (MQ=255)
                                                                                                                                                                                                          gtgtGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAgc                                               <  2:318760/240‑1 (MQ=255)
                                                                                                                                                                                                          gtgggCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGtt                                       <  1:209706/247‑1 (MQ=255)
                                                                                                                                                                                                                    tGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAAt                             >  1:107733/1‑251 (MQ=255)
                                                                                                                                                                                                                     gccgccAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAAt                             >  1:371042/1‑250 (MQ=255)
                                                                                                                                                                                                                     gccgccAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATa                            <  1:246056/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccgccAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCAGGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATAc                           <  1:331684/251‑1 (MQ=255)
                                                                                                                                                                                                                       cgccAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATAcc                          >  2:91676/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCGCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCATCGCCAGGTTGGCCTCTAATACCTGGCg                     <  2:219303/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cgggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTga                <  2:354415/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTg                 >  1:248051/1‑249 (MQ=255)
                                                                                                                                                                                                                                      aTGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGAtct           >  2:6818/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCCTGATCGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGAtcttc         >  1:132054/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cgagCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCCGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATg  >  2:253757/2‑251 (MQ=255)
                                                                                                                                                                                                                                               agagCAGTCGGTGAGGTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATg  <  1:287080/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATG  >  NC_000913/66062‑66551

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