Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,138,414 C→A 9.8% R173S (CGC→AGC)  yadG → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,138,4140CA 9.8% 174.1 / 9.0 61R173S (CGC→AGC) yadGABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base A (2/4);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

CTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAG  >  NC_002947/2138267‑2138553
                                                                                                                                                   |                                                                                                                                           
cTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGcgccg                                                                                                                                            >  8:290506/1‑149 (MQ=255)
   ggCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCt                                                                                                                                          <  1:166477/148‑1 (MQ=255)
    gCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCg                                                                                                                                        <  7:119315/149‑1 (MQ=255)
     ccTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGa                                                                                                                                       <  3:244084/149‑1 (MQ=255)
      cTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCGGCGCAGCTCGAt                                                                                                                                      <  8:7824/149‑1 (MQ=255)
      cTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGAGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCGGCGCAGCTCGAt                                                                                                                                      <  7:244050/149‑1 (MQ=255)
           ggACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGa                                                                                                                                 <  3:56797/149‑1 (MQ=255)
           ggACAAACGTGATGGGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCAGTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGa                                                                                                                                 <  5:70478/149‑1 (MQ=255)
            gACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAg                                                                                                                                <  6:143258/149‑1 (MQ=255)
                       tgtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTg                                                                                                                     <  6:346684/149‑1 (MQ=255)
                            aGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCt                                                                                                                 >  3:213999/1‑148 (MQ=255)
                            aGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCt                                                                                                                 >  2:285286/1‑148 (MQ=255)
                               cgcgTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAAc                                                                                                             <  3:189217/149‑1 (MQ=255)
                                      gttGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAgg                                                                                                      >  8:57411/1‑149 (MQ=255)
                                         gTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGc                                                                                                    <  4:213999/148‑1 (MQ=255)
                                          tCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCAt                                                                                                  <  5:187353/149‑1 (MQ=255)
                                           ccGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCAt                                                                                                  >  1:176490/1‑148 (MQ=255)
                                             ggcggcATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCAc                                                                                               >  5:39580/1‑149 (MQ=255)
                                                ggcATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACcat                                                                                            >  1:3055/1‑149 (MQ=255)
                                                   aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACcatcat                                                                                         >  1:192370/1‑149 (MQ=255)
                                                   aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACcatcat                                                                                         >  7:18320/1‑149 (MQ=255)
                                                   aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCAACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAACGGCATCACcatcat                                                                                         >  5:242379/1‑149 (MQ=255)
                                                    tGAAGCGGCGCCTGATGTTTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                        <  7:57411/149‑1 (MQ=255)
                                                    tGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                        <  7:290506/149‑1 (MQ=255)
                                                    tGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                        <  6:39580/149‑1 (MQ=255)
                                                     gAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCGGCGCAGCTCGAt                                                                                                                                      >  5:323818/1‑102 (MQ=255)
                                                     gAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCGGCGCAGCTCGAt                                                                                                                                      <  6:323818/102‑1 (MQ=255)
                                                            cgcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTcaccac                                                                                <  2:192370/149‑1 (MQ=255)
                                                                 gatgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCAccacccac                                                                            >  7:193234/1‑148 (MQ=255)
                                                                  atgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTa                                                                          >  7:179018/1‑149 (MQ=255)
                                                                    gatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCCCCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCAATGAGCTGAACCAGAAGGGCATCACCATCATCATAACCCACCACtctc                                                                        >  4:246584/1‑149 (MQ=255)
                                                                     atTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGGGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCt                                                                       <  5:330107/149‑1 (MQ=255)
                                                                        gCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCAGTGAGCTGAACGAGAAGGGCATCACCATCATCCTCACCACCCACTATCTgg                                                                     <  8:179018/148‑1 (MQ=255)
                                                                        gCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAGCCAGACCGGCATCACCATCCTCCTGCCCAGCCACTATCTgg                                                                     >  8:306872/1‑148 (MQ=255)
                                                                                aTTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCt                                                             >  8:336032/1‑148 (MQ=255)
                                                                                     tCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGGATAGCCATCATCCTCACCACCCACTATCTGGAAGAGGCTGagcag                                                       >  1:313003/1‑149 (MQ=255)
                                                                                     tCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGagcag                                                       <  7:336032/149‑1 (MQ=255)
                                                                                        aCGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtg                                                    >  7:79983/1‑149 (MQ=255)
                                                                                         cGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGTGCTGCGCCGCTCGCTGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTACCCACCCACGATCTGGAAGAGGCTGAGCAGCtgt                                                   >  1:92805/1‑149 (MQ=255)
                                                                                         cGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtgt                                                   <  1:146538/149‑1 (MQ=255)
                                                                                            aCCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGcc                                                >  1:299859/1‑149 (MQ=255)
                                                                                            aCCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGcc                                                >  2:103152/1‑149 (MQ=255)
                                                                                            aCCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGcc                                                <  7:66975/149‑1 (MQ=255)
                                                                                             ccGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCg                                               <  2:176490/149‑1 (MQ=255)
                                                                                                           ccTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGcatca                                 >  6:256818/1‑149 (MQ=255)
                                                                                                             tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCatcatc                               >  4:8008/1‑149 (MQ=255)
                                                                                                             tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCatcatc                               >  8:4888/1‑149 (MQ=255)
                                                                                                                   tatccACCCCGGCGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCAt                         >  6:179148/4‑149 (MQ=255)
                                                                                                                     cccACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCg                              >  1:256433/1‑142 (MQ=255)
                                                                                                                     cccACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCg                              <  2:256433/142‑1 (MQ=255)
                                                                                                                     cccACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATgg                       <  6:195368/149‑1 (MQ=255)
                                                                                                                      ccACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGag                                                          >  1:245744/1‑113 (MQ=255)
                                                                                                                      ccACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGag                                                          <  2:245744/113‑1 (MQ=255)
                                                                                                                            ccgGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCAtc                <  3:202225/149‑1 (MQ=255)
                                                                                                                                 gTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCga           <  7:306872/149‑1 (MQ=255)
                                                                                                                                  tccacatCGAGCGGCGCAGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCAtcgtc             >  4:200487/6‑146 (MQ=255)
                                                                                                                                      tatCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCga           <  7:320558/144‑1 (MQ=255)
                                                                                                                                      tatCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCga           >  8:320558/1‑144 (MQ=255)
                                                                                                                                          gAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCa   <  7:4888/148‑1 (MQ=255)
                                                                                                                                          gAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAg  <  2:299859/149‑1 (MQ=255)
                                                                                                                                          gAGCTGCGCCGCGCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCa   <  7:128184/148‑1 (MQ=255)
                                                                                                                                          gAGCGGCGCAGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCa   <  3:200487/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAG  >  NC_002947/2138267‑2138553

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