Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,871,629 G→C 5.6% L342V (CTC→GTC) ‡ ptsP ← phosphoenolpyruvate‑dependent regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,871,6290GC 5.6% 221.1 / 3.1 71L342V (CTC→GTC) ‡ptsPphosphoenolpyruvate‑dependent regulator
Reads supporting (aligned to +/- strand):  ref base G (36/31);  new base C (2/2);  total (38/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 = 2.20e-01

CAATACCGTCGACCTTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAGTCGTCGTCCATC  >  NC_002947/5871482‑5871777
                                                                                                                                                   |                                                                                                                                                    
cAATACCGTCGACCTTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGAATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAg                                                                                                                                                      <  8:133219/148‑1 (MQ=255)
        tCGACCGTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGCTGGCGCCATGGGAGGTGCCCGAGCCGAGTACCGAGACCAGGCCGCCCAGCTTGCCTTCCGGCCCTTCACCGAGCATGGCCgg                                                                                                                                             >  6:116810/1‑149 (MQ=255)
          gACCTTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGgtg                                                                                                                                           >  5:260098/1‑149 (MQ=255)
            ccTTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGGCCATGCTCGgtgaa                                                                                                                                         <  3:271110/149‑3 (MQ=255)
            ccTTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGGCCATGCTCGgtgaa                                                                                                                                         <  1:365174/149‑3 (MQ=255)
             cTTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGgtgt                                                                                                                                          >  2:157547/1‑147 (MQ=255)
             cTTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGACGATGGCGCCATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGgtgt                                                                                                                                          >  6:384127/1‑147 (MQ=255)
                ggAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGt                                                                                                                                      <  3:58584/148‑1 (MQ=255)
                   gTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCt                                                                                                                                  <  1:69527/149‑1 (MQ=255)
                       ggCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACt                                                                                                                              >  1:225088/1‑149 (MQ=255)
                        gCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCa                                                                                                                                        <  7:119469/138‑1 (MQ=255)
                        gCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCa                                                                                                                                        >  8:119469/1‑138 (MQ=255)
                        gCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACt                                                                                                                              >  2:225976/1‑148 (MQ=255)
                        gCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACt                                                                                                                              <  3:322676/148‑1 (MQ=255)
                             tCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAg                                                                                                                        >  5:338669/1‑149 (MQ=255)
                                 ccAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATg                                                                                                                    >  2:255175/1‑149 (MQ=255)
                                  cAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATgg                                                                                                                   <  8:72477/149‑1 (MQ=255)
                                  cAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATgg                                                                                                                   >  5:35211/1‑149 (MQ=255)
                                         aTCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTc                                                                                                            >  7:267686/1‑149 (MQ=255)
                                            aCCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGc                                                                                                         <  1:225976/149‑1 (MQ=255)
                                               gTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATa                                                                                                      >  6:336334/1‑149 (MQ=255)
                                                tGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGACTTCCGTCACTTCACCGAGCTTGGCCGTTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATac                                                                                                     <  8:125392/149‑1 (MQ=255)
                                                   ggATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATacac                                                                                                   >  2:87214/1‑148 (MQ=255)
                                                           aTGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTg                                                                                          >  4:342793/1‑149 (MQ=255)
                                                               cccGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAg                                                                                      <  3:342793/149‑1 (MQ=255)
                                                                   cgcgAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAgcgcg                                                                                  <  5:384127/149‑1 (MQ=255)
                                                                    gcgAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAgcgcgc                                                                                 <  3:339623/149‑1 (MQ=255)
                                                                       aGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTc                                                                              <  2:225088/149‑1 (MQ=255)
                                                                        ggATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTcc                                                                             <  1:157547/149‑1 (MQ=255)
                                                                         gATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACGAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCt                                                                            <  8:267686/149‑1 (MQ=255)
                                                                         gATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCt                                                                            <  1:87214/149‑1 (MQ=255)
                                                                                                agccgagTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGgcgc                                                     <  7:147998/149‑1 (MQ=255)
                                                                                                   cgagTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTCCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGccc                                                   >  3:349901/1‑148 (MQ=255)
                                                                                                    gagTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCa                                                 <  6:338669/149‑1 (MQ=255)
                                                                                                     agTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAg                                                <  6:35211/149‑1 (MQ=255)
                                                                                                     agTACCGAGACCAGGCCGACCAGCTTGCCTTCCCGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAg                                                <  6:372907/149‑1 (MQ=255)
                                                                                                         ccGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTcc                                            <  3:262035/149‑1 (MQ=255)
                                                                                                          cGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCt                                           >  1:71083/1‑149 (MQ=255)
                                                                                                          cGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCt                                           >  4:88404/1‑149 (MQ=255)
                                                                                                          cGAGACCAGGCCGACCAGCTTGCCTTCCGGCACATCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCt                                           >  5:336211/1‑149 (MQ=255)
                                                                                                          cGAGACCAGGCCGAACACCTTCCCTTCCGGCACATCAAAGAGCAAGGAAGGAGTAAGTTCCTCACTGACCAGCATGGTCTTGTCGGCATACACCACCGACTACGCCAGCGCATCCTGCACGTAGGCCAAAAGGCGGCGCCCCAGGACCt                                           >  8:230238/1‑149 (MQ=255)
                                                                                                            agaCCAGGCCGACCAGCTTGCCTTCCGTCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTg                                         <  1:380214/149‑1 (MQ=255)
                                                                                                             gaCCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGa                                        <  1:255175/149‑1 (MQ=255)
                                                                                                                cAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGTGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCtt                                          >  1:185191/1‑144 (MQ=255)
                                                                                                                cAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGTGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCtt                                          <  2:185191/144‑1 (MQ=255)
                                                                                                                cAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGt                                     >  8:327396/1‑149 (MQ=255)
                                                                                                                  ggCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCg                                   >  5:310174/1‑149 (MQ=255)
                                                                                                                  ggCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCg                                   >  3:168440/1‑149 (MQ=255)
                                                                                                                   gCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCtt                                          >  6:134121/1‑141 (MQ=255)
                                                                                                                   gCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCtt                                          <  5:134121/141‑1 (MQ=255)
                                                                                                                   gCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCg                                   >  3:289402/1‑148 (MQ=255)
                                                                                                                     cGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCaaa                                                              <  1:350451/119‑1 (MQ=255)
                                                                                                                     cGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCaaa                                                              >  2:350451/1‑119 (MQ=255)
                                                                                                                     cGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTc                                    >  4:313880/1‑145 (MQ=255)
                                                                                                                     cGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTc                                    <  3:313880/145‑1 (MQ=255)
                                                                                                                               gCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCAc                      >  1:349900/1‑149 (MQ=255)
                                                                                                                                 cTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCTGCATACACCAGCGACTGCGa                                                                                       <  2:151801/82‑1 (MQ=255)
                                                                                                                                 cTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGa                                                                                       >  1:151801/1‑82 (MQ=255)
                                                                                                                                 cTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACg                     >  1:290923/1‑148 (MQ=255)
                                                                                                                                  ttCCGGCACTTCACCGAGCATGGCCGGTGTCGGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACg                     >  3:181078/1‑147 (MQ=255)
                                                                                                                                    ccGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAgg                 >  2:162745/1‑149 (MQ=255)
                                                                                                                                    ccGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAg                  <  2:290923/148‑1 (MQ=255)
                                                                                                                                       gCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAg              <  5:112260/149‑1 (MQ=255)
                                                                                                                                          tgacACCGGCCATGCTCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAgtcg           >  2:366477/4‑149 (MQ=255)
                                                                                                                                          cTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAgtcg           <  2:233754/149‑1 (MQ=255)
                                                                                                                                           ttCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAgtcgt          >  8:318402/1‑149 (MQ=255)
                                                                                                                                           gacACCGGCCATGCTCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAgtcgt          >  6:78124/3‑149 (MQ=255)
                                                                                                                                             cACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAgtcgtcg        >  7:60826/1‑149 (MQ=255)
                                                                                                                                              aCCGAGCATGGCCGGTGTCAGTCCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAgtcgtcgt       <  4:349901/149‑1 (MQ=255)
                                                                                                                                                cGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAgtcgtcgtc      >  6:5965/1‑148 (MQ=255)
                                                                                                                                                 gAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAGTCGTCGTCCa    >  5:285550/1‑149 (MQ=255)
                                                                                                                                                  aGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAGTCGTCGTCCAt   <  2:349900/149‑1 (MQ=255)
                                                                                                                                                   gCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAGTCGTCGTCCATc  <  6:336211/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAATACCGTCGACCTTGGAGTATGGCAGGTCGACCAGGCCCATCACCGTGGGGATGCCCATGGCCCGCGCGAGGATGGCGACATGGGAGTTGCCCGAGCCGAGTACCGAGACCAGGCCGACCAGCTTGCCTTCCGGCACTTCACCGAGCATGGCCGGTGTCAGTTCCTCACTGACCAGGATGGTGTTGTCGGCATACACCAGCGACTGCGAGCGCGCTTCCTGCAGGTAGGCCAAAAGGCGGCGCCCCAGGTCCTTGACGTCGGAGGCACGCTCACGCAGGTAGTCGTCGTCCATC  >  NC_002947/5871482‑5871777

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: