Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 338,068 A→C 5.6% S167R (AGT→AGG)  PP_0280 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947338,0680AC 5.6% 281.6 / 5.1 89S167R (AGT→AGG) PP_0280amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (36/48);  new base C (3/2);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGC  >  NC_002947/337920‑338211
                                                                                                                                                    |                                                                                                                                               
gcagcTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCa                                                                                                                                                 >  5:86650/1‑149 (MQ=255)
    cTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTgg                                                                                                                                             >  5:75155/1‑149 (MQ=255)
        aaGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGcca                                                                                                                                          <  6:75155/148‑2 (MQ=255)
         aGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGccagt                                                                                                                                        <  7:337551/149‑4 (MQ=255)
          gCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCtt                                                                                                                                         >  8:234970/1‑147 (MQ=255)
          gCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                       <  7:175171/149‑1 (MQ=255)
           ccTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAg                                                                                                                                      >  5:278425/1‑149 (MQ=255)
                cACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCa                                                                                                                                 <  4:286154/149‑1 (MQ=255)
                 aCCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAg                                                                                                                                >  6:127591/1‑149 (MQ=255)
                  ccAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAATATTTCCACCGGCAGGTAGGTGCTGGCTATGATTGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCGTCAgg                                                                                                                               <  6:342227/149‑1 (MQ=255)
                   cAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGa                                                                                                                              <  4:346587/149‑1 (MQ=255)
                      cATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCa                                                                                                                           <  6:23127/149‑1 (MQ=255)
                        tGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACt                                                                                                                         >  8:182185/1‑149 (MQ=255)
                            cGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                       <  5:260414/131‑1 (MQ=255)
                            cGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                       >  6:260414/1‑131 (MQ=255)
                             gagaTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGTCTTGAGCATCAGGATCACTTCGtt                                                                                                                    >  1:307558/1‑149 (MQ=255)
                             gagaTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGtt                                                                                                                    <  2:59921/149‑1 (MQ=255)
                              agaTCACCAGGGAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTg                                                                                                                   <  8:271693/149‑1 (MQ=255)
                                aTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCt                                                                                                                 <  5:127591/149‑1 (MQ=255)
                                   aCCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                       <  6:193316/124‑1 (MQ=255)
                                   aCCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                       >  5:193316/1‑124 (MQ=255)
                                   aCCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTa                                                                                                              >  3:169898/1‑149 (MQ=255)
                                        gTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                         <  4:258982/149‑1 (MQ=255)
                                        gTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                         >  4:234198/1‑149 (MQ=255)
                                        gTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                         >  6:117019/1‑149 (MQ=255)
                                           gAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGc                                                                                                       <  3:353990/148‑1 (MQ=255)
                                              cacGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGACCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGc                                                                                                   <  1:311318/149‑1 (MQ=255)
                                               acGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGcc                                                                                                  >  4:122630/1‑149 (MQ=255)
                                                      cggcgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgcg                                                                                            <  6:224182/148‑1 (MQ=255)
                                                      cggcgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgc                                                                                             <  6:278425/147‑1 (MQ=255)
                                                         cgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCa                                                                                        <  5:117019/149‑1 (MQ=255)
                                                          gaagaaGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                       >  3:21317/1‑101 (MQ=255)
                                                          gaagaaGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                       <  4:21317/101‑1 (MQ=255)
                                                           aagaagATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGt                                                                                                                     >  5:304989/1‑118 (MQ=255)
                                                           aagaagATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGt                                                                                                                     <  6:304989/118‑1 (MQ=255)
                                                           aagaagATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgc                                                                                      >  2:238475/1‑149 (MQ=255)
                                                                   ttCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGt                                                                                                               <  6:321916/116‑1 (MQ=255)
                                                                   ttCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGt                                                                                                               >  5:321916/1‑116 (MQ=255)
                                                                   ttCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTgg                                                                                >  2:166435/1‑147 (MQ=255)
                                                                      cACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGgcagca                                                                           >  5:30472/1‑149 (MQ=255)
                                                                      cACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGgcagca                                                                           >  1:302568/1‑149 (MQ=255)
                                                                      cACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGgcagc                                                                            <  2:358378/148‑1 (MQ=255)
                                                                      cACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGgcagc                                                                            <  7:182185/148‑1 (MQ=255)
                                                                        ccGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCatg                                                                         <  3:234198/149‑1 (MQ=255)
                                                                            cAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGt                                                                     <  6:30472/149‑1 (MQ=255)
                                                                             aggtaggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTa                                                                    <  5:328316/149‑1 (MQ=255)
                                                                             aggtaggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCAGCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTa                                                                    <  7:234970/149‑1 (MQ=255)
                                                                              ggtaggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAg                                                                   <  3:316458/149‑1 (MQ=255)
                                                                               agaggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATc                                                                                                                                  <  4:211184/83‑1 (MQ=255)
                                                                                taggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGaa                                                                 <  1:192468/149‑1 (MQ=255)
                                                                                  ggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAAc                                                                <  7:39812/148‑1 (MQ=255)
                                                                                  ggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAAc                                                                <  3:330607/148‑1 (MQ=255)
                                                                                    tGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGc                                                             >  5:97492/1‑149 (MQ=255)
                                                                                       gggCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCa                                                                                        <  2:354595/119‑1 (MQ=255)
                                                                                       gggCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCa                                                                                        >  1:354595/1‑119 (MQ=255)
                                                                                         gCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGtttt                                                        <  3:122630/149‑1 (MQ=255)
                                                                                         gCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGtttt                                                        <  1:53891/149‑1 (MQ=255)
                                                                                         gCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGtttt                                                        <  6:86650/149‑1 (MQ=255)
                                                                                          cGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTg                                                       <  1:166435/149‑1 (MQ=255)
                                                                                          cGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTg                                                       >  5:16672/1‑149 (MQ=255)
                                                                                            atgatgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGcc                                                     <  2:302568/149‑1 (MQ=255)
                                                                                               atgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAcg                                                  >  5:171365/1‑149 (MQ=255)
                                                                                                tgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAcgc                                                 >  6:108589/1‑149 (MQ=255)
                                                                                                tgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAcgc                                                 >  8:192091/1‑149 (MQ=255)
                                                                                                tgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAcgc                                                 <  4:169898/149‑1 (MQ=255)
                                                                                                  gTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGTCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGa                                               <  2:307558/149‑1 (MQ=255)
                                                                                                        gCCATGCCGGTCAGTTCCAGCAGGGTAACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATAACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGAGGCAGCATGAAGTAGAACAGCGTTTTGCCACGCGACATGcc                                         >  7:297308/1‑149 (MQ=255)
                                                                                                          cATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGcc                                         >  4:202301/1‑147 (MQ=255)
                                                                                                           aTGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAg                                      >  2:269525/1‑149 (MQ=255)
                                                                                                                ggTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGccc                                 >  3:244501/1‑149 (MQ=255)
                                                                                                                 gTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCg                                >  7:187443/1‑149 (MQ=255)
                                                                                                                  tCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGt                               <  1:269525/149‑1 (MQ=255)
                                                                                                                  tCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGAGGCAGCATGATGTAGAGCAGCGtgt                                                         >  8:121124/1‑121 (MQ=255)
                                                                                                                    aGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGc                             >  6:124624/1‑149 (MQ=255)
                                                                                                                     gTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGcc                            >  4:100669/1‑149 (MQ=255)
                                                                                                                         cagcagGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCt                         <  1:176398/148‑1 (MQ=255)
                                                                                                                         cagcagGGTCACGGGACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCt                         <  5:51618/148‑1 (MQ=255)
                                                                                                                             agGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAAt                    <  3:100669/149‑1 (MQ=255)
                                                                                                                               ggTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAAtct                  >  7:256675/1‑149 (MQ=255)
                                                                                                                                      gTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCttt           <  6:97492/149‑1 (MQ=255)
                                                                                                                                      ctcaaggCCAGTGCCCTGGCCTTGCGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCCCGCGACATACCCAGCGCCAGTGCCGCTTCCATCTCCCCtg            >  1:364461/6‑147 (MQ=255)
                                                                                                                                       tACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTa          <  3:202301/149‑1 (MQ=255)
                                                                                                                                         cTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAgg        <  5:108589/149‑1 (MQ=255)
                                                                                                                                         cTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAgg        <  6:171365/149‑1 (MQ=255)
                                                                                                                                         cTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAg         <  8:187443/148‑1 (MQ=255)
                                                                                                                                           ggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGa      >  3:339554/1‑149 (MQ=255)
                                                                                                                                            gCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGAt     >  3:223737/1‑149 (MQ=255)
                                                                                                                                            gCCAGGGCACTGGCCTTGAGCATCAGGATCAATTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATCCCCACCGCCCGTGCCGCTTAAATCTCACATTTAGTGAt     >  8:266405/1‑149 (MQ=255)
                                                                                                                                              cAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCg   >  5:253639/1‑149 (MQ=255)
                                                                                                                                               aGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGc  <  5:124624/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGC  >  NC_002947/337920‑338211

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