Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,513,987 C→A 5.7% Y204* (TAC→TAA) ‡ 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,9870CA 5.7% 225.1 / 3.6 70Y204* (TAC→TAA) ‡urtEABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (30/36);  new base A (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

ATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTCGCGATGCCCTGC  >  NC_002947/5513841‑5514133
                                                                                                                                                     |                                                                                                                                                   
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                                                                                               <  3:185720/149‑1 (MQ=255)
                                                    gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                               <  2:400636/149‑1 (MQ=255)
                                                    gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                               <  7:155620/149‑1 (MQ=255)
                                                    gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                               <  7:161053/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                                                                              <  8:356553/149‑1 (MQ=255)
                                                                     cTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCgg                                                                              <  2:315311/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                                                         >  2:53744/1‑137 (MQ=255)
                                                                                                        gt‑ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGc                                                         <  1:53744/137‑1 (MQ=255)
                                                                                                           gtaGCAGTTCTACGACTTTGCCGAGGCGCTGTCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTaa                                                                      <  4:54414/120‑1 (MQ=255)
                                                                                                           ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGCCCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTaa                                                                      >  3:54414/1‑122 (MQ=255)
                                                                                                           ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                       <  4:344440/137‑1 (MQ=255)
                                                                                                           ggAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                       >  3:344440/1‑137 (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            >  7:9282/1‑148 (MQ=255)
                                                                                                                                          gctgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTcgc            <  2:107601/148‑1 (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)
                                                                                                                                                    aCCTGGTCATGGCCCGTTGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACTACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTCGCGATGCCCTGc  <  7:160943/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACAT‑‑GGCGATTTTGCTGGT‑GGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTGCGAGA‑CGCTGTCATGAGCTACGAAATTCGCGATGCCCTGC  >  NC_002947/5513841‑5514133

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