Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 498,050 T→A 6.3% L640Q (CTG→CAG)  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,0500TA 6.3% 243.2 / 6.7 81L640Q (CTG→CAG) PP_0409two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base A (2/3);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

TGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCTGTC  >  NC_002947/497906‑498199
                                                                                                                                                 |                                                                                                                                                      
tGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGACGTATCTAGGTGCGAGCTGGCATCTGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGa                                                                                                                                                   >  1:1207282/1‑150 (MQ=255)
  gagaGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGc                                                                                                                                                 <  2:2047210/150‑1 (MQ=255)
     ggcgtgcgtGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAtt                                                                                                                                              <  2:228022/149‑2 (MQ=255)
        gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAttc                                                                                                                                             >  2:2238898/1‑146 (MQ=255)
        gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAttc                                                                                                                                             <  1:2238898/148‑3 (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:1191048/1‑150 (MQ=255)
                    gAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGcc                                                                                                                               <  1:1941653/150‑1 (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                                                                                                                      <  2:1255054/129‑1 (MQ=255)
                                                 tgatgaACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAgttgtt                                                                                                                      >  1:1255054/1‑129 (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                                                                           >  2:307238/1‑150 (MQ=255)
                                                                       aCCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATc                                                                           >  1:170944/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                                         <  1:1528138/148‑1 (MQ=255)
                                                                                                           gCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGc                                         >  2:1528138/1‑148 (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:484894/1‑150 (MQ=255)
                                                                                                                           gACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAg                        <  1:307238/150‑1 (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:1008966/1‑150 (MQ=255)
                                                                                                                              ccTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCa                     >  2:1931422/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)
                                                                                                                                                 |                                                                                                                                                      
TGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCTGTC  >  NC_002947/497906‑498199

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: