Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 318,979 T→C 7.9% I255V (ATC→GTC) ‡ PP_0264 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947318,9790TC 7.9% 179.2 / 5.0 63I255V (ATC→GTC) ‡PP_0264sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base C (2/3);  total (30/33)
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.34e-01

GCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCA  >  NC_002947/318832‑319119
                                                                                                                                                   |                                                                                                                                            
gCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGCCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                             >  3:79374/1‑149 (MQ=255)
 cATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                             <  7:26553/148‑1 (MQ=255)
 cATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATc                                                                                                                                             <  6:298749/148‑1 (MQ=255)
     aCCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAgac                                                                                                                                        >  7:205179/1‑147 (MQ=255)
           ccTGCCATCAACATCACTCCCAGCAGCGTACCGGCACCGATCGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCcgcga                                                                                                                                  <  8:199629/149‑2 (MQ=255)
               ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGTCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTg                                                                                                                              <  2:118305/149‑1 (MQ=255)
               ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGtt                                                                                                                               >  4:32362/1‑148 (MQ=255)
               ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCcgcgacatc                                                                                                                              <  8:205179/149‑6 (MQ=255)
                               cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtgtg                                                                                                              >  4:209098/1‑149 (MQ=255)
                               cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtgtg                                                                                                              >  4:292135/1‑149 (MQ=255)
                               cagcacCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCtgtg                                                                                                              >  7:10717/1‑149 (MQ=255)
                                   agggTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGc                                                                                                          <  7:256694/146‑1 (MQ=255)
                                       tACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTgg                                                                                                      <  1:152785/149‑1 (MQ=255)
                                        aCCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTggg                                                                                                     >  1:149889/1‑149 (MQ=255)
                                         ccGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTggg                                                                                                     <  8:24141/148‑1 (MQ=255)
                                          cGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTc                                                                                                   >  6:155934/1‑149 (MQ=255)
                                          cGGCACCGATGGCCATTACCGTCTGCACCAAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGCGCACCAGCGGCTGTGGGGCCTGGGTc                                                                                                   >  4:223710/1‑149 (MQ=255)
                                            gCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCg                                                                                                  >  1:146236/1‑148 (MQ=255)
                                               ccGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTg                                                                                                       >  4:240085/1‑140 (MQ=255)
                                               ccGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTg                                                                                                       <  3:240085/140‑1 (MQ=255)
                                               ccGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGa                                                                                                >  1:142817/1‑147 (MQ=255)
                                                  aTGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAggg                                                                                           >  3:61285/1‑149 (MQ=255)
                                                   tGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCCGc                                                                                                                                    <  6:228782/107‑1 (MQ=255)
                                                   tGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCAGGCCTGGATCACCCAGACCCGc                                                                                                                                    >  5:228782/1‑107 (MQ=255)
                                                       cATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGa                                                                                      >  8:207389/1‑149 (MQ=255)
                                                        aTTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGAt                                                                                     >  3:157994/1‑149 (MQ=255)
                                                           aCCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATgg                                                                                   <  5:155934/148‑1 (MQ=255)
                                                           aCCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATgg                                                                                   <  4:289102/148‑1 (MQ=255)
                                                              gTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgat                                                                               <  3:283749/149‑1 (MQ=255)
                                                               tCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgat                                                                               >  4:165163/1‑148 (MQ=255)
                                                                cTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatga                                                                             >  1:163159/1‑149 (MQ=255)
                                                                   cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGNCCTGGGTCGGATAGGGCTGGATGGCgatga                                                                             <  1:30471/146‑1 (MQ=255)
                                                                   cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatga                                                                             <  3:292135/146‑1 (MQ=255)
                                                                   cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatga                                                                             <  4:61285/146‑1 (MQ=255)
                                                                      cGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCt                                                                       <  3:165163/149‑1 (MQ=255)
                                                                       gAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTg                                                                      <  2:163159/149‑1 (MQ=255)
                                                                        aGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTgg                                                                     <  4:79374/149‑1 (MQ=255)
                                                                            aTCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCt                                                                 <  1:123189/149‑1 (MQ=255)
                                                                               aaCCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCATGGCGCTCgg                                                              >  1:271360/1‑149 (MQ=255)
                                                                                 ccAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTTCGGGTTGAGAACAAGCGGCTGTCGGGACTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCg                                                            >  7:217566/1‑149 (MQ=255)
                                                                                  cAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGt                                                           >  8:241828/1‑149 (MQ=255)
                                                                                    gCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCg                                                         <  3:209098/149‑1 (MQ=255)
                                                                                       tGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTc                                                      <  2:142817/149‑1 (MQ=255)
                                                                                         cgcgGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAg                                                    <  7:241828/149‑1 (MQ=255)
                                                                                             gggCGAGGATGCTGACTAGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTg                                                >  4:60146/1‑149 (MQ=255)
                                                                                              ggCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGTGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGc                                               <  7:207389/149‑1 (MQ=255)
                                                                                               gCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCg                                              <  2:149889/149‑1 (MQ=255)
                                                                                                        cTGCCTTGCCAACCAGTTTCCTTTCTGTCGCGGGTCTGGGTGATCCAGGCCTGCGGTTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCAc                                      <  8:217566/148‑1 (MQ=255)
                                                                                                         tGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACg                                     >  8:14855/1‑148 (MQ=255)
                                                                                                           aCTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGa                                   <  7:14855/148‑1 (MQ=255)
                                                                                                               gCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgc                                <  3:60146/147‑1 (MQ=255)
                                                                                                                ccAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgcgc                              <  2:146236/148‑1 (MQ=255)
                                                                                                                 cAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAgcgc                              >  4:75537/1‑147 (MQ=255)
                                                                                                                  aaCCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCg                                                               >  8:248844/1‑113 (MQ=255)
                                                                                                                  aaCCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCg                                                               <  7:248844/113‑1 (MQ=255)
                                                                                                                   aCCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGt                          >  3:128913/1‑149 (MQ=255)
                                                                                                                        tttCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGCGTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGc                     >  1:88167/1‑149 (MQ=255)
                                                                                                                               gATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGaa               <  3:32362/148‑1 (MQ=255)
                                                                                                                               gATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAg              >  2:210225/1‑149 (MQ=255)
                                                                                                                                aTGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCAc                                      <  6:132933/124‑1 (MQ=255)
                                                                                                                                aTGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCAc                                      >  5:132933/1‑124 (MQ=255)
                                                                                                                                     gcgGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCAGCGAAGTCAGGa        >  3:187303/1‑149 (MQ=255)
                                                                                                                                              ggTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCa  <  4:157994/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGTAGCACGGAAGCGCCAGTCGGGCCGCGAAGTCAGGATGACCA  >  NC_002947/318832‑319119

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: ‑