Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,895,760 T→A 9.4% L57Q (CTG→CAG)  PP_4304 → VIC family cation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,895,7600TA 9.4% 191.0 / 11.2 64L57Q (CTG→CAG) PP_4304VIC family cation transporter
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base A (3/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

GCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACT  >  NC_002947/4895613‑4895909
                                                                                                                                                    |                                                                                                                                                     
gcgAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTgg                                                                                                                                                     >  1:621721/1‑150 (MQ=255)
       gCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCg                                                                                                                                                  <  1:825489/146‑1 (MQ=255)
       gCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCg                                                                                                                                                  >  2:825489/1‑146 (MQ=255)
                cATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCATCCAAAATA‑CTGGTGTGTGGATCATCTTGGCCAGCACGGTGACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTgggg                                                                                                                                     <  1:928818/150‑1 (MQ=255)
                   cGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCtt                                                                                                                                  <  2:1594525/150‑1 (MQ=255)
                    gTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCttt                                                                                                                                 >  2:407559/1‑150 (MQ=255)
                               cGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGAGCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCtt                                                                                                                      >  1:894844/1‑150 (MQ=255)
                                gACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACcg                                                                                                                             >  2:363428/1‑142 (MQ=255)
                                gACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACcg                                                                                                                             <  1:363428/142‑1 (MQ=255)
                                  cacCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCt                                                                                                                   <  2:1140877/150‑1 (MQ=255)
                                   agcGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTg                                                                                                                  >  2:1643493/3‑150 (MQ=255)
                                   acCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTATCGGTAATTCTCTACACCCTCCATGAGGGGCACCAGGGCTACGCCGGCCTGCTTGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTg                                                                                                                  <  1:115913/150‑1 (MQ=255)
                                      gTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGcc                                                                                                               <  1:506834/150‑1 (MQ=255)
                                         gCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTc                                                                                                                     <  1:312366/141‑1 (MQ=255)
                                         gCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTc                                                                                                                     >  2:312366/1‑141 (MQ=255)
                                           cGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGATTa                                                                                                          <  1:1643493/150‑1 (MQ=255)
                                            ggTCGCCGCTTCGACAAGATCCTGT‑TGCTGCTCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATAc                                                                                                         <  2:621721/150‑1 (MQ=255)
                                              tCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACAt                                                                                                       <  1:407559/150‑1 (MQ=255)
                                                        gACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGc                                                                                                >  2:1188633/1‑147 (MQ=255)
                                                        gACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGc                                                                                                <  1:1188633/147‑1 (MQ=255)
                                                                       gCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCaa                                                                               >  2:679297/1‑150 (MQ=255)
                                                                        cTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATc                                                                                                                        <  1:923752/108‑1 (MQ=255)
                                                                        cTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATc                                                                                                                        >  2:923752/1‑108 (MQ=255)
                                                                          gATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGcc                                                                            >  1:1218481/1‑150 (MQ=255)
                                                                           atcatcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTgg                                                                                                                                       >  1:1310782/1‑90 (MQ=255)
                                                                           atcatcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTgg                                                                                                                                       <  2:1310782/90‑1 (MQ=255)
                                                                           atcatcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCAGCAGGCCggcgt                                                                                                                                              >  2:1223633/1‑80 (MQ=255)
                                                                           atcatcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCAGCAGGCCggcgt                                                                                                                                              <  1:1223633/83‑4 (MQ=255)
                                                                               tcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTgc                                                                       >  2:1499997/1‑150 (MQ=255)
                                                                                ccTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTgcg                                                                      >  1:502296/1‑150 (MQ=255)
                                                                                           tGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGt                                                                                                                                         <  1:595520/72‑1 (MQ=255)
                                                                                           tGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGt                                                                                                                                         >  2:595520/1‑72 (MQ=255)
                                                                                                  ggTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAg                                                                                                                                          <  1:633256/64‑1 (MQ=255)
                                                                                                  ggTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAg                                                                                                                                          >  2:633256/1‑64 (MQ=255)
                                                                                                     gATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCa                                                           >  1:1353596/1‑140 (MQ=255)
                                                                                                     gATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCa                                                           <  2:1353596/140‑1 (MQ=255)
                                                                                                      aTCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCAGATAGCCGGCAACAAGTTTGGCGTGGTCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGc                                                <  1:603654/150‑1 (MQ=255)
                                                                                                       tCCTCGACAGCATCGATGGGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTCTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGcc                                               >  2:826040/1‑150 (MQ=255)
                                                                                                       tCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGcc                                               <  2:502296/150‑1 (MQ=255)
                                                                                                        ccTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGTCCt                                              >  1:543419/1‑150 (MQ=255)
                                                                                                        ccTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCt                                              >  1:1166072/1‑150 (MQ=255)
                                                                                                          tCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTgg                                            >  1:976195/1‑150 (MQ=255)
                                                                                                               aGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGAt                                       <  2:1218481/150‑1 (MQ=255)
                                                                                                                 cATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTc                                                                                                                     <  1:1148492/70‑1 (MQ=255)
                                                                                                                 cATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTc                                                                                                                     >  2:1148492/1‑70 (MQ=255)
                                                                                                                  aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGc                                                                                   <  1:889627/103‑1 (MQ=255)
                                                                                                                  aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGc                                                                                   >  2:889627/1‑103 (MQ=255)
                                                                                                                  aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTa                                    >  2:902432/1‑150 (MQ=255)
                                                                                                                     gATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTg                                 <  2:976195/150‑1 (MQ=255)
                                                                                                                     gATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTg                                 <  1:1499997/150‑1 (MQ=255)
                                                                                                                      aTGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCt                                                                                                                   <  1:650407/67‑1 (MQ=255)
                                                                                                                      aTGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCt                                                                                                                   >  2:650407/1‑67 (MQ=255)
                                                                                                                               aCCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGc                       >  2:224075/1‑150 (MQ=255)
                                                                                                                                ccAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGcc                      >  2:1156768/1‑150 (MQ=255)
                                                                                                                                  aGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCgg                    <  1:902432/150‑1 (MQ=255)
                                                                                                                                    ggCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTgc                                                                       <  2:1339885/97‑1 (MQ=255)
                                                                                                                                    ggCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTgc                                                                       >  1:1339885/1‑97 (MQ=255)
                                                                                                                                      cgataccGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCggg                   <  2:1299963/142‑1 (MQ=255)
                                                                                                                                      cgataccGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCggg                   >  1:1299963/6‑147 (MQ=255)
                                                                                                                                       tACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCgg                    >  1:468532/1‑145 (MQ=255)
                                                                                                                                       tACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCgg                    <  2:468532/145‑1 (MQ=255)
                                                                                                                                         taccGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGt                                                         <  1:1315182/104‑1 (MQ=255)
                                                                                                                                         taccGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGt                                                         >  2:1315182/3‑106 (MQ=255)
                                                                                                                                              gCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCgcgc        <  2:894844/150‑1 (MQ=255)
                                                                                                                                                    tgGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTgg                                            >  2:1139992/1‑108 (MQ=255)
                                                                                                                                                    tgGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTgg                                            <  1:1139992/108‑1 (MQ=255)
                                                                                                                                                    tgGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGtact  >  1:212462/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGT‑TGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACT  >  NC_002947/4895613‑4895909

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑