Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,513,986 A→T 8.9% Y204F (TAC→TTC) ‡ urtE → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,513,9860AT 8.9% 245.5 / 10.5 87Y204F (TAC→TTC) ‡urtEABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (43/36);  new base T (5/3);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CTGATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTCGCGATG  >  NC_002947/5513838‑5514127
                                                                                                                                                    |                                                                                                                                              
ctgATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTa                                                                                                                                                <  7:287534/149‑1 (MQ=255)
 tgATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTAc                                                                                                                                               <  7:304942/149‑1 (MQ=255)
   aTTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCt                                                                                                                                             >  2:98167/1‑149 (MQ=255)
       tcGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGt                                                                                                                                          >  2:236745/1‑148 (MQ=255)
              gCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGGTACTGGTCagccagc                                                                                                                                  <  4:334423/149‑7 (MQ=255)
                 gaccgaAGGCATCCTGCCCTCGGTGATCTAGGAAATCGGCGCGGTGATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTAGGACTTTGCCGAGGCGCTGGCTGACCTGTACCTGGTCATGGCCCGTg                                                                                                                               <  6:299168/149‑1 (MQ=255)
                  accgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGGTACTGGTCag                                                                                                                                       >  4:289595/1‑139 (MQ=255)
                  accgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGGTACTGGTCag                                                                                                                                       <  3:289595/140‑2 (MQ=255)
                    cgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                             >  6:366548/1‑148 (MQ=255)
                    cgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                             >  8:187750/1‑148 (MQ=255)
                     gaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTgg                                                                                                                              >  7:294030/1‑146 (MQ=255)
                     gaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTgg                                                                                                                              <  8:294030/146‑1 (MQ=255)
                      aaGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGa                                                                                                                           >  1:140716/1‑148 (MQ=255)
                       aGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGaaa                                                                                                                         >  2:233890/1‑149 (MQ=255)
                                    tCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCCTCCAGCGCgg                                                                                                            >  7:269680/1‑149 (MQ=255)
                                      ggTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGGTACTGGTCagccag                                                                                                                                   >  6:104625/1‑119 (MQ=255)
                                           tCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCgcgg                                                                                                      <  2:140716/148‑1 (MQ=255)
                                              aGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                             <  6:19689/122‑1 (MQ=255)
                                              aGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGa                                                                                                   >  1:276170/1‑148 (MQ=255)
                                              aGGAAATCGGCGCGGTCATCCGCCGCCTGACCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGa                                                                                                   >  3:172989/1‑148 (MQ=255)
                                               ggAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGaaa                                                                                                 >  7:318064/1‑149 (MQ=255)
                                                   aTCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACAt                                                                                             <  4:213817/149‑1 (MQ=255)
                                                    tCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATg                                                                                            >  5:34588/1‑149 (MQ=255)
                                                       gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTAcc                                                                                                                                              >  2:321699/1‑96 (MQ=255)
                                                       gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTAcc                                                                                                                                              >  4:256677/1‑96 (MQ=255)
                                                       gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTAcc                                                                                                                                              <  3:256677/96‑1 (MQ=255)
                                                       gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTAcc                                                                                                                                              <  1:321699/96‑1 (MQ=255)
                                                       gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGAACCTCGTCATGGCCCGTGGCGAAATCATCCAGCGCGGACGCGGCGAAAACaaccat                                                                                         >  8:165615/1‑144 (MQ=255)
                                                         gcgGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGc                                                                                       <  7:187750/149‑1 (MQ=255)
                                                             tCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGaa                                                                                   <  4:330735/149‑1 (MQ=255)
                                                               aTCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGc                                                                                           <  1:186372/139‑1 (MQ=255)
                                                               aTCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGc                                                                                           >  2:186372/1‑139 (MQ=255)
                                                               aTCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAgg                                                                                 <  2:276170/149‑1 (MQ=255)
                                                               aTCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAg                                                                                  <  6:34588/148‑1 (MQ=255)
                                                               aTCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAg                                                                                  <  5:197765/148‑1 (MQ=255)
                                                               aTCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAg                                                                                  <  1:98167/148‑1 (MQ=255)
                                                               aTCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGgt                                                                                <  5:142983/149‑1 (MQ=255)
                                                                tCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGgt                                                                                <  4:100644/149‑1 (MQ=255)
                                                                tCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGgt                                                                                <  4:76317/149‑1 (MQ=255)
                                                                 ccgccgCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGgtg                                                                               >  5:319942/1‑149 (MQ=255)
                                                                  cgccgcCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGgtg                                                                               >  4:314906/1‑148 (MQ=255)
                                                                    ccgccTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTAc                                                                            >  8:58794/1‑149 (MQ=255)
                                                                             cGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTgg                                                                  <  3:331741/149‑1 (MQ=255)
                                                                                   tGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCa                                                            <  2:159060/149‑1 (MQ=255)
                                                                                      cGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCAtt                                                          >  2:121775/1‑148 (MQ=255)
                                                                                        aCATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGGCCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTaa                                                       >  3:221536/1‑149 (MQ=255)
                                                                                           tGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCt                                                    >  8:273840/1‑149 (MQ=255)
                                                                                               gATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                <  5:366548/149‑1 (MQ=255)
                                                                                               gATTTTGCTGGTGGAGCAGTTCTACGACTTTGACGAGTCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGACGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGc                                                  >  1:38747/1‑147 (MQ=255)
                                                                                                aTTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAg                                               >  7:282956/1‑149 (MQ=255)
                                                                                                        ggtggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGAGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACAATTTAATCTGAGAGCAAGGAcc                                       >  8:354961/1‑149 (MQ=255)
                                                                                                         gtggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCCGGGCTGGTAACCATTTAATCTGCGAGCAACGAc                                        <  1:233890/148‑1 (MQ=255)
                                                                                                              gCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCg                                  <  5:58971/148‑1 (MQ=255)
                                                                                                              gCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCg                                  >  6:58971/1‑148 (MQ=255)
                                                                                                                aGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCgaga                               <  3:314906/149‑1 (MQ=255)
                                                                                                                  ttCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACg                             >  8:310708/1‑149 (MQ=255)
                                                                                                                    cTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAg                                                                                  <  2:305368/95‑1 (MQ=255)
                                                                                                                    cTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAg                                                                                  >  1:305368/1‑95 (MQ=255)
                                                                                                                     tACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCtt                                     >  6:186073/1‑138 (MQ=255)
                                                                                                                     tACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCtt                                     <  5:186073/138‑1 (MQ=255)
                                                                                                                     tACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTg                          >  8:164581/1‑149 (MQ=255)
                                                                                                                     tACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTg                          >  7:292431/1‑149 (MQ=255)
                                                                                                                      aCGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCACGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGt                         >  2:303564/1‑149 (MQ=255)
                                                                                                                      aCGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGt                         >  5:231414/1‑149 (MQ=255)
                                                                                                                      aCGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGt                         >  7:331663/1‑149 (MQ=255)
                                                                                                                       cGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTg                          <  1:236745/147‑1 (MQ=255)
                                                                                                                       cGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTc                        <  6:319942/149‑1 (MQ=255)
                                                                                                                        gACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCa                       >  3:345812/1‑149 (MQ=255)
                                                                                                                         aCTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCAt                      <  4:172989/149‑1 (MQ=255)
                                                                                                                           tttGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACAATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCAt                      >  5:250345/1‑147 (MQ=255)
                                                                                                                           tttGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACAATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCAt                      >  7:91092/1‑147 (MQ=255)
                                                                                                                              gCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCg                                                 <  2:171896/117‑1 (MQ=255)
                                                                                                                              gCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCg                                                 >  1:171896/1‑117 (MQ=255)
                                                                                                                                cGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTAc               >  2:310586/1‑149 (MQ=255)
                                                                                                                                 gAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACg              <  7:58794/149‑1 (MQ=255)
                                                                                                                                      gctggctgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACg                             <  6:240915/129‑1 (MQ=255)
                                                                                                                                      gctggctgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACg                             >  5:240915/1‑129 (MQ=255)
                                                                                                                                       ctggctgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTc        <  7:273840/149‑1 (MQ=255)
                                                                                                                                       ctggctgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTc        <  6:231414/149‑1 (MQ=255)
                                                                                                                                         ccatgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTaa                                                               >  4:184081/4‑92 (MQ=255)
                                                                                                                                         ccatgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTaa                                                               <  3:184081/89‑1 (MQ=255)
                                                                                                                                         ccatgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGGCCTTGCGAGACGCTGTCATTCGCTACGAAATTcgc      >  7:129094/4‑149 (MQ=255)
                                                                                                                                         ccatgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTcgc      >  4:191330/4‑149 (MQ=255)
                                                                                                                                          catgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGTAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTcgcg     >  3:248256/3‑149 (MQ=255)
                                                                                                                                          catgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGCCCTTGCGAGACGCTGTCATGAGCTACGAAATTcgc      >  3:253067/3‑148 (MQ=255)
                                                                                                                                           ctgACCCGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAACACATGCAAGCCGAAGGGGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTCGCGa    <  8:318064/149‑1 (MQ=255)
                                                                                                                                            tgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTCGCGAt   >  7:68843/1‑149 (MQ=255)
                                                                                                                                              aCCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTCGCGATg  <  1:310586/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTGATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACGCTGTCATGAGCTACGAAATTCGCGATG  >  NC_002947/5513838‑5514127

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: