Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,820 C→T G129R (GGA→AGA)  thiE ← thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,8200CT98.1% 164.6 / ‑5.2 52G129R (GGA→AGA) thiEthiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/20);  minor base G (0/1);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

ATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTC  >  NC_000913/4193575‑4194060
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
atCGCCAACTCCTGCAATTTCCAGCAACTGTGCCTTGGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGCCGGCGATACTGCCGACACCCGTTGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGGGGGATAATCCGCCAGGCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAgcg                                                                                                                                                                                                                                              <  2:96702/251‑1 (MQ=255)
atCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGACGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACCTGTCTCAgcg                                                                                                                                                                                                                                              >  1:184678/1‑251 (MQ=255)
     cAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATCCTGACGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGat                                                                                                                                                                                                                                          >  1:305161/1‑250 (MQ=255)
        cTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGAtata                                                                                                                                                                                                                                       >  2:39603/1‑250 (MQ=255)
          ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGAtata                                                                                                                                                                                                                                       >  2:76824/1‑248 (MQ=255)
          ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGCGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGCAGGCATCTGTTTGGTTTGCGTCGGGCACACATGTCTCCGCGCGagctaaa                                                                                                                                                                                                                                     >  2:186962/1‑244 (MQ=255)
          ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAAAGGCGCTGAAGACGGCGATCCTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCCGACTGATGACGCCAATCGCCACGGTGGGGTAATCCGCCAGTCGCTAAGCATGGCGTGCCAGCTGTTACAGCCCCTGAGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGAtata                                                                                                                                                                                                                                       >  1:26007/1‑248 (MQ=255)
                    ccAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCCGCCCCCGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGAGGCGtgct                                                                                                                                                                                                                          >  2:86795/1‑251 (MQ=255)
                          aCTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTGGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAgc                                                                                                                                                                                                                    <  1:39603/251‑1 (MQ=255)
                             gtgCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGCCGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGcgcg                                                                                                                                                                                                                 >  1:108678/1‑251 (MQ=255)
                              tgCCGTTGCCAAACGCCAGTCCGCGGCTTAAGTAATGGAGCTGACGACGACGATAATGCCGACACCCGTTGCTATCACCGCAGGCGCACGTGCCAGACTGATTCCGCCAATCGCCACGGTGGGATAATCCGCCAGCCCCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGCGCAGAAGGCATCTGTTCGGTTTGCGTCGGGAACACATGTCTCAGCCCGATAAAAACGGGGCGcgctgacaacccg                                                                                                                                                                                                                 >  2:280752/1‑242 (MQ=255)
                                  gTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTc                                                                                                                                                                                                            <  1:189432/251‑1 (MQ=255)
                                   ttGCCAACCGCCAGTCTGCGGCTTGAGTAATGGCGCTGCCGACGGCGATACTTCCAACACCCGTTCCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGTGAACACATGTCGCATCGCGATATAAGAGGGGCGTTCTGCCAGCGCGACGTc                                                                                                                                                                                                            <  2:306926/250‑1 (MQ=255)
                                    tGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCg                                                                                                                                                                                                           <  1:202062/250‑1 (MQ=255)
                                       cAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAttt                                                                                                                                                                                                       >  1:69343/1‑251 (MQ=255)
                                          aCGCCAGTCTGCGGCTTGCGTAATGGCGCTGACGACGAAGATACTGCCGCCACACGTTGCTATCACCGCAGGCGCGCGTACCACACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGCGCCAGCTGTTCCAGCCCCTGCGGCGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGCGCTGCCAGCGCGACGTCGATTTccc                                                                                                                                                                                                    >  1:201318/1‑250 (MQ=255)
                                                gTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGt                                                                                                                                                                                              >  1:41364/1‑251 (MQ=255)
                                                             gTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAAc                                                                                                                                                                                <  1:293306/251‑1 (MQ=255)
                                                                     gcTGACGCCGGCGCTACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATCTAAGAGTGGCGTGCTGCCAACGCGACGTCGATTTCCATATCGTCATGGTGCGAAACGCCCAGcc                                                                                                                                                                        >  2:273903/1‑251 (MQ=255)
                                                                       tGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGc                                                                                                                                                                      >  1:243159/1‑251 (MQ=255)
                                                                        gacgacGGCGATACTGCCGACCCCCGTTGCTATCACCCCAGACGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGcctcc                                                                                                                                                                     >  2:132385/1‑248 (MQ=255)
                                                                                    aCTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCcgcgc                                                                                                                                                         >  1:218008/1‑251 (MQ=255)
                                                                                         cGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATg                                                                                                                                                    <  2:277719/251‑1 (MQ=255)
                                                                                          gACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATgg                                                                                                                                                   <  1:263616/251‑1 (MQ=255)
                                                                                                   ggCTATCACCGCAGGCGCGCGTGACAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGATTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTgaga                                                                                                                                           <  2:211321/249‑1 (MQ=255)
                                                                                                    gCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTg                                                                                                                                              >  2:281326/1‑246 (MQ=255)
                                                                                                      tATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCCTTg                                                                                                                                              >  2:61407/1‑244 (MQ=255)
                                                                                                                            aGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGCTTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGccc                                                                                                                 >  2:74942/1‑251 (MQ=255)
                                                                                                                             gACTGATGCCGCCAACCGCCCCGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGcccc                                                                                                                <  1:163950/251‑1 (MQ=255)
                                                                                                                                    gccGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGa                                                                                                         <  2:302848/250‑1 (MQ=255)
                                                                                                                                       gCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGc                                                                                                      <  2:127433/251‑1 (MQ=255)
                                                                                                                                         cAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGCGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGcc                                                                                                     >  1:259286/1‑250 (MQ=255)
                                                                                                                                               ccACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATAGGcc                                                                                              >  1:302337/1‑251 (MQ=255)
                                                                                                                                                           aaTCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGAt                                                                                   <  1:69057/250‑1 (MQ=255)
                                                                                                                                                                      cGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAgccgcc                                                                        >  1:64362/1‑250 (MQ=255)
                                                                                                                                                                         tCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGt                                                                     >  1:11160/1‑250 (MQ=255)
                                                                                                                                                                           aaCATGCCGTGCCAGCTGTTCCAGCCACTGCGGAGCAGACGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGCGGGGCGTGCTGCAAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGAAAt                                                                  >  2:127475/1‑251 (MQ=255)
                                                                                                                                                                            aCATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATc                                                                 <  1:185693/251‑1 (MQ=255)
                                                                                                                                                                                  cGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAAGCATCTCTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGAAGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCAATACGACTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGAt                                                           >  1:142302/1‑251 (MQ=255)
                                                                                                                                                                                              tCCAGCCCCGGCGGTGCAGAAGGCATTTGTTTGGTTTGCGTCGGGAACACAGTTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATcgcg                                                <  2:243159/250‑1 (MQ=255)
                                                                                                                                                                                                     ccTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAg                                        <  1:137083/251‑1 (MQ=255)
                                                                                                                                                                                                           gTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGCGGGGCGTGCTGCCAGCGCGACGTCGCTTTCCATATCGTCATGTGTCGAAACGCCCCGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCCAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGcggcgg                                  >  2:200698/1‑251 (MQ=255)
                                                                                                                                                                                                           gTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGcggcgg                                  >  1:96203/1‑251 (MQ=255)
                                                                                                                                                                                                              cAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGccc                               >  1:130799/1‑251 (MQ=255)
                                                                                                                                                                                                                gaggtattCTGTTTGGTTTGCGTCGTGACCGAATGTCTCGTCCCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATGTCGGCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGGGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCAGACGCCTGAGGCTTGATCGCCAGCCGCCAGTAAACGTTGATAAACAATCCCGCGTTATAGCGGCGGCCCAg                             <  2:291031/244‑1 (MQ=255)
                                                                                                                                                                                                                          gtttggtttgCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAt                       <  1:76824/247‑1 (MQ=255)
                                                                                                                                                                                                                          gtttggttggCGTCGGGAACACATGTCTCAGCGCGATATAAGGGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAtt                      <  2:251053/248‑1 (MQ=255)
                                                                                                                                                                                                                                  tgCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGAc           >  1:281889/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATCGCGGCGCCCCAGCGCAATTGCCGCCACGCCATCg       >  2:195743/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aaCACATGTCTCAGCGCGCTATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGc                                                                                                                                                  >  2:27259/1‑107 (MQ=255)
                                                                                                                                                                                                                                           aaCACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGc                                                                                                                                                  <  1:27259/107‑1 (MQ=255)
                                                                                                                                                                                                                                           aaCACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTc  <  1:58982/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
ATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTC  >  NC_000913/4193575‑4194060

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: