Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,489,775 C→T V97M (GTG→ATG)  yjgR ← DUF853 domain‑containing protein YjgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,489,7750CT98.2% 178.9 / NA 57V97M (GTG→ATG) yjgRDUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/20);  minor base G (0/1);  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

ATGCGCTACTGATATTACCGTACTGATTCTGGAAGGATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCC  >  NC_000913/4489529‑4490007
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
aTGCGCTACTGATATTACCGTACTGATTCTGGAAGGATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCCcatcacc                                                                                                                                                                                                                                      >  2:455501/1‑251 (MQ=255)
         tGATATTACCGTACTGATTCTGGAAGGATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTTAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGAtt                                                                                                                                                                                                                              <  2:537279/250‑1 (MQ=255)
          gATATTCCCGTCCTCATTCTGGAAGGATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAACGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGACGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCCACAGTCGTGCCAGCAACCGCGGCCCCAGATCCGCAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAACGCTATCCCACATCACCACCGGATTgg                                                                                                                                                                                                                            >  2:367207/1‑251 (MQ=255)
                  cGTACTGATTCTGGAAGGATTTGGCGTTATCGCCGATGTACTGGGTAATTGCACGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGt                                                                                                                                                                                                                    >  1:373150/1‑251 (MQ=255)
                   gTACTGGTGCTGGAAGGATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAAAAGCGGCCCCAGATCCGAAACCGTCACCCGCACCGGATGGCCTTTATCGCCAAAGATATCCCACATCACCACCGGATTGCCATGCGGtt                                                                                                                                                                                                                   >  2:7933/1‑251 (MQ=255)
                               gAAGGATTTGGCGTTATCGCCGCTGTACTGGGAAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTg                                                                                                                                                                                                        >  1:107517/1‑250 (MQ=255)
                                 aGGATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGCAATGGCATGCGGCTGCACGTCCTTGACg                                                                                                                                                                                                     >  1:378022/1‑251 (MQ=255)
                                  ggATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGc                                                                                                                                                                                                    >  1:537112/1‑251 (MQ=255)
                                   gATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGAGTGCACAGCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGcc                                                                                                                                                                                                   <  2:63747/251‑1 (MQ=255)
                                         gCGTTATCGCCGATGTACTTGTCTATTTCCCGCAGATCTTTAAAGTCGAGCAGTAACATCCCCTGATCGTCAGCAATGCGGAGGATGATATTCAGCACGCCGGATTGCACATCGTTGAGACTCAACAGTCGTGCCAGCACCAGCGGCCCCAAATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCTCCACCGGATAGGCAGGCGGTTGCCAGTCAGTGACACCGATAtt                                                                                                                                                                                             <  2:125285/251‑1 (MQ=255)
                                            ttATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATAttttt                                                                                                                                                                                          >  1:401111/1‑251 (MQ=255)
                                                   cGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCtt                                                                                                                                                                                   >  2:298142/1‑251 (MQ=255)
                                                          cTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGc                                                                                                                                                                             >  1:370534/1‑250 (MQ=255)
                                                           tGGGTAATTGACCGCAGATCTTTAAATTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCAAAATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGCAACCGACGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCCCCGGATTGGAATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCgcgc                                                                                                                                                                             >  2:309532/1‑246 (MQ=255)
                                                                               gttaaaggtgAGCAGTAACAGCCCCTGAGCGGCATCAACGCGGAAGATGATAGTCAGCACGCCAGAGTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATAACCCACATCACCGCCGGATTGGCAGGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCATGTTTTCCGACACCGTTcct                                                                                                                                                        <  2:36253/241‑1 (MQ=255)
                                                                                        gAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGcgcg                                                                                                                                               >  1:276493/1‑250 (MQ=255)
                                                                                        gAGCAGTAACAGCACCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGcgcgc                                                                                                                                              >  1:103366/1‑250 (MQ=255)
                                                                                                 cAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCa                                                                                                                                     <  2:233941/251‑1 (MQ=255)
                                                                                                    ccccTGATCGTCAGCAATGCGGCAGATGATATTCCGCACGCCAGATAGCACATCGTTGCGATTCCACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCACGCACCGGATGGCCTTTCTCGACAACGATATCCCACATCACCACCGGATTGGCATGCGGGTGCCAGTCATTGCCGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGCTCCTGCCTGCGCGCCGCCGGTc                                                                                                                                      >  2:265302/1‑247 (MQ=255)
                                                                                                         gATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCt                                                                                                                             >  2:76590/1‑251 (MQ=255)
                                                                                                            cGTCAGCAAGGCGTAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGGTCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCCCCGGATGGCCTTTCTCTCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGGCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGGTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTc                                                                                                                          <  2:537112/251‑1 (MQ=255)
                                                                                                                  caatcaagAAGATGATATTCAGCACGCCAGATATCACATCGTTGAGAATCAACAGCCGTGCCAGCAACAGCGGCCCCAGGTCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCa                                                                                                                    <  2:164559/244‑1 (MQ=255)
                                                                                                                     tGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCAGGTCAGATCGCCTTTCACATCAGcc                                                                                                                 >  1:483770/1‑251 (MQ=255)
                                                                                                                         gAAGATGATATTCAGCACGCCAGATTTCACATTGTTGAGATTCAACAGTCTTGCCAGCAACAGCTGCCCCAGATCCGAACCCGGCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCGGCGCGACGCCGGTCAGTTCTCCTTTCACATCAGCCATaa                                                                                                             <  2:53046/251‑1 (MQ=255)
                                                                                                                         gAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATaa                                                                                                             <  1:481351/251‑1 (MQ=255)
                                                                                                                             atgatATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAacac                                                                                                         >  2:368087/1‑251 (MQ=255)
                                                                                                                             atgatATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGAGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAacac                                                                                                         >  1:279318/1‑251 (MQ=255)
                                                                                                                                    tCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACg                                                                                                  >  1:417849/1‑251 (MQ=255)
                                                                                                                                                 ttGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGccc                                                                                                   >  2:428347/1‑235 (MQ=255)
                                                                                                                                                 ttGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCAc                                                                                                   <  1:428347/237‑1 (MQ=255)
                                                                                                                                                    cacaTCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCATGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATg                                                                                  >  1:306575/1‑251 (MQ=255)
                                                                                                                                                      caTCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGAc                                                                                >  2:175317/1‑251 (MQ=255)
                                                                                                                                                         cggtGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGGCGCCCGCACCGGATGGCCTTTCCCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTcc                                                                             <  1:455501/248‑1 (MQ=255)
                                                                                                                                                         cGTTGTGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTcc                                                                             >  2:122383/1‑251 (MQ=255)
                                                                                                                                                         cGTTGATATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATTGCCTTTCTCGCCAAAGATAGCCCACATCACCACCGGATTTGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTTCGACAATGACTcc                                                                             <  1:106469/251‑1 (MQ=255)
                                                                                                                                                          gTTGAGATTCCACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCg                                                                            >  1:186422/1‑251 (MQ=255)
                                                                                                                                                           ttGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCCTAAACACCGGCACGCCGATTTCCGACAATGACTCCGc                                                                           >  1:517617/1‑251 (MQ=255)
                                                                                                                                                           ttGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGAATGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTACGACAATGACTCCGc                                                                           >  1:525756/1‑251 (MQ=255)
                                                                                                                                                            ggagaTTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCGCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTGAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTGTCACATCAGCCATAAACACCGGCACGCCGATGTCCGTCAATGACTCCGcc                                                                          <  2:276493/250‑1 (MQ=255)
                                                                                                                                                                      cAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGc                                                                 >  1:64805/1‑250 (MQ=255)
                                                                                                                                                                            ggcccgCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGTATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTGTCCGACACCGTCCCTTCCCGCGCGACGCCGTTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCATTTTTTGCAGCGTaa                                                          <  2:269072/246‑1 (MQ=255)
                                                                                                                                                                                  cAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGttt                                                    <  2:339451/251‑1 (MQ=255)
                                                                                                                                                                                   aaCAGCGGCCCCAGATCCGACACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCTACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGtttt                                                   >  1:188334/1‑251 (MQ=255)
                                                                                                                                                                                   aaCAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCCCCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTGACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGtttt                                                   <  1:76590/251‑1 (MQ=255)
                                                                                                                                                                                    aCAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGtttt                                                   >  1:503127/1‑250 (MQ=255)
                                                                                                                                                                                       gCGGCCCAAGATCCGAACCCGTCGCCCGCACCGGTTGGCCTTTCTCGCCAAAGATATCCCACAGCCCCACCGGATTGGCATGCGGGTGCCAGTCATGGACGCCGATATTTTGAAGCCTGGCGAGCAGTTTTTCCGACACCGTTCCTGCCGGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGAGTCCGCCAGTTTTTGCAGCGGAACGGTTTTAccc                                               <  2:567802/251‑1 (MQ=255)
                                                                                                                                                                                                   ccGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTgc                                    <  1:440052/250‑1 (MQ=255)
                                                                                                                                                                                                         ccGTCGCTCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCCCCGGATTGGCATGCGGTTGCCAGTCATTGCCGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAACCACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTTCGCCGGt                              >  2:95588/1‑250 (MQ=255)
                                                                                                                                                                                                         ccGTCGCTCGCACCGGAAGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGt                              >  2:204217/1‑250 (MQ=255)
                                                                                                                                                                                                                  gCACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGcccc                    >  2:52582/1‑250 (MQ=255)
                                                                                                                                                                                                                   cACCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCAt                   >  2:314697/1‑251 (MQ=255)
                                                                                                                                                                                                                    aCCGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATg                  <  1:121473/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGGATGGCCTTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGAc                <  2:370534/251‑1 (MQ=255)
                                                                                                                                                                                                                         aTGGCATTTCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCAGTAATCAGCCCATGACGGt             >  2:293672/1‑251 (MQ=255)
                                                                                                                                                                                                                               tttCTCGCCACAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGc                                                                                                                  >  2:397733/1‑144 (MQ=255)
                                                                                                                                                                                                                               tttCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTACGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGcc        >  1:86495/1‑250 (MQ=255)
                                                                                                                                                                                                                               tttCTCGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGc                                                                                                                  <  1:397733/144‑1 (MQ=255)
                                                                                                                                                                                                                                    cGCCAAAGATATCCCACATCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTccc  >  1:149088/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
ATGCGCTACTGATATTACCGTACTGATTCTGGAAGGATTTGGCGTTATCGCCGATGTACTGGGTAATTGCCCGCAGATCTTTAAAGTCGAGCAGTAACAGCCCCTGATCGTCAGCAATGCGGAAGATGATATTCAGCACGCCAGATTGCACATCGTTGAGATTCAACAGTCGTGCCAGCAACAGCGGCCCCAGATCCGAAACCGTCGCCCGCACCGGATGGCCTTTCTCGCCAAAGATATCCCACACCACCACCGGATTGGCATGCGGTTGCCAGTCATTGACGCCGATATTTTTAAGCCTTGCGAGCAGTTTTTCCGACACCGTTCCTGCCTGCGCGACGCCGGTCAGATCGCCTTTCACATCAGCCATAAACACCGGCACGCCGATTTCCGACAATGACTCCGCCAGTTTTTGCAGCGTAACGGTTTTACCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCC  >  NC_000913/4489529‑4490007

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