Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,563,567 T→C I13V (ATT→GTT)  yjiK ← PF06977 family protein YjiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,563,5670TC94.3% 152.6 / 0.1 54I13V (ATT→GTT) yjiKPF06977 family protein YjiK
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (32/18);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

GTGCGGATAAAATCACCATTGGTGGTCATTTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCC  >  NC_000913/4563323‑4563811
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gttcGGATAAAATCACCATTGGTGGTCATTTCAACAATGGCGGCAGGGTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGAt                                                                                                                                                                                                                                                <  2:661539/248‑1 (MQ=255)
gtgtggATAAAATAACCATTGGTGGGCATGTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGGGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGAGTGCACGAAGAAAGTACACCCTGCCACAATGGCATAAAGAATCACAATAACGGAGAt                                                                                                                                                                                                                                                <  2:374275/247‑1 (MQ=255)
    ggATAAAATCACCATTGGTGGTCATTTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGt                                                                                                                                                                                                                                            >  1:474012/1‑251 (MQ=255)
     gATAAATTCACCATTGTTGGTCATTTCAACAATGGCGGCAGGATTGTTGATAGTGCTAAAAAGAGTATTTCTTTTTGCCGACCAGGTCAGCGAAGAGATCTTATTCGTTCTGCCTGCAATTTCTTTGCCGACAATGGTGGCGTGGTCATTCTGAAATGCTGCTGCATGTTTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACCATGCCACAATGGCAAAAAGAATAACAATCACGGAGCTGCGtt                                                                                                                                                                                                                                           >  2:632789/1‑251 (MQ=255)
             cACCATTGGTGGTCATTTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAAc                                                                                                                                                                                                                                   >  1:263679/1‑251 (MQ=255)
                   tgggggtCATTTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACACTGGCAAACAGAATCACAATAACGGAGATACGTTAACTTAACGAAAt                                                                                                                                                                                                                              >  2:500993/5‑250 (MQ=255)
                             ttCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTc                                                                                                                                                                                                                    <  1:187159/250‑1 (MQ=255)
                                        cggcAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCt                                                                                                                                                                                                        >  1:89882/1‑251 (MQ=255)
                                          gcAGATTTGTTGATAGAGCTAAAAAGAGGATTAATTTGTGCCGCACAGATCAGAGTAGAGATATTAATCGTTATGCATGCAATTTCTTTGCCGTCAAAGGCGGCGTGGTAATTCTGAAATGATGCTGCATGATTGATTTTACGGGCGCAAGATTGACCGAAGAAAGTACAACTTGCCACAATGTCAAAAAGAATCACAATAAGGGAGATGCGATTGATTCACGAAGTACTTTTTGtaccggaactcccttt                                                                                                                                                                                                      >  2:35643/1‑236 (MQ=255)
                                                 tgttgCTAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGAGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCAAGGTAATCCCTTTAAGGCa                                                                                                                                                                                                >  2:437669/1‑250 (MQ=255)
                                                         gTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATa                                                                                                                                                                                        >  1:563075/1‑250 (MQ=255)
                                                            cTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATa                                                                                                                                                                                        >  1:327593/1‑247 (MQ=255)
                                                                    agTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGTTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGaa                                                                                                                                                                            <  1:95070/251‑1 (MQ=255)
                                                                       aTTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTACGGCAAATCAATAAGATATCGAGAAGt                                                                                                                                                                          >  1:100064/1‑250 (MQ=255)
                                                                              tgtgCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCACTGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACAGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTATAGGCAAATCAATAAGATATCGAGAGGTATACGt                                                                                                                                                                    >  2:6583/1‑249 (MQ=255)
                                                                               gtgCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGc                                                                                                                                                                 <  2:176255/251‑1 (MQ=255)
                                                                                   cGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGAcc                                                                                                                                                             <  1:618592/251‑1 (MQ=255)
                                                                                    gACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCt                                                                                                                                                            >  1:139913/1‑251 (MQ=255)
                                                                                    acatgacgctcctccgtcctagtattCGTTTTTCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGTTAATTCTTAAATGATGCTGCAATATTGCTTTTACGGGCGCAAGATTGTACGAAGAAAGTACTCACTGCCACAATGGCAAAAATAATCACAATAACGGAGATGCGTTTACTGAACGAAATACTTTTTGGCACGGTAAGCCCTTTAAGGCAAATCACTAAGATTTCGAGAAGTAGACGTAGCGAcc                                                                                                                                                             <  2:675898/228‑1 (MQ=255)
                                                                                                      atatTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATGCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATAGCGAGAAGTATACGTAGCGACCTTAATGAATCCGTAAcgcg                                                                                                                                          <  2:611270/251‑1 (MQ=255)
                                                                                                        attattCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGAcc                                                                                                                                                             >  1:675898/1‑230 (MQ=255)
                                                                                                                ttATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTTCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAg                                                                                                                                                                  <  2:496083/217‑1 (MQ=255)
                                                                                                                ttATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAg                                                                                                                                                                  >  1:496083/1‑217 (MQ=255)
                                                                                                                               cTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGc                                                                                                                  >  2:214505/1‑250 (MQ=255)
                                                                                                                               cTTTGCCGACAATGGTGGCGTGGTAATTCTGAAATGCTGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTAACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGAACTTAATGAATCATTAACGCGGAGAAACCCCGTAATGTGCAATGcc                                                                                                                 >  2:534431/1‑250 (MQ=255)
                                                                                                                                    ccGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCCAGCTTGCACGAAGAAAGTAAACACTGCCACAATGGCAAACAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCACATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGcc                                                                                                            >  1:82675/1‑251 (MQ=255)
                                                                                                                                              gtgCGTGGTAAATCAGAAATGATGCTGCATGATTGCTTTTAAGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGAAAAAATAATCCCAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTGTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAAAGAAGCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACTCCAta                                                                                                  <  2:364881/249‑2 (MQ=255)
                                                                                                                                                  gTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGc                                                                                                                       <  1:184346/226‑1 (MQ=255)
                                                                                                                                                  gTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGc                                                                                                                       >  2:184346/1‑226 (MQ=255)
                                                                                                                                                             tacatgatgCTGCATGATTGCTTTAAAGGGCGCAAGCTTGCAAGAAGAAAGTACACACTGCCACAATGGCAAACAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTAACGGTAATCCATTTAAGGCAAATCAATAAGATATCGAGAAGAATACGTAGCGACATTAATGAAACCTTAACGCGGAGAAACACCGTAATGTGCAATGCAGCGCCAGCCACCATCCTGATGGaaataaat                                                                                   >  2:52476/4‑251 (MQ=255)
                                                                                                                                                             gAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAAt                                                                                   <  1:214505/251‑1 (MQ=255)
                                                                                                                                                               aaTGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGAACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTg                                                                                 >  1:232272/1‑251 (MQ=255)
                                                                                                                                                                     gctgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACtcag                                                                           <  2:319825/251‑1 (MQ=255)
                                                                                                                                                                         cATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAAcgc                                                                                                                                           <  1:551023/183‑1 (MQ=255)
                                                                                                                                                                         cATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAAcgc                                                                                                                                           >  2:551023/1‑183 (MQ=255)
                                                                                                                                                                         cATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACtcagtc                                                                         >  1:361618/1‑249 (MQ=255)
                                                                                                                                                                         cATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACtcagtc                                                                         >  1:630323/1‑249 (MQ=255)
                                                                                                                                                                                        gggCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCAccc                                                         >  1:149801/1‑250 (MQ=255)
                                                                                                                                                                                               aGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTAtt                                                 >  1:436781/1‑251 (MQ=255)
                                                                                                                                                                                                   tGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCt                                             >  1:198617/1‑251 (MQ=255)
                                                                                                                                                                                                     cACGAAGCAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTAACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCg                                           >  2:59523/1‑251 (MQ=255)
                                                                                                                                                                                                     cACGAAGAAAGGACACACTGCCACAATGGCAAAAAGAATCACCATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATCCGTAGCGACCTTAATGAATCCTTAACGCGGAGCAACCCCGTAATGTGCAATGCAGCGCCAGACACAATCCTGCTGGCCATAAATTGAATCCATCAGCCCAGGGATCACCCATTGTATTTCCTCg                                           >  2:406/1‑251 (MQ=255)
                                                                                                                                                                                                       cGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGCCACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGcc                                         >  2:598463/1‑251 (MQ=255)
                                                                                                                                                                                                       cGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGcc                                         <  1:437669/251‑1 (MQ=255)
                                                                                                                                                                                                        gaagaaAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGAAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGCGAAACCCCATAATGTGCAATGCAGCGCCAGACACCATCCTGCAGGCAAAAAATTGACTAAGGCAGCCAAGGAAGCACCCATTGTCTTTCCTCCCCg                                        >  2:366875/1‑251 (MQ=255)
                                                                                                                                                                                                                    cacTGCCACGATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTATCGAAATACTTTTAGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCATTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAACTTGACTCAGTCAGCCCAGGGATCACCCAATGTATTTCATCGCCGGGGGGACAGGTc                            >  2:462843/1‑251 (MQ=255)
                                                                                                                                                                                                                     acTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCa                           <  1:115817/251‑1 (MQ=255)
                                                                                                                                                                                                                       tGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAAt                                                                                                                    <  2:52363/160‑1 (MQ=255)
                                                                                                                                                                                                                       tGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAAt                                                                                                                    >  1:52363/1‑160 (MQ=255)
                                                                                                                                                                                                                          ccACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGACAGcc                        >  2:606383/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAACTCAATAAGATATCGAGAAGTATAAGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATAACCCATTGTATTTCCTCGCCGGGGTGACAGGTCAGCCAAATGCATACCCCGc         >  2:443498/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAgcc     >  1:614328/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cacaATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAgccg    <  1:675034/251‑1 (MQ=255)
                                                                                                                                                                                                                                               caATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAgccgcc  <  1:60743/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 ataacGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAg                               >  1:323994/1‑220 (MQ=255)
                                                                                                                                                                                                                                                 ataacGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAAGAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAg                               <  2:323994/220‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GTGCGGATAAAATCACCATTGGTGGTCATTTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCC  >  NC_000913/4563323‑4563811

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