Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,426,929 C→T C165Y (TGC→TAC)  yjfZ ← DUF2686 domain‑containing protein YjfZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,426,9290CT97.8% 145.9 / NA 46C165Y (TGC→TAC) yjfZDUF2686 domain‑containing protein YjfZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/23);  minor base G (1/0);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGA  >  NC_000913/4426689‑4427170
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
ggCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGTGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGAGGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCAt                                                                                                                                                                                                                                         <  2:243196/251‑1 (MQ=255)
     tagatagaTACGAGCGCACCTGATGACTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTggg                                                                                                                                                                                                                                    >  2:102640/1‑251 (MQ=255)
              ccgatcgcAACTGATGTCGCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAAGACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAgg                                                                                                                                                                                                                            <  1:386117/245‑1 (MQ=255)
               cGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAAGTATGCATTTGGGCAGGAAGGGt                                                                                                                                                                                                                          <  1:141242/251‑1 (MQ=255)
                   cgcAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTAc                                                                                                                                                                                                                       >  1:205775/1‑250 (MQ=255)
                     cAACTGATGTCTCGCTCCTGCGGCACTCGCTCGCAACCCAACCCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCCGACCAATTCTCAAATCATGGGGTAAATCAGAATCCGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGACGTTGGTCGCTATAACGATCGATTGTTGGTAGTTCAGATGACCAGGAATTAGGCATTTGGGCCGGAAggagatcggg                                                                                                                                                                                                                    >  2:386117/1‑243 (MQ=255)
                     cAACTGATGTCTCGCTCCTGCGGCACTCGCTCGCAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATACGGAACGAAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATCAGGCATTTGGGCAGGAAGGGTTACAAc                                                                                                                                                                                                                    >  2:88089/1‑251 (MQ=255)
                             gTCTCGCTCCGGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGAGCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCt                                                                                                                                                                                                             <  2:534737/250‑1 (MQ=255)
                                   cTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGttt                                                                                                                                                                                                      >  1:111944/1‑251 (MQ=255)
                                    tCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCACTTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTg                                                                                                                                                                                                     >  1:264837/1‑251 (MQ=255)
                                          ggCACTCGCTCGAAACCCACACCCGAGTAGCTTGCGGCTGCACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGACTCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGt                                                                                                                                                                                               >  2:214402/1‑251 (MQ=255)
                                                 gctcgAAGCCCAGACCCGCGTAGCTTGCGGCTGCTCTTGGTAAGGCACGCAGACAAATTTTACCATCAGGGGGTAAATCAGAGTACGCGTTTAGCAATGCAGCAATAATGAATCCGGGAGAGTTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGAAATTAGGCATTTGGGCGGGGAGGGTTCGAACGTATccagatcggaagagcgtggttagg                                                                                                                                                                                        >  2:76662/1‑229 (MQ=255)
                                                 gctcgAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTATCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCGGTACCAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGTTAGTCCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTGACAACGTATcc                                                                                                                                                                                                              <  1:76662/229‑1 (MQ=255)
                                                  ctcgAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTg                                                                                                                                                                                       >  2:323230/1‑251 (MQ=255)
                                                               cccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGTTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTTGAGTACCATAAcc                                                                                                                                                                          <  2:565335/251‑1 (MQ=255)
                                                                 cGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAg                                                                                                                                                                        <  1:193642/251‑1 (MQ=255)
                                                                             ggCTGAACTAGAGAAGGAACGCAGACCAATGTTCAAATCATGGGGTAAATCAGAATACGCGTTGAGCAATGCAGCAGTGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGAGCGATTGTTGGGATTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTAGCCTCTCGTCTGGCCGGTAGAGTCTGGAGTACCATAACCACGTCCCAAAAccg                                                                                                                                                            <  2:122085/251‑2 (MQ=255)
                                                                                       gATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgta                                                                                                                                                   <  2:111944/250‑1 (MQ=255)
                                                                                              aaCGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtattgtgtat                                                                                                                                           >  1:428223/1‑251 (MQ=255)
                                                                                                      ccAATTTTCAAATCATGGGCTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCAATAGAACAGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGCTCAGATGCCCAGTAATTAGGCATTTGGGCCGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGtt                                                                                                                                   >  2:568271/1‑251 (MQ=255)
                                                                                                           tttCAAATCATGGGGAAAATCATAATACGCGTTTAGCAATGCCGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAATGCGAAGTTGGTCGATATAAAGATAGATTGTTGGTAGATCAGATGACCAGTAATTAGGAATTTGGGCAGGAAGGGTTACAACGTATCCTATCATTTGGCCGGCATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTTTATTGTGTATTTACTGTTGGGtt                                                                                                                              >  2:107180/1‑251 (MQ=255)
                                                                                                             tCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAACCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTGGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGTAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCGCATCGTGTATTGTGTATTTACTGTTGGGTTc                                                                                                                             <  2:586983/250‑1 (MQ=255)
                                                                                                               aactCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGttt                                                                                                                                                                                                      >  2:306916/4‑175 (MQ=255)
                                                                                                               aaaTCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGttt                                                                                                                                                                                                      <  1:306916/175‑1 (MQ=255)
                                                                                                                    aTGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGAtt                                                                                                                     >  2:453102/1‑251 (MQ=255)
                                                                                                                               cAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGc                                                                                                          >  1:340679/1‑251 (MQ=255)
                                                                                                                                            tAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTc                                                                                              <  2:315045/250‑1 (MQ=255)
                                                                                                                                                      gcaATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGt                                                                                   <  1:568271/251‑1 (MQ=255)
                                                                                                                                                           gatgaATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATTAATATTTCAGTTCTGCGGTTAAc                                                                               <  1:430193/250‑1 (MQ=255)
                                                                                                                                                           gatgaATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAAcc                                                                              >  2:164855/1‑251 (MQ=255)
                                                                                                                                                               aaTCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGc                                                                          <  1:504680/251‑1 (MQ=255)
                                                                                                                                                                 tCCGAGGGCATTAGAAACGCCTGTAACAGTGCGACGATGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATCAATATTTCTGTTCTGCGGTTAACCTCGCAc                                                                        >  2:3934/1‑251 (MQ=255)
                                                                                                                                                                               aaaCGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTa                                                          <  1:140431/251‑1 (MQ=255)
                                                                                                                                                                                   gCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGt                                                      <  2:184267/251‑1 (MQ=255)
                                                                                                                                                                                     cTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTAt                                                    >  1:102914/1‑251 (MQ=255)
                                                                                                                                                                                         aaCAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTc                                                 <  1:287099/250‑1 (MQ=255)
                                                                                                                                                                                         aaCAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCt                                                <  1:291667/251‑1 (MQ=255)
                                                                                                                                                                                                            atatAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTa                              >  2:504719/1‑250 (MQ=255)
                                                                                                                                                                                                                aaaGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGt                                          <  2:195624/234‑1 (MQ=255)
                                                                                                                                                                                                                aaaGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGt                                          >  1:195624/1‑234 (MQ=255)
                                                                                                                                                                                                                        aTTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATg                 >  1:86350/1‑251 (MQ=255)
                                                                                                                                                                                                                          tgttgGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATgg                <  1:214402/250‑1 (MQ=255)
                                                                                                                                                                                                                             tgGTAGTTCAGATAACAAATACTTAGGCATTAGGGCAGGCATGGTTAAAACGTATCCTATAGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTACCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAAGATTCGACATCTTGCAATAAATTTTTCAGTTCTGTGGTTAACCTCGCACCCGCGTAAGATTTATTGTATTTCTGGGCGTGACTATTATCTCATGTCGCTTCATGGTtcc            >  2:341353/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAAAGTATCCGATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTACAAAAACCACTTCGTGTAGTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGGTCTGCGGTTAACCTCGCACCCGCGTAAGCGTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCAAGGTtccatcca       <  2:96308/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTAGTCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGGCCAAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTTGAAGGATTCGACATCTTGCAATACATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAATCTTTATTGTATTTCTTAGCGTGAATATTATCTAATGTAGTTTAATGGTTCCATCCaaaa    <  2:583216/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTCTCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGa  >  2:370181/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
GGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGA  >  NC_000913/4426689‑4427170

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