Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 498,044 A→C 7.6% D638A (GAT→GCT)  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,0440AC 7.6% 215.0 / 7.4 80D638A (GAT→GCT) PP_0409two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base C (3/4);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CTCGGCCCACCTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAAT  >  NC_002947/497895‑498191
                                                                                                                                                      |                                                                                                                                                    
ctcGGCCCACTTTTAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGa                                                                                                                                                      <  2:508888/150‑1 (MQ=255)
ctcGGCCCACCTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTCGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGa                                                                                                                                                      <  1:1455894/150‑1 (MQ=255)
           tGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGACGTATCTAGGTGCGAGCTGGCATCTGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGa                                                                                                                                           >  1:1207282/1‑150 (MQ=255)
             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                                                                                                                                          >  2:1785302/1‑143 (MQ=255)
                   gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAg                                                                                                                                          <  1:1785302/143‑1 (MQ=255)
                   gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGaaa                                                                                                                                   >  1:282597/1‑150 (MQ=255)
                   gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGaaa                                                                                                                                   <  1:1085252/150‑1 (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:1941653/150‑1 (MQ=255)
                               gAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGcc                                                                                                                       >  1:1191048/1‑150 (MQ=255)
                               gAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGcc                                                                                                                       >  1:631350/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                                                                                              >  2:2318188/1‑147 (MQ=255)
                                                          cATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGAt                                                                                              <  1:2318188/147‑1 (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                                                                                                                    <  2:1451107/111‑1 (MQ=255)
                                                                        gCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCa                                                                                                                    >  1:1451107/1‑111 (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:487306/1‑150 (MQ=255)
                                                                                                                                      gACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAg                <  1:307238/150‑1 (MQ=255)
                                                                                                                                      gACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAg                >  1:484894/1‑150 (MQ=255)
                                                                                                                                        cccTGGCCTGAATGCTCGCCTGGCGAGCATGAATCGCCTGATCGCGCAGTTGTTCCAGTTGGTAGGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCTTCGGCCATCCAATGGCAGGCGCGGCGTTTCGAGGc              <  2:812492/150‑1 (MQ=255)
                                                                                                                                         ccTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCa                                                      >  2:2126269/1‑109 (MQ=255)
                                                                                                                                         ccTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCa                                                      <  1:2126269/109‑1 (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                                 >  2:1712103/1‑124 (MQ=255)
                                                                                                                                               cTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGc                                 <  1:1712103/124‑1 (MQ=255)
                                                                                                                                                tGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCGCGTTGTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGc      <  1:23780/150‑1 (MQ=255)
                                                                                                                                                tGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGc      >  2:254775/1‑150 (MQ=255)
                                                                                                                                                    tGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAAt  >  1:924268/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTCGGCCCACCTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAAT  >  NC_002947/497895‑498191

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: