Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,770,539 C→T V281M (GTG→ATG)  yibH ← inner membrane protein YibH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,770,5390CT97.4% 122.0 / ‑4.5 38V281M (GTG→ATG) yibHinner membrane protein YibH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/14);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

ATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTC  >  NC_000913/3770301‑3770785
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
aTCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGcctc                                                                                                                                                                                                                                            >  1:44147/1‑248 (MQ=255)
       aaaCATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCAAGACAGGCAAAATACTAGTCAGTTTGCCGTGAAACACCTCCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGcc                                                                                                                                                                                                                                      >  2:128644/1‑250 (MQ=255)
          cATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCATGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGt                                                                                                                                                                                                                                  >  1:40528/1‑251 (MQ=255)
          cATGGCGGAAAGGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGt                                                                                                                                                                                                                                  <  1:111646/251‑1 (MQ=255)
                         cGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCCCGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGc                                                                                                                                                                                                                   >  1:52992/1‑251 (MQ=255)
                                       aCCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGaaa                                                                                                                                                                                                     <  2:100887/251‑1 (MQ=255)
                                          tGGGCGGAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTACTGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTCCCGTGAAACACCTGCCCCGGTAGCGCGTTAAACATCCCTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTTCCGCAATTg                                                                                                                                                                                                  >  2:11024/1‑251 (MQ=255)
                                           gggCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCAAACCGTCCGTGCAGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTgg                                                                                                                                                                                                 >  1:281325/1‑251 (MQ=255)
                                                    ggCCGGCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTTAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCAGCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAAtt                                                                                                                                                                                        <  2:52992/250‑1 (MQ=255)
                                                             gTAAGGCATCGATATCAGCGTTAGGGTCCAGTTCAATGGTGCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGtttt                                                                                                                                                                               <  2:241331/251‑1 (MQ=255)
                                                              tAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGttttt                                                                                                                                                                              >  1:84504/1‑251 (MQ=255)
                                                               aaGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCACCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTg                                                                                                                                                                             >  1:134455/1‑251 (MQ=255)
                                                                 ggCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTTCGCCTGCTAAGAACCGCCTGGCAAGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCt                                                                                                                                                                           >  1:19323/1‑251 (MQ=255)
                                                                               cGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTCTTTGCTCGGGGATGAAGa                                                                                                                                                               <  1:96945/249‑1 (MQ=255)
                                                                                            ttCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGc                                                                                                                                                 >  1:73051/1‑250 (MQ=255)
                                                                                                      cccatcccacCGTCCGTGCCGGGCAGGACCGTTAAGTATGGCAATACCCCCTGCGCCTTATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTGAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACGTCCGCATCATCGCCGGGTTTCAGACGTGACAGCGAGTTTTGCCGAAATGGTGCGACGAGTTGGCGTTGTTGCTCGGGGATGAAGACCGTCGCCGGACGCAGCGGCAAGGc                                                                                                                                      <  2:297203/244‑1 (MQ=255)
                                                                                                            ccacCGGCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGc                                                                                                                                <  2:258268/250‑1 (MQ=255)
                                                                                                                cgtccgtGCCGGGCACGACCGTTAATGAGTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTAt                                                                                                                            <  2:9081/251‑1 (MQ=255)
                                                                                                                                gACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATc                                                                                                             >  1:176469/1‑250 (MQ=255)
                                                                                                                                                   aCCCCCTGCGCCTTATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGTTAGCGCGTTAAACATCACTTCCTCATCAGCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGTCGGATCAGTACCTGAGTAACGTAGcc                                                                                         <  2:40528/251‑1 (MQ=255)
                                                                                                                                                   aCCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGcc                                                                                         <  2:25859/251‑1 (MQ=255)
                                                                                                                                                    cccccTGCGCCCGAAAAGAACCGCCTGGCACGACAGGTAAAATACTAGGCAGTTTGCCGTGAAACACCCGACCAGGTAGCGCGTTAAACATTCCTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAAATGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGGATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCa                                                                                        <  2:44147/251‑1 (MQ=255)
                                                                                                                                                           cgcCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCg                                                                                 >  1:205575/1‑251 (MQ=255)
                                                                                                                                                             ccTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGgc                                                                               <  1:38366/251‑1 (MQ=255)
                                                                                                                                                                       aCCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATg                                                                      >  1:12896/1‑250 (MQ=255)
                                                                                                                                                                          gCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACCATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGAAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGAc                                                                    >  2:111215/1‑249 (MQ=255)
                                                                                                                                                                           ccTGGAACGACAGGTAAAATACTAGTCAGTTTGCCGTGACCCACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCAGTTTTTGCTGGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGATCGGCGCGCGAATGACAGt                                                                 >  2:15226/1‑251 (MQ=255)
                                                                                                                                                                                cACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGCATGACAGTCTGCt                                                            >  2:13113/1‑251 (MQ=255)
                                                                                                                                                                                  cGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATTAAGACCATCACCGGACGCAGCg                                                                                                                                             <  2:228551/168‑1 (MQ=255)
                                                                                                                                                                                  cGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCg                                                                                                                                             >  1:228551/1‑168 (MQ=255)
                                                                                                                                                                                       ggTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGtt                                                     >  2:119181/1‑251 (MQ=255)
                                                                                                                                                                                        gTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTACACATCACTTCCGCATCATCGCCAGGTTTCAGCCGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGCATGACAGTCTGCTCACGGtt                                                     >  1:277160/1‑250 (MQ=255)
                                                                                                                                                                                             aTACTAGTCAGTTTGCCGTGAAACACCTGCCCAGATAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGCTGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGCCAGTCTGCTCAAGGTTATAttt                                               >  2:44239/1‑251 (MQ=255)
                                                                                                                                                                                                cTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGc                                            >  1:223160/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTTTGCCGTGTCACCCCTGCCCAGGTAGCGCGTTACACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGCGTTTTGCCGAAATTGGGCGCCAATTTGCCGTTTTTGCTCGGGGATGACGACCATCACCGGCCGCAGCGGCAAGGCAGCTGCGTCTGTACCTGGGCGGCTCAGCACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCCAGGTTATATTTTGCTTccgt                                      >  2:41308/1‑248 (MQ=255)
                                                                                                                                                                                                      aGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGt                                      >  1:154058/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccTGCCCAGGTAGCGCGTTAAACATCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTa                        >  2:125974/1‑249 (MQ=255)
                                                                                                                                                                                                                                          aaCATCACTTCCGCATCATCGCCAGGTATCAGACTGAACAGCGAGTTTTGCAGAAATAGGGCGACAATTTGCCGTTTGTGCTCGGGGATGACGACCATCACCGGACGCAGCGGCAAGGGAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGTCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGGAAGTTGCGCTCTTAGGCTCACAATCTGCTATTGCTc  <  2:185441/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      
ATCACCGAAACATGGCTGAAATGGTCGGAGTAAACCGCCACCTGGGCGTAGATGCCGTCGGGTAAGGCATCGATATCATCGTTAGGGTCCAGTTCAATGGTTCCCAGCACACCGTCCGTGCCGGGCACGACCGTTAATGATTGCAATACCCCCTGCGCCTGATAAGAACCGCCTGGCACGACAGGTAAAATACTAGTCAGTTTGCCGTGAAACACCTGCCCAGGTAGCGCGTTAAACACCACTTCCGCATCATCGCCAGGTTTCAGACGTAACAGCGAGTTTTGCCGAAATTGGGCGACAATTTGCCGTTTTTGCTCGGGGATGAAGACCATCACCGGACGCAGCGGCAAGGCAGCTGCGTATGTACCTGGGCGGATCAGTACCTGAGTAACGTAGCCATTGCTCGGCGCGCGAATGACAGTCTGCTCAAGGTTATATTTTGCTTCAGTAAGTTGCGCTCTTAAGCTCACAATCTGCGATTGCTC  >  NC_000913/3770301‑3770785

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: