Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,894,839 T→C L12P (CTC→CCC)  pabB → aminodeoxychorismate synthase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,894,8390TC97.4% 133.4 / ‑4.7 38coding (35/1362 nt)pabBaminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/16);  minor base . (0/1);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGTGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATT  >  NC_000913/1894594‑1895080
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
ggAGTGGTTTGCACTCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTTATTACTTTACCCTgtc                                                                                                                                                                                                                                               <  1:481252/251‑3 (MQ=255)
ggAGTGGTTGGCACGCTGCTCCTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTTATCCCTTGTGCCATCAGTTCCTGTATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTgga                                                                                                                                                                                                                                               <  1:481265/251‑2 (MQ=255)
  aatGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGCTTGCCTGGCCGCTGATCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGGTTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGt                                                                                                                                                                                                                                            >  2:247617/3‑251 (MQ=255)
        tGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGAc                                                                                                                                                                                                                                      >  1:508063/1‑251 (MQ=255)
          gCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGAcgc                                                                                                                                                                                                                                    >  1:314338/1‑251 (MQ=255)
            aCGCAGCTCTTCTGACACAAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTCTTTCATGCTACCCGGCACTTAGCCAGTAGAGTCAGGACTGGTGAAGACGTTATTTCCCGCTGTGATTACTTTACCCTGGCGTCAGGAcgccg                                                                                                                                                                                                                                  >  2:328812/1‑251 (MQ=255)
                                              ccGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGcc                                                                                                                                                                                                <  2:506858/251‑1 (MQ=255)
                                                  tgctCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTCAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGTCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTCTTTCTCCCGCTTAAGCCACCTGCCGTgg                                                                                                                                                                                            >  2:90739/1‑251 (MQ=255)
                                                   gctCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTAATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTggg                                                                                                                                                                                           >  1:567902/1‑251 (MQ=255)
                                                          ccTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCt                                                                                                                                                                                    <  1:474264/251‑1 (MQ=255)
                                                             tgtgCCATCATTTCCGGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCtt                                                                                                                                                                                   <  2:433867/249‑1 (MQ=255)
                                                                     cAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTcc                                                                                                                                                                          <  1:410017/250‑1 (MQ=255)
                                                                      aGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCgg                                                                                                                                                                        >  1:223431/1‑251 (MQ=255)
                                                                                 ccGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGatc                                                                                                                                                             <  1:241607/251‑1 (MQ=255)
                                                                                      cgcaaGCTTTTTGCTGCTGTTAGTGGGTGAGTGAAGGTAAACCTGCAAACATTGTTACCTCCTGCTAAATTGTGGGCGCTAATTATTTCACGCTACCCGGCACATAGCCAGTGGAGTCAGGACTGATGAAGACTTTTTCTTCCGCTTTGATTATTTTTACCTGGCGTCAGGAGGCTGCTGAATTTTATTTCTCCTGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATcc                                                                                                                                                         <  2:147694/246‑1 (MQ=255)
                                                                                        aCAGCTTTTTGCTGCTGTTCATGGGGGAGTTAAGGTAAACCTGCAAACATTGTGAACTCCTGCTAAAGTGTTGGCGCTAATTATTTCAGGCTACCCGGCCCATAGCCAGTAGAGACAGGACTTATGAAGACGGTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCAACCTGCCTTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGta                                                                                                                                                      <  2:559121/251‑1 (MQ=255)
                                                                                                            aTGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTgg                                                                                                                                  >  2:240231/1‑251 (MQ=255)
                                                                                                               gaagagagaAGGAAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCg                                                                                                                               >  2:116798/8‑251 (MQ=255)
                                                                                                                      gaAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCCCCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCgccgagcc                                                                                                                        <  1:152795/251‑1 (MQ=255)
                                                                                                                       aaGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGAGGCTTTTACACTCCGGCTATGCCGATCATCCGTATATCCGCTTTGATATTGGTGTCgccgagccg                                                                                                                       <  1:408044/251‑1 (MQ=255)
                                                                                                                       aaGGAAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGATAACTATCTCATGCTACCCGGCACATAGCCAGTAGATTCAGGCCTGATGAAGACGTTATCCCCCGCTGTGATTCCTTTCCCCTGGCGTCAGGACGCGGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGGGGGCGATGCTTTTACACTCCGGCTACGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTTgctgagccg                                                                                                                       <  2:518644/251‑1 (MQ=255)
                                                                                                                                   cAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACttt                                                                                                            >  1:508735/1‑250 (MQ=255)
                                                                                                                                              aaCTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCgg                                                                                                >  1:87765/1‑251 (MQ=255)
                                                                                                                                                   cTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCggccaac                                                                                           >  2:307641/1‑246 (MQ=255)
                                                                                                                                                    tGCTAAATTGTTGGCGCTAATTATTTCATACTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTCTTTCTCCCGCTTAAGCCACCTTCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTaa                                                                                             >  1:443449/1‑248 (MQ=255)
                                                                                                                                                     gCTAAATTGTTGGCGCTAATTATTGCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGa                                                                                          <  1:435999/250‑1 (MQ=255)
                                                                                                                                                             gttgGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCgttgt                                                                                 >  2:119915/1‑251 (MQ=255)
                                                                                                                                                                  cgtcaaTTATCTCATGCTCCCCGGCACATGGCCAGTAGAGTCAGGACTGATGAAGACTTTATCTCCCGTTGTGATTTCTTTGCCCGGGCGTCGGGACGCCGCTGAATTTTATTTCCCCCGCTTAAGCCACCGGCCTTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTAGAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTg                                                                            <  2:532699/247‑1 (MQ=255)
                                                                                                                                                                  cgcTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGTCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTg                                                                            >  1:114585/1‑251 (MQ=255)
                                                                                                                                                                      aaTTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAg                                                                        <  1:213994/251‑1 (MQ=255)
                                                                                                                                                                        ttattTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCg                                                                      >  2:98708/1‑251 (MQ=255)
                                                                                                                                                                                   tACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGGATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAATAAACCGTTGTTAGAGAAAGCGAAAAACGcacc                                                           >  2:232233/1‑250 (MQ=255)
                                                                                                                                                                                               aGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACtgat                                               >  1:12958/1‑251 (MQ=255)
                                                                                                                                                                                               aGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACtga                                                <  1:346238/250‑1 (MQ=255)
                                                                                                                                                                                               aGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGAGATTACTTTACCCTGGCGTAAGTACGCCGCTGAATTTTATTTCTCCCGCTTAAGAAACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCAGGATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGAAAAGACACCGTTGTTAGGGAAAGCGAAAAACGCACAACGACCACtga                                                >  2:353648/1‑250 (MQ=255)
                                                                                                                                                                                                              gACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTg                                >  2:193566/1‑251 (MQ=255)
                                                                                                                                                                                                                cgtatgaAGACTTTATCTCCCGCTGTATTTACTTTCCCTGGGCGTCAGGCCGCCGTGGAATTTTTTTTCCCCCTCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCt                              <  2:308336/248‑1 (MQ=255)
                                                                                                                                                                                                                           gTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGc                   >  1:269739/1‑251 (MQ=255)
                                                                                                                                                                                                                                             tactttaccCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACAtt  <  1:549171/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
GGAGTGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATT  >  NC_000913/1894594‑1895080

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