Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,987,992 C→A 7.6% R152R (CGC→CGA PP_1779 → lipopolysaccharide ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,987,9920CA 7.6% 185.8 / 4.0 63R152R (CGC→CGAPP_1779lipopolysaccharide ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base A (3/2);  total (33/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 = 9.71e-01

GTGTACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGATAGTGT  >  NC_002947/1987855‑1988137
                                                                                                                                         |                                                                                                                                                  
gtgtACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGACGCGCATTTGc                                                                                                                                         <  8:261275/149‑1 (MQ=255)
gtgtACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGACGCGCATTTGc                                                                                                                                         <  3:176708/149‑1 (MQ=255)
      cACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTTAAGCCTTTGACGAGATTGGCGAGTACTTTGATGAGGAAGTTCGTACCTATTCCAGCGGTATCCAAATCCGCGTCGCATTTCCCGtcca                                                                                                                                   >  8:254033/1‑146 (MQ=255)
         gCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCATTGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATCCGCGTCGCATTTGCCGTGGCGAc                                                                                                                                <  1:125634/149‑1 (MQ=255)
             cGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGc                                                                                                                             <  1:74721/148‑1 (MQ=255)
              gggTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGt                                                                                                                           >  7:32635/1‑149 (MQ=255)
                             cGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTcc                                                                                                                    >  8:76400/1‑141 (MQ=255)
                             cGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTcc                                                                                                                    <  7:76400/141‑1 (MQ=255)
                               aTGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCt                                                                                                          >  3:351370/1‑148 (MQ=255)
                                  tACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTAtc                                                                                                      >  8:45178/1‑149 (MQ=255)
                                    cagcagATCAATGATGTCATTGATGACATTGAAGCCTTTGACGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTAtcgt                                                                                                    >  2:41188/1‑149 (MQ=255)
                                     agcagATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTAtcgtc                                                                                                   >  4:261818/1‑149 (MQ=255)
                                        agATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGAt                                                                                                <  4:351370/149‑1 (MQ=255)
                                          aTCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGAATTTGCCGTGGCGCATGCGTTCCGTCCAGAAATCCTTATCGTCGAt                                                                                                >  3:119354/1‑147 (MQ=255)
                                              atgatgTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGgcg                                                                                          >  2:205483/1‑149 (MQ=255)
                                                    tCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGCAAATCCTTATCGTCGATGAGGCGCTGt                                                                                      <  1:205483/148‑1 (MQ=255)
                                                    tCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGt                                                                                      >  1:257127/1‑147 (MQ=255)
                                                          atgaCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGACGCGCATTTGc                                                                                                                                         >  2:196579/1‑91 (MQ=255)
                                                           tgaCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAggg                                                                              <  7:45178/148‑1 (MQ=255)
                                                           tgaCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAggg                                                                              <  2:102254/148‑1 (MQ=255)
                                                                ttGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGc                                                                        >  6:205216/1‑149 (MQ=255)
                                                                   aaGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTAt                                                                     >  2:75373/1‑149 (MQ=255)
                                                                    aGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTAtt                                                                    >  4:31897/1‑149 (MQ=255)
                                                                      ccTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTc                                                                  <  3:31897/149‑1 (MQ=255)
                                                                        tttGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGAGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCa                                                                >  3:233050/1‑149 (MQ=255)
                                                                          tGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCa                                                                <  2:257127/147‑1 (MQ=255)
                                                                             cGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGt                                                                                  <  4:114333/126‑1 (MQ=255)
                                                                             cGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCa                                                                >  8:121787/1‑144 (MQ=255)
                                                                             cGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCa                                                                <  7:121787/144‑1 (MQ=255)
                                                                                       gtgTACTTTGATGAGTGAGTGCGTACCTATTCAAGCGGTATGCAAATGCGCTTCGCATTTTCCTTTGCGACTGCGTTCCGTCCAGAAATACTTATAGTCTATGAGGCGCTGTCATTAGGTTACAGCTATTTCCAGCACAAGAGTTTTGa                                                 <  4:119354/147‑1 (MQ=255)
                                                                                       gAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGtttt                                                   <  1:250380/147‑1 (MQ=255)
                                                                                       gAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGtttt                                                   >  2:250380/1‑147 (MQ=255)
                                                                                       gAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGa                                                 <  1:41188/149‑1 (MQ=255)
                                                                                         gTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAgc                                               <  5:205216/149‑1 (MQ=255)
                                                                                           aCTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACaa                                                          <  8:15117/136‑1 (MQ=255)
                                                                                           aCTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACaa                                                          >  7:15117/1‑136 (MQ=255)
                                                                                             tttGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAaga                                                        <  8:155738/136‑1 (MQ=255)
                                                                                             tttGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAaga                                                        >  7:155738/1‑136 (MQ=255)
                                                                                             tttGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCAt                                           >  6:282407/1‑149 (MQ=255)
                                                                                                gatgaGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCg                                        >  7:106238/1‑149 (MQ=255)
                                                                                                    aGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAg                                    >  3:37751/1‑149 (MQ=255)
                                                                                                    aGGCAGTGAGTACCTATTCCAGCGGTATGCAAATGAGCGTCGCATTTGCAGTGGCGACTGCGTTACGGCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGa                                     <  7:327877/148‑1 (MQ=255)
                                                                                                     ggCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGGGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGtttt                                                   >  1:247137/1‑133 (MQ=255)
                                                                                                     ggCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGtttt                                                   <  2:247137/133‑1 (MQ=255)
                                                                                                        aGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTc                                 >  6:117860/1‑148 (MQ=255)
                                                                                                         gTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAgc                                               >  6:252475/1‑133 (MQ=255)
                                                                                                             gTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATcc                                         >  7:293043/1‑135 (MQ=255)
                                                                                                             gTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATcc                                         <  8:293043/135‑1 (MQ=255)
                                                                                                             gTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCgg                           >  1:255507/1‑149 (MQ=255)
                                                                                                                ccTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTg                        <  1:75373/149‑1 (MQ=255)
                                                                                                                 cTATTCCAGCTTTATGGAAATGCTCTTCGGATTTGGCGTGGCGACTGCGTTCCGTCCATAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTTTTTCCAGCACAATAGTGTTTAGCGCATCCGTGAGTTCCAGCGGCTgg                       <  7:254033/149‑1 (MQ=255)
                                                                                                                  tATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCAAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTgg                       >  1:102408/1‑148 (MQ=255)
                                                                                                                       cAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCa                               >  1:318773/1‑135 (MQ=255)
                                                                                                                       cAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCa                               <  2:318773/135‑1 (MQ=255)
                                                                                                                         gCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCATTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGATTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCg               >  2:97064/1‑149 (MQ=255)
                                                                                                                            gTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGt                                       <  4:280367/122‑1 (MQ=255)
                                                                                                                            gTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGt                                       >  3:280367/1‑122 (MQ=255)
                                                                                                                            cacggCAAATGCGACGCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCtt            >  7:86997/5‑149 (MQ=255)
                                                                                                                             tATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCttt           <  8:106238/149‑1 (MQ=255)
                                                                                                                             tATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCttt           <  5:282407/149‑1 (MQ=255)
                                                                                                                             tATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCttt           <  3:261818/149‑1 (MQ=255)
                                                                                                                                gCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGa        >  5:44607/1‑149 (MQ=255)
                                                                                                                                 cAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGa        <  8:32635/148‑1 (MQ=255)
                                                                                                                                       cgcgTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGATAgtgt  <  6:44607/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
GTGTACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGATAGTGT  >  NC_002947/1987855‑1988137

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