Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,513,986 A→T 7.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 7.9% 180.5 / 5.5 63Y204F (TAC→TTC) ‡urtEABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base T (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 = 5.85e-01

GATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGCCCTGCA  >  NC_002947/5513840‑5514134
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                    >  3:351166/1‑148 (MQ=255)
                  cgaAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGc                                                                                                                                    >  6:30872/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                   <  5:372758/140‑1 (MQ=255)
                                                                                                                                       ccatgACCAGGTACTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGaa                   >  6:372758/4‑143 (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          >  3:64190/1‑149 (MQ=255)
                                                                                                                                          tgACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGAt          >  2:274793/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)
                                                                                                                                                  aCCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGCCCTGca  >  2:59603/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATTCTCGACGAGCCGACCGAAGGCATCCAGCCCTCGGTGATCAAGGAAATCGGCGCGGTCATCCGCCGCCTGGCCGAGCGTGGCGACATGGCGATTTTGCTGGTGGAGCAGTTCTACGACTTTGCCGAGGCGCTGGCTGACCAGTACCTGGTCATGGCCCGTGGCGAAATCATCCAGCGCGGTCGCGGCGAAAACATGCAAGCCGAAGGTGTACGCGGGCTGGTAACCATTTAATCTGCGAGCAACGACCTT‑GCGAGACGCTGTCATGAGCTACGAAATTCGCGATGCCCTGCA  >  NC_002947/5513840‑5514134

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