Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,831,878 A→C 15.8% V89G (GTC→GGC)  ftsY ← signal recognition particle protein translocation receptor GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,831,8780AC15.8% 124.7 / 2.7 58V89G (GTC→GGC) ftsYsignal recognition particle protein translocation receptor GTPase
Reads supporting (aligned to +/- strand):  ref base A (28/20);  new base C (2/8);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 3.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGACTCCACAGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTCGACAGCGGGGGGTT  >  NC_002947/5831736‑5832025
                                                                                                                                               |                                                                                                                                                   
ggCAGACTCCACAGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCg                                                                                                                                               <  4:53458/149‑1 (MQ=255)
  cAGACTCCACAGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCg                                                                                                                                               >  5:236580/1‑147 (MQ=255)
  cAGACTCCACAGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGa                                                                                                                                             <  8:114051/149‑1 (MQ=255)
  cAGACTCCACAGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGa                                                                                                                                             <  3:349295/149‑1 (MQ=255)
  cAGACTCCACAGGGGCCGCCTGCAAGGGCACGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGa                                                                                                                                             <  5:202701/149‑1 (MQ=255)
       tCCACAGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCg                                                                                                                                        >  4:318767/1‑149 (MQ=255)
         cacaGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCggg                                                                                                                                      >  4:163829/1‑149 (MQ=255)
           caGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACAGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCggggg                                                                                                                                    >  4:167187/1‑149 (MQ=255)
            aGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCgggggg                                                                                                                                   <  8:346219/149‑1 (MQ=255)
               ggCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTg                                                                                                                                <  4:161016/149‑1 (MQ=255)
                 ccgccTGCACGGGCCGGGACGCATCGGGCTGGGGGGCAGGTGCCACCAGGGGCTGGGTTGCGACGGGGT‑CTGGCGCTGGCGCCGGCATCGGCTGCGCCGCAACCGGTTCTGGGAGGGGGGCCTCGACTACCGGAGGCCGGGGGGGTgcg                                                                                                                              <  4:263638/149‑1 (MQ=255)
                       gCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGAcc                                                                                                                        >  1:229365/1‑149 (MQ=255)
                           gggCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCt                                                                                                                    >  1:124279/1‑149 (MQ=255)
                                   aCGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCGGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACGGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGGGCGGCGACCGGCTGCGGCTca                                                                                                            <  2:249746/149‑1 (MQ=255)
                                     gCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGCCTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTcagc                                                                                                          >  5:215607/1‑149 (MQ=255)
                                       aaCCGGCTGCGGGGCAGGTGCCCCCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCACCGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACGACCGGAGGCCGGGGGGGTGCGGCGCCCGGCGGCGGCTccgacc                                                                                                        >  4:251512/1‑144 (MQ=255)
                                       aaCCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTcagc                                                                                                          >  5:104303/1‑147 (MQ=255)
                                                  gggCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAGCGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGcc                                                                                             >  2:238395/1‑149 (MQ=255)
                                                   ggCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCg                                                                                            >  8:177093/1‑149 (MQ=255)
                                                    tgggggcccacccGGGGCTGGGCTGCGCCCGGCG‑CTGGCGCAGGCGCCTCCCCCGGCGGCGCCGCCCCCGGTGCCGGGCCGGGGGCCTCGACTCCGGGCGGCCGGGGGGGGGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGc                                                                                           <  6:271072/142‑1 (MQ=255)
                                                      aGGTGCCACCAGCGGCTGGGATGCGACCGGCTGATGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGCCGGGGGCCGCGCCTACCGGAGGCCGGGGGGGTGCGGCGACCGGATGCGGGTCAGCCAGGTGTGGGGCCGc                                                                                           >  1:233996/1‑148 (MQ=255)
                                                        gTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGc                                                                                           >  5:271072/1‑145 (MQ=255)
                                                        gTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACg                                                                                       >  7:298035/1‑149 (MQ=255)
                                                          gCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGc                                                                                     >  5:194936/1‑149 (MQ=255)
                                                                 gCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACg                                                                              >  1:123345/1‑149 (MQ=255)
                                                                          aggcgcCCGGCG‑CGGGCTCCGGCGCCGGCCCCGGCTGCGCCGCCCCCGGGTCCGGGCCCGGGGCCGCGCCTCCCGGCGGCCGGGGGGGGGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCgggg                                                                     <  3:251512/147‑1 (MQ=255)
                                                                            gcgcCCGGCC‑CTGGCGCCGCCGCCCGCAACGGCCGCGACGCCACCGGTTCAGGGACGGGGGCCTCGACTACCGGCGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGc                                                                   <  4:238440/149‑1 (MQ=255)
                                                                             cgcccGGCG‑CTGGCTCAGCCGCCTGCACCGGCTGCGCCGCAACCGGTTCAGGGCCGGGGGCCTCCACACCCGGCGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGc                                                                   <  2:233996/145‑1 (MQ=255)
                                                                                ccGGCG‑CTGGCGCCGGCGCCTGCCCCGGCGGCGCCGCAACCGGTGCCGGGCCGGGGCCCCCGACCCCCGCCGGCCGGGGGGGGGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCtgt                                                               <  1:238395/149‑1 (MQ=255)
                                                                                    gcgcgggCTCCGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGCCGGGGGCCTCGACTCCCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGgc                                                            <  6:186910/144‑1 (MQ=255)
                                                                                     t‑ctGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGgcgc                                                          >  8:204448/1‑149 (MQ=255)
                                                                                            tCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGgcgcctgagc                                                    >  6:359951/1‑143 (MQ=255)
                                                                                                       cAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCt                                         >  7:46665/1‑149 (MQ=255)
                                                                                                       cAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCt                                         >  3:331328/1‑149 (MQ=255)
                                                                                                       cAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGCGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCt                                         >  7:270108/1‑149 (MQ=255)
                                                                                                        cccGGCTGCGCCTCACCCGGGGCCGGGCCGGGGGCCTCGACTACCGGCGGCGGGGGGGGTGCGGCGACCGGCTGCGGCGCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCt                                         <  8:283000/146‑1 (MQ=255)
                                                                                                        aaCGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCt                                         >  7:283000/1‑148 (MQ=255)
                                                                                                          cGGCTGCGACGCAACCGGTTCCGGGACGGGGGCCTCGACTACCGGCGGCGGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGa                                      <  8:46665/149‑1 (MQ=255)
                                                                                                               gCGACGCACCCGGTGCCGGGCCGGGGGCCGCGACTCCCGGAGGCCGGGGGGGGGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGa                                 <  2:229365/149‑1 (MQ=255)
                                                                                                                 gACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGc                               >  3:342132/1‑149 (MQ=255)
                                                                                                                    gCAACCGGTTCCGGGCCGGGGGCCGCGACTACCGGCGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGc                            <  6:236580/149‑1 (MQ=255)
                                                                                                                    gCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGc                            >  5:101001/1‑149 (MQ=255)
                                                                                                                     cacccGGTTCAGGGACGGGGGCCGCGCCTACCGGCGGCGGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCt                           <  5:5305/146‑1 (MQ=255)
                                                                                                                     cAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCt                           >  6:248062/1‑149 (MQ=255)
                                                                                                                      aaCCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTg                          >  3:59368/1‑149 (MQ=255)
                                                                                                                        ccGGTGCCGGGACGGGGCCCGCGCCTCCCGGCGGCCGGGGGGGCGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCa                        <  7:215416/149‑1 (MQ=255)
                                                                                                                          gttgcAGGGCGGGGGGCCTCGCCCACCGGTGGCGGGGGGGGGGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGa                      <  2:124279/145‑1 (MQ=255)
                                                                                                                          ggtgccgggCCGGGGGCCGCCCCTCCCGGCGGCCGGGGGGGGGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCCGCAGa                      <  6:194936/143‑1 (MQ=255)
                                                                                                                          ggTTCCGGGCCGGGGGCCTCGACTCCCGGCGGCCGGGGGGGGGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGa                      <  3:356260/149‑1 (MQ=255)
                                                                                                                             tCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTc                                                       >  3:341182/1‑113 (MQ=255)
                                                                                                                                gggACGGGGGCCTCTCCTCCCGGCGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTc                <  3:287126/149‑1 (MQ=255)
                                                                                                                                gggACGGGGGCCGCGCCTACCGGTGGCGGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTc                <  8:298035/149‑1 (MQ=255)
                                                                                                                                  gACGGGGGCCTCGCCTCCCGGCGGCCGGGGGGGGGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTCGa              <  3:318767/149‑1 (MQ=255)
                                                                                                                                      ggggCCTCGACTCCCGGCGGGGGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTCGACAGc          <  6:101001/149‑1 (MQ=255)
                                                                                                                                       gggCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTCGACAGCg         >  3:45084/1‑149 (MQ=255)
                                                                                                                                       gggCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTCGACAGCg         >  8:158027/1‑149 (MQ=255)
                                                                                                                                              gACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTCGACAGCGGGGGGtt  >  3:250432/1‑149 (MQ=255)
                                                                                                                                              gACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTCGACAGCGGGGGGt   <  4:45084/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGCAGACTCCACAGGGGCCGCCTGCAAGGGCATGGACGCAACCGGCTGCGGGGCAGGTGCCACCAGCGGCTGGGATGCGACCGGCT‑CTGGCTCAGGCGCCTGCAACGGCTGCGACGCAACCGGTTCAGGGACGGGGGCCTCGACTACCGGAGGCCGGGGGGGTGCGGCGACCGGCTGCGGCTCAGCCAGCTGTGGGGCCGCAACGGCCTGAACGGGTTCGGGGGCCTGTGGCGCTTCGGCGGCAGATCCCTCGACCGGAGCGGCTGCAGACTGCTCGACAGCGGGGGGTT  >  NC_002947/5831736‑5832025

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