Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,513,985 T→G 7.9% 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 7.9% 186.9 / 7.3 63Y204D (TAC→GAC) ‡urtEABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base G (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGCCCTG  >  NC_002947/5513840‑5514132
                                                                                                                                                 |                                                                                                                                                    
gATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTAcc                                                                                                                                                   >  8:59446/1‑149 (MQ=255)
   tctcGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTgg                                                                                                                                                >  6:368608/1‑149 (MQ=255)
    ctcGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTgg                                                                                                                                                <  6:219783/148‑1 (MQ=255)
     tcGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCCATATTCCTCGAGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCAGGt                                                                                                                                               >  4:318262/1‑148 (MQ=255)
       gacgaGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCa                                                                                                                                             <  1:351733/148‑1 (MQ=255)
            gCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGccc                                                                                                                                       <  1:221763/149‑1 (MQ=255)
                 ccgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                                  <  2:127167/149‑1 (MQ=255)
                  cgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                                  >  6:30872/1‑148 (MQ=255)
                  cgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                                  >  3:351166/1‑148 (MQ=255)
                  cgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                                  >  5:245790/1‑148 (MQ=255)
                        cATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGATCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAatca                                                                                                                           <  2:207050/149‑1 (MQ=255)
                           ccAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAatcatc                                                                                                                         >  6:192525/1‑148 (MQ=255)
                                ccTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAgcgc                                                                                                                   <  4:14666/149‑1 (MQ=255)
                                        aTCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCgcgg                                                                                                           <  3:178830/149‑1 (MQ=255)
                                        aTCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGATCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCgcgg                                                                                                           >  1:355357/1‑149 (MQ=255)
                                                     gcgcGGTCATCCGCCGCCTGGCCGAGCGTGGAGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                              <  4:46849/149‑1 (MQ=255)
                                                          gTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGa                                                                                         <  1:233820/149‑1 (MQ=255)
                                                             aTCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTTCGACTTTGCCGTGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAg                                                                                       <  3:318262/148‑1 (MQ=255)
                                                               ccgccgCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCCTCCCGCGCGGTCTCGGCGAAAACATGCAAGCCGAAGgtg                                                                                    >  7:156168/1‑149 (MQ=255)
                                                               ccgccgCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGgtg                                                                                    >  5:211916/1‑149 (MQ=255)
                                                                 gccgccTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                              >  2:76740/1‑137 (MQ=255)
                                                                 gccgccTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCaa                                                                                              <  1:76740/137‑1 (MQ=255)
                                                                    gccTGGCCGAGCGTGGCGACATGACAATTTTGCAGGTAGAGCAATTCTACAACTTTGCCTAAGCGCTGACTGACAAGTACCTGGTCATGGACAGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACAAGCAAGCGGAAGGAGTAcgc                                                                               >  8:56769/1‑149 (MQ=255)
                                                                          ccGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTg                                                                         <  7:59446/149‑1 (MQ=255)
                                                                                   gCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTTCGCGGGATGGTAACCAtt                                                                <  8:71064/149‑1 (MQ=255)
                                                                                   gCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCAt                                                                 <  5:368608/148‑1 (MQ=255)
                                                                                   gCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGGCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCAtt                                                                <  8:116900/149‑1 (MQ=255)
                                                                                     gACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTa                                                              <  4:351166/149‑1 (MQ=255)
                                                                                      aCATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTaa                                                             <  5:30872/149‑1 (MQ=255)
                                                                                       cATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAAt                                                            <  5:192525/149‑1 (MQ=255)
                                                                                         tGGCGATTTTGCTGGTGGAGCAGATCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCt                                                          <  2:355357/149‑1 (MQ=255)
                                                                                           gCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAAGCATTTAATCTGc                                                        >  8:275294/1‑149 (MQ=255)
                                                                                             gATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGc                                                        >  8:334813/1‑147 (MQ=255)
                                                                                             gATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGc                                                        >  3:52392/1‑147 (MQ=255)
                                                                                             gATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGa                                                      <  7:7837/149‑1 (MQ=255)
                                                                                               ttttGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGc                                                    <  7:334813/149‑1 (MQ=255)
                                                                                                   gCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAAc                                                 >  7:280110/1‑148 (MQ=255)
                                                                                                   gCTGGGGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCTAGCCGAAGTTGTACGCGGGCTGGTTACCATTTAATCTGCGAGCAAc                                                 <  2:70876/148‑1 (MQ=255)
                                                                                                          gAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑Gc                                        <  6:245790/149‑1 (MQ=255)
                                                                                                               gTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGAc                                   <  8:280110/149‑1 (MQ=255)
                                                                                                               gTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGAc                                   <  1:2501/149‑1 (MQ=255)
                                                                                                                     cGACTTTGCCGAGGCGCTGGAGGACCAGTACCTGTGCATGGCCAGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACCTGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTACTCTGCGAGCAACGACCTT‑GCGAGACGCTGTc                             <  8:156168/149‑1 (MQ=255)
                                                                                                                      gACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCa                            >  8:83574/1‑149 (MQ=255)
                                                                                                                      gACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCa                            >  5:134206/1‑149 (MQ=255)
                                                                                                                       aCTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCAt                           <  3:126812/149‑1 (MQ=255)
                                                                                                                         tttGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCAt                           >  6:54347/1‑147 (MQ=255)
                                                                                                                           tGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAg                        <  2:344426/148‑1 (MQ=255)
                                                                                                                                    gggccatgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAAtt              >  2:202288/7‑149 (MQ=255)
                                                                                                                                       ccatgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGaa                 >  6:372758/4‑143 (MQ=255)
                                                                                                                                       ccatgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGaa                 <  5:372758/140‑1 (MQ=255)
                                                                                                                                        gctgACCAGTACCTGGTCATGGCCCGGGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTcgcg          <  2:112932/149‑1 (MQ=255)
                                                                                                                                         atgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCt                                <  4:135669/125‑1 (MQ=255)
                                                                                                                                         atgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCt                                >  3:135669/2‑126 (MQ=255)
                                                                                                                                          tgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGAt        >  2:274793/1‑149 (MQ=255)
                                                                                                                                          tgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGAt        >  3:64190/1‑149 (MQ=255)
                                                                                                                                           gACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATg       >  6:11103/1‑149 (MQ=255)
                                                                                                                                             ccAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGcc     <  7:83574/149‑1 (MQ=255)
                                                                                                                                             ccAGTACCTGGTCATGGCCCGTGGCGAAATAATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGcc     >  8:23440/1‑149 (MQ=255)
                                                                                                                                              cAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGccc    >  6:239460/1‑149 (MQ=255)
                                                                                                                                              cAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGGCAGGAGCTACGAAATTCGCGATGccc    >  4:103977/1‑149 (MQ=255)
                                                                                                                                                gTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGTTCGCTGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTTCGCGAGACGCTGTCATGAGCTACGAAATTcg            <  1:34411/140‑1 (MQ=255)
                                                                                                                                                gTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGCCCTg  >  4:12873/1‑149 (MQ=255)
                                                                                                                                                gTACCTGGTCATGGCCCGTGGCGAAATAATCCAGCGCGGTCGCGTCGAAAACATGCAAGCAGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTcg            >  2:34411/1‑139 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGCCCTG  >  NC_002947/5513840‑5514132

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑