Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,529,039 C→T G402D (GGT→GAT)  sgcC ← putative PTS enzyme IIC component SgcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,529,0390CT98.2% 185.9 / NA 57G402D (GGT→GAT) sgcCputative PTS enzyme IIC component SgcC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (36/20);  total (36/21)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGT  >  NC_000913/4528797‑4529284
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
ggTCGATAACCCAGTTCACCCCCAACTGCGCACCAAAGCCCGGGGCGCCGGGCAATGCACGGAGATGGCACATGGCGATGACCGCTTTTTCCGCCCCAATAACCTCTTTCAGCCACCCCATTGAAAGTCCTCCTTATCCCGATAATTCCCTTTTATGTAAGAAAAGCACACCGCTGGCAACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACat                                                                                                                                                                                                                                                <  2:88640/249‑1 (MQ=255)
    gATAACCCAGTTCATCCCCAACTGCGCTTCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGc                                                                                                                                                                                                                                           >  1:69220/1‑251 (MQ=255)
         cccAGTTCATCCACACCTGCGCATAAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCACTAACCTATTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACCATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGCCCAATCGAACATATTGCCCACg                                                                                                                                                                                                                                      >  2:332495/1‑251 (MQ=255)
           cAGTTCACCCCCAACGGCAAATCAAAACTCGGGACGCCGGTCAATGCCCGTAGAGGGCACATGGCGAGGACCGCTTTCCCCGCCCCAATAACCTCTTTCAGCAAACTCATGGAAAGTACTCCTGACACCGAAAATCCACTTTTCCGTAAGAACAGCACGACGCGGGCTACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACAGATTGCCCACGGa                                                                                                                                                                                                                                    <  2:92585/251‑1 (MQ=255)
                      cAACTGCGAATCAAACCTGGGGTCGCCGGCCAAGGCACGTAGATGGCACATGGCGGTGACCGTTGTTTCCGCACCAATACCCTCTTTCAGCCAACGAATGGAAAGTACTCCTTATTCCGATAATCCCCGTTTACGTAAGACCAGCACGACGCTGGCGACAATTCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTgaga                                                                                                                                                                                                                         <  2:364848/251‑1 (MQ=255)
                       aaCTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGAt                                                                                                                                                                                                                        <  1:6127/251‑1 (MQ=255)
                        aCTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAATGTACTCCTTATTCCGATAATTCCCTTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGAtt                                                                                                                                                                                                                       <  1:9563/251‑1 (MQ=255)
                                        cGGGTCGCAGGGCCATGCACGAAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACACCGCCAATAATCCCAAGCGACATCAATTCAGAAATCGCCCAATCGAACATATTGCCAACGGACAACGCTGCGATTTGCGCACTTTCACCgg                                                                                                                                                                                                       >  2:320495/1‑251 (MQ=255)
                                           gTCGCCGGGCAAAACCCGCAGAGGGCCCATGGCGTTGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACATAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGATATTTGCGCACTTTCACCgcca                                                                                                                                                                                                     >  2:314215/1‑247 (MQ=255)
                                                  ggCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCACCTCATTGAAAGTCCTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCGCCGGCAAAACTGaa                                                                                                                                                                                             <  1:84032/251‑1 (MQ=255)
                                                           cGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAACACTGAAACCGCCttt                                                                                                                                                                                    >  1:138583/1‑251 (MQ=255)
                                                                 tGGCACCTGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCACACTTTCACCGGCAAAACTGAAACCGCCTTTGAgcgc                                                                                                                                                                              >  1:357060/1‑251 (MQ=255)
                                                                                                cAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCtt                                                                                                                                                                                     >  1:258393/1‑213 (MQ=255)
                                                                                                aaaaacccTCTTCCAGCAAACCCATAGAAAGTACTCCTTATCCCAATAATTCACGTTTACGAAAAAAAAGCACGACGCGGGCAACAAACCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACTCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCtt                                                                                                                                                                                     <  2:258393/207‑1 (MQ=255)
                                                                                                      ccTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATc                                                                                                                                         >  1:159111/1‑251 (MQ=255)
                                                                                                      ccTCTTTCAGCCAACGCATGGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAAAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATc                                                                                                                                         <  1:7153/251‑1 (MQ=255)
                                                                                                         ctTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATcagc                                                                                                                                      >  1:334054/1‑251 (MQ=255)
                                                                                                               gCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGAAGGACg                                                                                                                                >  1:76740/1‑251 (MQ=255)
                                                                                                                    cTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAat                                                                                                                           >  2:127754/1‑251 (MQ=255)
                                                                                                                        ttGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAatcatc                                                                                                                       >  1:157641/1‑251 (MQ=255)
                                                                                                                              gTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATg                                                                                                                 >  1:430225/1‑251 (MQ=255)
                                                                                                                              gTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATg                                                                                                                 >  1:452222/1‑251 (MQ=255)
                                                                                                                                 cTCCTTATTCCGCTAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGAcc                                                                                                              >  2:54944/1‑251 (MQ=255)
                                                                                                                                          ccGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTAAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATc                                                                                                     >  1:315242/1‑251 (MQ=255)
                                                                                                                                                aTTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGACGTAAGGCGCGAACTGGGTGGCGATCAGCAGGAAGGTAATCATCACGATGACCCCGCTTATCAACGTc                                                                                               >  2:18928/1‑251 (MQ=255)
                                                                                                                                                         tACGTAAGAACAGCCCGCCGCTGGCAACAAACCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAg                                                                                       <  1:438762/250‑1 (MQ=255)
                                                                                                                                                             tAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGtctc                                                                                  >  1:8621/1‑251 (MQ=255)
                                                                                                                                                             tAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGAAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACAATGACCCCGCTTATCAACGTCCTGACCAGGtctc                                                                                  >  1:369247/1‑251 (MQ=255)
                                                                                                                                                                 aaCAGCACGACGCTGACGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGACATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGCCGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCcgcg                                                                              >  2:183502/1‑251 (MQ=255)
                                                                                                                                                                   cAGCACGACGCTGGCGACACTACCGACCGCCACGACAACGCCAATAATCCCACGCGACATCAATTCAGCAATCGCCCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCCCGATGACCCCGCTTATCAACGTCCAGACCAGATCTCCGCggg                                                                            >  2:340609/1‑250 (MQ=255)
                                                                                                                                                                    aGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTg                                                                           >  2:276022/1‑251 (MQ=255)
                                                                                                                                                                      cACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTgg                                                                          >  1:228540/1‑250 (MQ=255)
                                                                                                                                                                         gacgCTGGCGACAATACCGCCCGACACGACAACGCCAATAATCCCACGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATa                                                                       >  2:67138/1‑250 (MQ=255)
                                                                                                                                                                         gacgCCGGCGACAATACCGACCGCAACGACACCGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGCACATATTGCCCCCGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTACGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATa                                                                       >  2:448710/1‑250 (MQ=255)
                                                                                                                                                                          ccgtagcccacaaaccccaccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCAAACATATTGCCCACGGACACCGCTGAGATTTGCGCACTTTCCCCGGCAAAACTGAAACCGCCTGTGAGCGCCAGTTCAGTGAAGGAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACTATGACCCCGCTGATCAACGTCCTGACCAGGTCTCCGCGGTGGATAAc                                                                     <  2:398180/233‑1 (MQ=255)
                                                                                                                                                                           cgCTGGCGACAATCCCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAAcc                                                                    <  1:18928/251‑1 (MQ=255)
                                                                                                                                                                               ggCGAAAAAACCGACCCCCACGACAACGCCAAAAACCCCAAGCGACATGAATTCAGAAATCGACCAATCGAACAGATTGCCCACGTACATCGGGGAGATTTGCGCGCTTTCCCCGGCAAAACTGAAACCGCCTTTGAGCGCCATGTCATTGAAGTAAGGCGCGAGCTGTGTGTCGATCGGCGGGTCGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAg                                                                <  1:434801/251‑1 (MQ=255)
                                                                                                                                                                                gagacaaccccgcccgCCACGACAACGCCAATAAACCCAAGCGACAACAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGTTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCCATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGGCGCCGCGGTAGATATCCGTAGc                                                               <  2:315242/238‑1 (MQ=255)
                                                                                                                                                                                  gACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGcc                                                              >  1:92274/1‑250 (MQ=255)
                                                                                                                                                                                     aaTACCGAACGCCACGACAACGCCCATAATCCCAAGCGACATCAATTCCGAAATCGACCAATAGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGCAACCGCCTTTGAGAGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGCCGGCAATCATCCCGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGAAGCCATAc                                                          >  2:397374/1‑251 (MQ=255)
                                                                                                                                                                                      accccacccgCCACCAAAACGCCAATAACCCCAAGCGACATCAATTCAGAAATCGACCAACCAAACAAATGGCCCCCGGACAACGTTGAGATTTGCTCACTTTCTCCGTCAAAACTGAAACCGCCGTTGAGCGCCTTTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCTGCATTACGGTAATCATCACGATGACCCCGCTTAGCAACGTCCTGACCAGGTCGCCGCTGGGGAGAACCGTAGTCATACa                                                         <  2:413337/244‑1 (MQ=255)
                                                                                                                                                                                       tACCGACCGCCACGACAACGACAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATAc                                                          >  1:369508/1‑249 (MQ=255)
                                                                                                                                                                                           gaccgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATaa                                                    <  2:213417/251‑1 (MQ=255)
                                                                                                                                                                                            accgCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATaaa                                                   <  1:276022/251‑1 (MQ=255)
                                                                                                                                                                                                    gACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGACCGGTg                                           >  1:2906/1‑251 (MQ=255)
                                                                                                                                                                                                               taatCCCAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGCCCAGGTCTCCGCGGTGGATAACCGTAGCCATCCAGATAAAGAACGGTGCTACCGGCAgg                                >  2:170478/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGc                           <  1:157314/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccAAGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGGCTgg                           >  1:386532/1‑250 (MQ=255)
                                                                                                                                                                                                                      aaGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGGCTGCgg                         >  2:335516/1‑250 (MQ=255)
                                                                                                                                                                                                                       aGCGACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGCTAAAGAACGGTGCTACCGGCAgggttgcggg                        >  2:291707/1‑243 (MQ=255)
                                                                                                                                                                                                                          gACATCAATTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAgggg                     >  1:368902/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAATACAGAAAGCGACCAATCGAACAGATGGCCCACGGACACCGCTGAGATTGGCGCACTCTCACCGTCAAAACAGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTACGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAac                 <  1:314215/250‑1 (MQ=255)
                                                                                                                                                                                                                                aaTTCAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAgggggaaacc               >  2:147673/1‑245 (MQ=255)
                                                                                                                                                                                                                                    cAGAAATCGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACttt            <  1:30887/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cGACCAATCGAACATATTGCCCGCGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCg    <  1:211053/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cGACCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCg    <  1:248408/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aCCAATCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGCGTGTGGGGGGCAACACTTTGTTGCCCggg  >  2:137070/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
TGTCGATAACCCAGTTCATCCCCAACTGCGCATCAAAACTCGGGTCGCCGGGCAATGCACGTAGATGGCACATGGCGATGACCGCTTTTTCCGTCCCAATAACCTCTTTCAGCCAACTCATTGAAAGTACTCCTTATTCCGATAATTCACGTTTACGTAAGAACAGCACGACGCTGGCGACAATACCGACCGCCACGACAACGCCAATAATCCCAAGCGACATCAATTCAGAAATCGACCAACCGAACATATTGCCCACGGACAACGCTGAGATTTGCGCACTTTCACCGGCAAAACTGAAACCGCCTTTGAGCGCCATTTCAGTGAAGTAAGGCGCGAACTGGGTGGCGATCAGCAGGACGGTAATCATCACGATGACCCCGCTTATCAACGTCCTGACCAGGTCTCCGCGGTGGATAACCGTAGCCATACAGATAAAGAACGGTGCTACCGGCAGGTCTGCGAGGGGTAACACTTTGTTGCCCGGT  >  NC_000913/4528797‑4529284

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