Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,555,060 A→G 6.4% Y41H (TAT→CAT)  eutR ← eut operon transcriptional activator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,555,0600AG 6.4% 175.9 / 2.7 62Y41H (TAT→CAT) eutReut operon transcriptional activator, AraC family
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base G (2/2);  total (32/30)
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.48e-01

GGCACGACTGACGCAGCGCCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGTTAATGCCGGGTGTTGTAGGAC  >  NC_000913/2554913‑2555206
                                                                                                                                                   |                                                                                                                                                  
gtcctgACTGACGCACGGCCAGACCGGTGTATTCACGAAAAACCTGAATATCCTCGAGTAGAATTTCCGTAAACTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCTAGGCGGTAATCGTTAAAGCATGTTCATAg                                                                                                                                                   <  2:1150404/144‑1 (MQ=255)
  cACGACTGAGGCGGCGCCAGACCGTTGTATTCACGAAAAACCTGAAGATCATCGAGTAGAATTGCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGAc                                                                                                                                                 <  2:872615/149‑1 (MQ=255)
    cGACTGACGCAGCGCCAGACCGGTGGATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACAt                                                                                                                                               >  1:396682/1‑149 (MQ=255)
      aCTGAGGCAGCGCCAGACCGGTGTATGCAGGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATcc                                                                                                                                             <  1:281761/149‑1 (MQ=255)
       cTGACGCAGCGCCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCg                                                                                                                                            >  2:386171/1‑149 (MQ=255)
        tGACGCAGCGCCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGATTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGt                                                                                                                                           >  1:471727/1‑149 (MQ=255)
           cGCAGCGCCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTtgt                                                                                                                                        <  2:933239/149‑1 (MQ=255)
             cAGCGCCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTtggct                                                                                                                                      >  1:849374/1‑149 (MQ=255)
               gcgcCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCGGGCGGTAATCGTTAAAGCATGTTCATGGCCATCCGTtgtctgt                                                                                                                                    >  1:1129756/1‑149 (MQ=255)
               gcgcCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATGGACATCCGTtgtctgt                                                                                                                                    >  1:486603/1‑149 (MQ=255)
               gcgcCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTtgtctgt                                                                                                                                    >  1:103074/1‑149 (MQ=255)
               gcgcCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTtgtctgt                                                                                                                                    >  1:3952/1‑149 (MQ=255)
                tgcgtgACCGGTGTATCCACGAAACACGTGTATGTCGTGGGGGATCATTTCCGTAAAGTCGCCATGAAACTTGCCCGGATGCAGGTGACCGTAGTGCTGCGGCCAGGCGGTAATCGTGAAAGCATGTTCCTAGACCTCTGTtgtctgtt                                                                                                                                   <  2:341150/144‑2 (MQ=255)
                cgcaggACAGGTGTATTCACGCAACCCCAGAAAATCATCCATTAGAATTTCCGTAAATTACCCATGAAAATAGCCCGGATGCTGCTGATTGTAAATCTTCTGCCACGCTTGAATTGTTAAAGCATGTTCATAGACATCCGTtgtctgtc                                                                                                                                   <  2:1081082/144‑1 (MQ=255)
                  cctgACCGGTGTATTCCCGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGt                                                                                                                                 <  1:1093640/146‑1 (MQ=255)
                  ccAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGt                                                                                                                                 >  1:769872/1‑149 (MQ=255)
                      aCCCGTGTATTCCCGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTtgat                                                                                                                             <  2:1223083/149‑1 (MQ=255)
                         ggTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGaa                                                                                                                          <  1:910757/149‑1 (MQ=255)
                           tgtATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGTTGAACa                                                                                                                        >  1:860126/1‑149 (MQ=255)
                           tgtATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACa                                                                                                                        >  1:489639/1‑149 (MQ=255)
                                tCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACatcatc                                                                                                                   <  1:585239/149‑1 (MQ=255)
                                   cGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGGTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCAc                                                                                                                >  1:706385/1‑149 (MQ=255)
                                    gAAAAACCTGAATATCATGGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGATGCCAGGCGGGAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACt                                                                                                               >  2:65231/1‑149 (MQ=255)
                                                    atcGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGa                                                                                               >  1:422686/1‑149 (MQ=255)
                                                     tcGAGTAGAATTTCCGTAAATTCACCATTAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAg                                                                                              >  1:779912/1‑149 (MQ=255)
                                                          tAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCa                                                                                         >  2:679525/1‑149 (MQ=255)
                                                                  ccGTAAATTCACCATGAACCTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCg                                                                                 <  1:1218375/149‑1 (MQ=255)
                                                                    gTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGCGCTTCCGGTTTTCggg                                                                               >  2:310909/1‑149 (MQ=255)
                                                                                tGAGACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCCTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATg                                                                   >  1:267691/1‑149 (MQ=255)
                                                                                tGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTAATg                                                                   >  1:52743/1‑149 (MQ=255)
                                                                                tGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGATTTCGGGCAAGGTTTCATg                                                                   >  2:418136/1‑149 (MQ=255)
                                                                                tGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCAGAAAGCGGTAAAGCATGTTAATAGAAATCAGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGCGATTCAGGTTTTCGGGTACGGGTTCaag                                                                   >  2:218113/1‑147 (MQ=255)
                                                                                      ttGCCCGGATGCAGCTGATCGTACATCTGCAGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTGCGACCTTCGGCGTGGGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAg                                                             >  2:891004/1‑149 (MQ=255)
                                                                                       tGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGa                                                            <  2:471727/149‑1 (MQ=255)
                                                                                           cGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGt                                                        <  2:953558/149‑1 (MQ=255)
                                                                                                  aGCTGATCGTAAATCTGCTGCCAGGCGGTAGTCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCGGGATTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAAtt                                                 >  1:678469/1‑149 (MQ=255)
                                                                                                  aGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAAtt                                                 >  1:705878/1‑149 (MQ=255)
                                                                                                  aGCTGATCGTAAAGCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAAtt                                                 <  2:155149/149‑1 (MQ=255)
                                                                                                   gCTGTTCTTAAATCTTCGTCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATTCTTTGTCTGTCGTTAATGAACTTTTTCCACGTCGAGCTTCGGCGTGAGCTTCAGCTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTg                                                <  2:533952/149‑1 (MQ=255)
                                                                                                     tGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGc                                              >  1:238607/1‑149 (MQ=255)
                                                                                                        tCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATGGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGt                                           <  2:486603/149‑1 (MQ=255)
                                                                                                        tCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGt                                           <  2:3952/149‑1 (MQ=255)
                                                                                                        ccGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATGGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGt                                           <  2:1129756/148‑1 (MQ=255)
                                                                                                         cGTAAATCTGGCTCCAGGTGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTa                                          <  2:524914/149‑1 (MQ=255)
                                                                                                         cGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAGGTGGGTGTa                                          <  2:1189099/149‑1 (MQ=255)
                                                                                                         cGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTa                                          <  2:489639/149‑1 (MQ=255)
                                                                                                            aaaTCTGCTGCCAGGCCGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATTATAAAGATGGTGCAAATTGGCTGTACgg                                       <  1:695727/149‑1 (MQ=255)
                                                                                                                      gcAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCa                             <  2:1202948/148‑1 (MQ=255)
                                                                                                                      ccAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCa                             <  2:238607/149‑1 (MQ=255)
                                                                                                                        aGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCatg                           >  1:56560/1‑149 (MQ=255)
                                                                                                                            ggTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGt                       <  1:896957/149‑1 (MQ=255)
                                                                                                                             gTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGGGCTTCAGGTTTTTGGGTAAGGGTTCAGGACGTAGAGGGTGCAAATTGGTTGTTCGGGTCTTTTTCATGATGtt                      >  2:1175044/1‑149 (MQ=255)
                                                                                                                             gTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGtt                      >  2:709138/1‑149 (MQ=255)
                                                                                                                             gTAATCGCTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGtt                      <  1:485074/149‑1 (MQ=255)
                                                                                                                                 tCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCGAATTGGCTGTACGGGTCTTTTTCATGATGTTAATg                  <  2:779912/149‑1 (MQ=255)
                                                                                                                                   gTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTGATGATGTTACTGcc                >  1:337023/1‑149 (MQ=255)
                                                                                                                                    ttAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGTTAATGCCg               <  1:171816/149‑1 (MQ=255)
                                                                                                                                       aaGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGCCGGTACGGGTCTTTTTCATGATGTTAAGGCCgggg            >  2:560965/1‑148 (MQ=255)
                                                                                                                                            tGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGCGCAAATTGGCTGTACGGGTCTTTTTCATGATGTTAATGCCGGGtgttgt       >  2:851387/1‑149 (MQ=255)
                                                                                                                                                cATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGTTAATGCCGGGTGTTGTAGGa   <  2:242479/149‑1 (MQ=255)
                                                                                                                                                cATAGACATCCGTTGGCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGTTAATGCCGGGTGTTGTAGGa   <  2:242132/149‑1 (MQ=255)
                                                                                                                                                 aTAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTTATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGTTAATGCCGGGTGTTGTAGGac  >  2:1105573/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCACGACTGACGCAGCGCCAGACCGGTGTATTCACGAAAAACCTGAATATCATCGAGTAGAATTTCCGTAAATTCACCATGAAACTTGCCCGGATGCAGCTGATCGTAAATCTGCTGCCAGGCGGTAATCGTTAAAGCATGTTCATAGACATCCGTTGTCTGTCGTTGATGAACATTATCCACTTCGACCTTCGGCGTGAGCTTCAGGTTTTCGGGTAAGGGTTCATGATAAAGATGGTGCAAATTGGCTGTACGGGTCTTTTTCATGATGTTAATGCCGGGTGTTGTAGGAC  >  NC_000913/2554913‑2555206

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