Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 498,056 G→T 7.1% S642I (AGC→ATC)  PP_0409 → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947498,0560GT 7.1% 248.6 / 4.3 79S642I (AGC→ATC) PP_0409two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (36/37);  new base T (3/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCTGTCTGGTGC  >  NC_002947/497908‑498205
                                                                                                                                                     |                                                                                                                                                      
gagaGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGc                                                                                                                                                       <  2:2047210/150‑1 (MQ=255)
   ggcgtgcgtGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAtt                                                                                                                                                    <  2:228022/149‑2 (MQ=255)
      gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAttc                                                                                                                                                   <  1:2238898/148‑3 (MQ=255)
      gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAttc                                                                                                                                                   >  2:2238898/1‑146 (MQ=255)
      gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAg                                                                                                                                                        <  1:1785302/143‑1 (MQ=255)
      gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAg                                                                                                                                                        >  2:1785302/1‑143 (MQ=255)
      gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGaaa                                                                                                                                                 <  1:1085252/150‑1 (MQ=255)
      gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGaaa                                                                                                                                                 >  1:282597/1‑150 (MQ=255)
         cgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGc                                                                                                                                              <  2:282597/150‑1 (MQ=255)
             aaGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGa                                                                                                                                          <  2:1333937/150‑1 (MQ=255)
                agGAAAAAGCCCGTACCGCCCGGGAAGTGCATGATGAACTGGGGCAGCTGCTGACCGTGCTCAAGCTGGACGTGTCGATGTGCGAGCTGGCCTATGCCGCGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAGAAGCCTGCTCg                                                                                                                                       >  2:728492/1‑150 (MQ=255)
                  gAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGcc                                                                                                                                     >  1:631350/1‑150 (MQ=255)
                  gAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGcc                                                                                                                                     <  1:1941653/150‑1 (MQ=255)
                  gAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGcc                                                                                                                                     >  1:1191048/1‑150 (MQ=255)
                   aaaacGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGccc                                                                                                                                    <  1:829339/150‑1 (MQ=255)
                    aaaaGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCa                                                                                                                                  <  2:1938582/150‑1 (MQ=255)
                       aGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAg                                                                                                                                                        <  2:2322756/126‑1 (MQ=255)
                       aGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAg                                                                                                                                                        >  1:2322756/1‑126 (MQ=255)
                          ccGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAgttgtt                                                                                                                            >  1:2175522/1‑150 (MQ=255)
                          ccGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCCCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAgttgtt                                                                                                                            >  1:1753004/1‑150 (MQ=255)
                           cGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTc                                                                                                                           <  1:1129237/150‑1 (MQ=255)
                              ttCGCCCAGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGCGCTGGCCTATGCCGAGCTGGCCCCTGGTCTGCACGCACGACCGCCATTCATTAAACTCCTGATCTCCCAGTTGTTCCAg                                                                                                                        >  2:23780/2‑150 (MQ=255)
                               tCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGt                                                                                                                       <  1:2149476/150‑1 (MQ=255)
                                cGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGtt                                                                                                                      <  2:872684/150‑1 (MQ=255)
                                cGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGtt                                                                                                                      <  1:2352930/150‑1 (MQ=255)
                                 gCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTg                                                                                                                     <  1:387426/150‑1 (MQ=255)
                                  ccggggAGGTGCATGGGGAACTGGGGGAGATGCTGACCGGGCTCAGGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGGGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTgg                                                                                                                    <  2:757262/147‑1 (MQ=255)
                                   ccGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGt                                                                                                                   >  1:796295/1‑150 (MQ=255)
                                     gggAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGACACGCCTGATCGCCCAGTTGTTCCAGTTGGTAc                                                                                                                 <  1:1625605/150‑1 (MQ=255)
                                       gAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTAcgc                                                                                                               >  2:121968/1‑150 (MQ=255)
                                          gTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGAt                                                                                                            <  1:825574/150‑1 (MQ=255)
                                            gCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGAGATCAAAACCGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGt                                                                                                          >  1:1970757/1‑150 (MQ=255)
                                             cATGATGAACTGGGGCAGCTGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTc                                                                                                         <  1:121968/150‑1 (MQ=255)
                                             cATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGAt                                                                                                            <  1:2318188/147‑1 (MQ=255)
                                             cATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGAt                                                                                                            >  2:2318188/1‑147 (MQ=255)
                                              atgatgAACTGGGGCAGATGCTGACTGCCAACAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCg                                                                                                        >  1:1244520/1‑150 (MQ=255)
                                               tgatgaACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAgttgtt                                                                                                                            >  1:1255054/1‑129 (MQ=255)
                                               tgatgaACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAgttgtt                                                                                                                            <  2:1255054/129‑1 (MQ=255)
                                                gatgaACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGCTCGCCCAGTTGTTCCAGGTGGTACGCGAGGTCGcc                                                                                                      >  1:812492/1‑150 (MQ=255)
                                                       tGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAtt                                                                                                                                                    <  1:1320025/98‑2 (MQ=255)
                                                       tGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAtt                                                                                                                                                    >  2:1320025/1‑97 (MQ=255)
                                                        ggggCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGccc                                                                                                                                    <  1:1593087/113‑1 (MQ=255)
                                                        ggggCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGccc                                                                                                                                    >  2:1593087/1‑113 (MQ=255)
                                                          ggCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGc                                                                                            <  1:728492/150‑1 (MQ=255)
                                                           gCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCa                                                                                                                                  >  1:1451107/1‑111 (MQ=255)
                                                           gCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCa                                                                                                                                  <  2:1451107/111‑1 (MQ=255)
                                                                tGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTccgc                                                                                      <  2:631350/150‑1 (MQ=255)
                                                                 gCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTccgcc                                                                                     >  2:1944237/1‑150 (MQ=255)
                                                                   tGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGa                                                                                   >  1:255931/1‑150 (MQ=255)
                                                                     aCCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATc                                                                                 >  1:170944/1‑150 (MQ=255)
                                                                     aCCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATc                                                                                 >  2:307238/1‑150 (MQ=255)
                                                                        gTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTc                                                                              <  2:1787556/150‑1 (MQ=255)
                                                                                  gggATGTGTCGATGTGCGAGCTGGCATATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCCTCCGCCGATCCTCGATGCCGGCa                                                                    <  2:1191048/149‑1 (MQ=255)
                                                                                            gATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGc                                                          >  2:477387/1‑150 (MQ=255)
                                                                                                    ggCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGAGAGCATGAAACGCCTGCTCGCCCAGTTGTTCCAGTTCGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGact                                                  <  2:1394030/149‑3 (MQ=255)
                                                                                                         gCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGc                                               >  2:1528138/1‑148 (MQ=255)
                                                                                                         gCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGc                                               <  1:1528138/148‑1 (MQ=255)
                                                                                                           cTATGCCGAGCTGGACCCTGGCCTGAATGAGCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGgc                                           >  2:1851375/1‑150 (MQ=255)
                                                                                                            tATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCCTTGCCTCGGCCATCGAATGGCAGgcg                                          >  2:136594/1‑150 (MQ=255)
                                                                                                                  gAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGt                                    >  2:1014714/1‑150 (MQ=255)
                                                                                                                    gCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGttt                                  >  2:102242/1‑150 (MQ=255)
                                                                                                                    gCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCACTGCCTCGGCCATCGAATGGCAGGCGCGGCGttt                                  >  2:744531/1‑150 (MQ=255)
                                                                                                                     cTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTc                                 <  2:159447/150‑1 (MQ=255)
                                                                                                                       ggACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGa                               <  1:477387/150‑1 (MQ=255)
                                                                                                                         gACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAg                              >  1:484894/1‑150 (MQ=255)
                                                                                                                         gACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAg                              <  1:307238/150‑1 (MQ=255)
                                                                                                                         gACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAg                              >  1:487306/1‑150 (MQ=255)
                                                                                                                           cccTGGCCTGAATGCTCGCCTGGCGAGCATGAATCGCCTGATCGCGCAGTTGTTCCAGTTGGTAGGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCTTCGGCCATCCAATGGCAGGCGCGGCGTTTCGAGGc                            <  2:812492/150‑1 (MQ=255)
                                                                                                                            ccTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCa                                                                    <  1:2126269/109‑1 (MQ=255)
                                                                                                                            ccTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCa                                                                    >  2:2126269/1‑109 (MQ=255)
                                                                                                                            ccTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCa                           >  2:1931422/1‑150 (MQ=255)
                                                                                                                            ccTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCa                           >  2:1008966/1‑150 (MQ=255)
                                                                                                                             cTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCAc                          >  1:395552/1‑150 (MQ=255)
                                                                                                                                 ccTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTAc                      <  2:1753004/150‑1 (MQ=255)
                                                                                                                                  cTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGc                                               <  1:1712103/124‑1 (MQ=255)
                                                                                                                                  cTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGc                                               >  2:1712103/1‑124 (MQ=255)
                                                                                                                                   tGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCGCGTTGTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGc                    <  1:23780/150‑1 (MQ=255)
                                                                                                                                   tGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGc                    >  2:254775/1‑150 (MQ=255)
                                                                                                                                       tGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAAt                >  1:924268/1‑150 (MQ=255)
                                                                                                                                            gCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCt           >  1:569519/1‑150 (MQ=255)
                                                                                                                                               tGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGCGGCACGTACGCAAAGCCCCTGTc        <  2:1970757/150‑1 (MQ=255)
                                                                                                                                               tGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCTGTc        <  1:1096953/150‑1 (MQ=255)
                                                                                                                                               tGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTGCGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCTGTc        <  1:733620/150‑1 (MQ=255)
                                                                                                                                                     gCATGAACCGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGGCATCGGATTTCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCTGTCTGGTGc  <  2:1529599/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCTGTCTGGTGC  >  NC_002947/497908‑498205

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