Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,059,957 A→C 8.1% intergenic (+31/+7) asnU → / ← cbl tRNA‑Asn/ssuEADCB/tauABCD operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,059,9570AC 8.1% 215.0 / 9.9 62intergenic (+31/+7)asnU/cbltRNA‑Asn/ssuEADCB/tauABCD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base C (2/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

AAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGTCCCAACC  >  NC_000913/2059812‑2060103
                                                                                                                                                 |                                                                                                                                                  
aaaaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAttt                                                                                                                                                 >  3:138736/1‑149 (MQ=255)
 aaaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAATCCCGCTATAAAAGCGGGCtttt                                                                                                                                                >  2:21462/1‑149 (MQ=21)
 aaaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAtttt                                                                                                                                                >  1:66451/1‑149 (MQ=255)
   aGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAATCCCGCTATAAAAGCGGGCttttca                                                                                                                                              <  1:21462/149‑3 (MQ=21)
   aGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAATCCCGCTATAAAAGCGGGCttttca                                                                                                                                              <  2:388108/149‑3 (MQ=18)
      tttttCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCta                                                                                                                                           >  3:187577/1‑149 (MQ=255)
            gCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAATCCCGCTATAAAAGCGGGCttttca                                                                                                                                              >  1:109541/1‑138 (MQ=18)
            gCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAATCCCGCTATAAAAGCGGGCttttca                                                                                                                                              <  2:109541/140‑3 (MQ=18)
                        tgCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTc                                                                                                                         <  1:253508/149‑1 (MQ=255)
                            tCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctct                                                                                                                     >  3:389850/1‑149 (MQ=255)
                            tCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctct                                                                                                                     >  1:158007/1‑149 (MQ=255)
                            tCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctct                                                                                                                     >  4:304458/1‑149 (MQ=255)
                             caacaaCGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCtt                                                                                                                    <  4:138736/149‑1 (MQ=255)
                               acaacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCa                                                                                                                  >  4:415579/1‑149 (MQ=255)
                                       tCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTc                                                                                                          >  3:79755/1‑149 (MQ=255)
                                          tctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCAt                                                                                                       >  1:412468/1‑149 (MQ=255)
                                          tctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCAt                                                                                                       <  1:233243/149‑1 (MQ=255)
                                            tGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                                     >  2:195742/1‑149 (MQ=255)
                                                       cGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTg                                                                                          >  3:424426/1‑149 (MQ=255)
                                                       cGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTg                                                                                          <  3:304458/149‑1 (MQ=255)
                                                       cGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTg                                                                                          <  2:158007/149‑1 (MQ=255)
                                                          tAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGata                                                                                       >  1:33107/1‑149 (MQ=255)
                                                            gAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATc                                                                                     <  3:415579/149‑1 (MQ=255)
                                                              aCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCt                                                                                   <  2:412468/149‑1 (MQ=255)
                                                                     aCTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                            <  4:424426/149‑1 (MQ=255)
                                                                     aCTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                            >  3:3077/1‑149 (MQ=255)
                                                                                   tGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTAc                                                              >  4:24682/1‑149 (MQ=255)
                                                                                      cACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACaaa                                                           <  1:195742/149‑1 (MQ=255)
                                                                                       aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAg                                                          >  1:63561/1‑149 (MQ=255)
                                                                                        cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctc                                                                                                                      >  4:476070/1‑88 (MQ=255)
                                                                                        cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctc                                                                                                                      <  3:476070/88‑1 (MQ=255)
                                                                                            ttCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACaa                                                            <  4:122253/142‑1 (MQ=255)
                                                                                            ttCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACaa                                                            >  3:122253/1‑142 (MQ=255)
                                                                                            ttCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCa                                                     >  4:375290/1‑149 (MQ=255)
                                                                                                    cAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcgc                                             >  3:144874/1‑149 (MQ=255)
                                                                                                    cAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcgc                                             <  3:442008/149‑1 (MQ=255)
                                                                                                       tCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAg                                          >  2:53214/1‑149 (MQ=255)
                                                                                                       tCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAg                                          >  3:97346/1‑149 (MQ=255)
                                                                                                           aggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCTCTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCAtt                                                                <  2:141734/123‑1 (MQ=255)
                                                                                                           aggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCAtt                                                                >  1:141734/1‑123 (MQ=255)
                                                                                                             gagCCAAATTCCTGGAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAg                                    >  3:223695/1‑149 (MQ=255)
                                                                                                              agCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt                                   <  2:168263/149‑1 (MQ=255)
                                                                                                                 cAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTAc                                <  4:3077/149‑1 (MQ=255)
                                                                                                                  aaaTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACg                               >  3:92741/1‑149 (MQ=255)
                                                                                                                     ttCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTg                            >  1:220696/1‑149 (MQ=255)
                                                                                                                       ccTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGaa                          <  4:389850/149‑1 (MQ=255)
                                                                                                                       ccTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGaa                          <  2:66451/149‑1 (MQ=255)
                                                                                                                       ccTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGaa                          <  4:187577/149‑1 (MQ=255)
                                                                                                                       ccTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGCTCCAGAAAGCGCCAGACATAGTTACGCTGaa                          >  2:471661/1‑149 (MQ=255)
                                                                                                                        cTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAg                         >  4:224954/1‑149 (MQ=255)
                                                                                                                             aaGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTc                    <  3:224954/149‑1 (MQ=255)
                                                                                                                               gCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGCTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCt                  <  1:471661/149‑1 (MQ=255)
                                                                                                                                cccGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTAATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTTTCCTc                 >  4:311593/1‑149 (MQ=255)
                                                                                                                                 ccGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACAt                                      <  3:229484/127‑1 (MQ=255)
                                                                                                                                 ccGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACAt                                      >  4:229484/1‑127 (MQ=255)
                                                                                                                                  cGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGc               <  1:397253/149‑1 (MQ=255)
                                                                                                                                  cGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGc               <  4:97346/149‑1 (MQ=255)
                                                                                                                                    cTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCt                                                                                            >  2:421774/1‑70 (MQ=255)
                                                                                                                                    cTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCt                                                                                            <  1:421774/70‑1 (MQ=255)
                                                                                                                                           aGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGTccc      <  3:24682/149‑1 (MQ=255)
                                                                                                                                              gggATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGTCCCAAc   >  3:421297/1‑149 (MQ=255)
                                                                                                                                               ggATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGTCCCAAcc  <  4:79755/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCTTCAGTCCCAACC  >  NC_000913/2059812‑2060103

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