Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,513,985 T→G 5.8% Y204D (TAC→GAC) ‡ 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,9850TG 5.8% 214.5 / 4.6 69Y204D (TAC→GAC) ‡urtEABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base G (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

ATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTCGCGATGCCC  >  NC_002947/5513841‑5514130
                                                                                                                                                   |                                                                                                                                                  
aTTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGCTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGCCCGAGCCTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTATGCCGAGGCGCTGGCTGACCAGTACCt                                                                                                                                                >  1:400636/1‑149 (MQ=255)
 ttCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTg                                                                                                                                               <  4:312037/149‑1 (MQ=255)
      gacgaGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTc                                                                                                                                            <  5:153830/147‑1 (MQ=255)
         gaGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATTACGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGGTACTGGTCagc                                                                                                                                         <  1:325116/149‑3 (MQ=255)
         gaGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGGTACTGGTCagcca                                                                                                                                       <  2:230783/149‑5 (MQ=255)
              gaccgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTg                                                                                                                                  >  4:314628/1‑149 (MQ=255)
               accgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTgg                                                                                                                                 <  7:359065/149‑1 (MQ=255)
                ccgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                                <  6:400629/149‑1 (MQ=255)
                 cgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                                >  4:120816/1‑148 (MQ=255)
                 cgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                                >  2:223837/1‑148 (MQ=255)
                  gaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCg                                                                                                                               >  3:61008/1‑148 (MQ=255)
                     ggCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAat                                                                                                                           >  5:88425/1‑149 (MQ=255)
                        aTCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAatcat                                                                                                                        <  1:68384/149‑1 (MQ=255)
                        aTCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAatcat                                                                                                                        <  1:223837/149‑1 (MQ=255)
                          ccAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAatcatc                                                                                                                       >  7:356553/1‑148 (MQ=255)
                                   ggTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTc                                                                                                             >  5:326252/1‑149 (MQ=255)
                                          aaGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCgcggcg                                                                                                       >  3:230383/1‑148 (MQ=255)
                                            ggAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGATGCGCTGGCTGACCAGTACATGGTCATGGCCCGTGGCGAAATCATCAAGCGCGGTCGCGGCGaa                                                                                                     >  8:160943/1‑148 (MQ=255)
                                            ggAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGaa                                                                                                     >  6:341099/1‑148 (MQ=255)
                                            ggAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGAACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGaa                                                                                                     >  4:154454/1‑148 (MQ=255)
                                                aTCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACAt                                                                                                <  8:371019/149‑1 (MQ=255)
                                                    gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                            <  7:161053/149‑1 (MQ=255)
                                                    gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                            <  7:155620/149‑1 (MQ=255)
                                                    gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                            <  3:185720/149‑1 (MQ=255)
                                                    gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                            <  2:400636/149‑1 (MQ=255)
                                                            aTCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAg                                                                                     <  5:341099/148‑1 (MQ=255)
                                                            aTCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAg                                                                                     <  1:178794/148‑1 (MQ=255)
                                                             tCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGgt                                                                                   <  6:326252/149‑1 (MQ=255)
                                                                   gccTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTAcgc                                                                             >  1:142668/1‑149 (MQ=255)
                                                                    ccTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTAagcg                                                                            >  8:318592/1‑149 (MQ=255)
                                                                     cTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCgg                                                                           <  2:315311/149‑1 (MQ=255)
                                                                     cTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCgg                                                                           <  8:356553/149‑1 (MQ=255)
                                                                        gCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCt                                                                        <  8:51449/149‑1 (MQ=255)
                                                                         ccGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTg                                                                       >  1:369103/1‑149 (MQ=255)
                                                                          cGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGTGGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCggg                                                                          >  1:243063/1‑146 (MQ=255)
                                                                          cGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCggg                                                                          <  2:243063/145‑1 (MQ=255)
                                                                                  gCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTTGCGAAATCATCCAGCGCGGTCGCGGCGAATACATGCAAGCCGAAGGTGTACGGGGTCTGGTAACCAtt                                                              <  4:61008/149‑1 (MQ=255)
                                                                                    gACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTa                                                            <  3:120816/149‑1 (MQ=255)
                                                                                           ggCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGTCGAAGGTGTACGCGGGCTGGTaa                                                                   <  1:209450/137‑1 (MQ=255)
                                                                                           ggCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGTCGAAGGTGTACGCGGGCTGGTaa                                                                   >  2:209450/1‑137 (MQ=255)
                                                                                              gATGTTGCGGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGGGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                    <  6:344468/149‑1 (MQ=255)
                                                                                                     cTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGa                                             <  3:314628/149‑1 (MQ=255)
                                                                                                      tggt‑ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGc                                                      >  5:106592/1‑139 (MQ=255)
                                                                                                      tggt‑ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGc                                                      <  6:106592/139‑1 (MQ=255)
                                                                                                      tggt‑ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGAc                                            >  8:272295/1‑149 (MQ=255)
                                                                                                        gt‑ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGc                                                      <  1:53744/137‑1 (MQ=255)
                                                                                                        gt‑ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGc                                                      >  2:53744/1‑137 (MQ=255)
                                                                                                           gtaGCAGTTCTACGACTTTGCCGAGGCGCTGTCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTaa                                                                   <  4:54414/120‑1 (MQ=255)
                                                                                                           ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGCCCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTaa                                                                   >  3:54414/1‑122 (MQ=255)
                                                                                                           ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                    >  3:344440/1‑137 (MQ=255)
                                                                                                           ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                    <  4:344440/137‑1 (MQ=255)
                                                                                                            gAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGc                                       >  3:193401/1‑149 (MQ=255)
                                                                                                               cAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCga                                     <  2:142668/148‑1 (MQ=255)
                                                                                                                  ttCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑Cg                                >  2:326973/1‑149 (MQ=255)
                                                                                                                     tACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTg                             >  1:107601/1‑149 (MQ=255)
                                                                                                                       cGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCAtt                                                              <  3:246439/115‑1 (MQ=255)
                                                                                                                         aCTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCAt                         <  1:326973/149‑1 (MQ=255)
                                                                                                                             tGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGa                       <  4:230383/147‑1 (MQ=255)
                                                                                                                                    gcgcTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGaaa              <  6:88425/149‑1 (MQ=255)
                                                                                                                                          gctgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTcgc         <  2:107601/148‑1 (MQ=255)
                                                                                                                                          gctgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTcgc         >  7:9282/1‑148 (MQ=255)
                                                                                                                                          catgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTCCGAAATTcgc         >  2:339031/3‑148 (MQ=255)
                                                                                                                                          catgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTcgc         >  2:60613/3‑148 (MQ=255)
                                                                                                                                            tgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTCGCGAt      >  2:309770/1‑149 (MQ=255)
                                                                                                                                               ccAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGACCGCTGTCATGAGCTACGAAATTCGCGAt      <  1:230648/147‑1 (MQ=255)
                                                                                                                                               ccAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTCGCGAt      >  2:230648/1‑146 (MQ=255)
                                                                                                                                                cAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTCGCGATGccc  >  6:302463/1‑149 (MQ=255)
                                                                                                                                                  gTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                    <  8:114578/98‑1 (MQ=255)
                                                                                                                                                  gTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                    >  7:114578/1‑98 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTCGCGATGCCC  >  NC_002947/5513841‑5514130

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