Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,224,464 T→A 6.1% T72S (ACC→TCC)  gltK ← glutamate/aspartate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,224,4640TA 6.1% 198.8 / 2.8 65T72S (ACC→TCC) gltKglutamate/aspartate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (28/33);  new base A (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GACTGCACACCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAG  >  NC_002947/1224319‑1224611
                                                                                                                                                 |                                                                                                                                                     
gACTGCACACCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGAt                                                                                                                                                    >  5:166597/1‑149 (MQ=255)
 aCTGCACACCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATc                                                                                                                                                   <  8:235428/149‑1 (MQ=255)
   tGCACACCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGAGATCAc                                                                                                                                                 >  7:300859/1‑149 (MQ=255)
   tGCACACCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGAGATCAc                                                                                                                                                 >  3:232642/1‑149 (MQ=255)
     cacaCCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCAc                                                                                                                                                 >  3:266541/1‑147 (MQ=255)
     cacaCCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGAGATCACca                                                                                                                                               <  8:300859/149‑1 (MQ=255)
      acacCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGTCCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGAGATCACca                                                                                                                                               <  4:232642/148‑1 (MQ=255)
        acCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACcag                                                                                                                                              >  4:108252/1‑147 (MQ=255)
             cgctGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGTTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGAGCACCagcagcagc                                                                                                                                       <  8:210709/149‑1 (MQ=255)
              gcgGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCg                                                                                                                                      <  5:310181/149‑1 (MQ=255)
                   cGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATc                                                                                                                                 <  6:260491/149‑1 (MQ=255)
                    gATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGACGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATc                                                                                                                                 >  1:130629/1‑148 (MQ=255)
                     aTTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGa                                                                                                                               <  1:173920/149‑1 (MQ=255)
                      tttCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAg                                                                                                                              >  7:100772/1‑149 (MQ=255)
                          gCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCg                                                                                                                                      >  1:312246/1‑137 (MQ=255)
                          gCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCg                                                                                                                                      <  2:312246/137‑1 (MQ=255)
                          gCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCgg                                                                                                                           >  3:45160/1‑148 (MQ=255)
                               aaTACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGCATCGAGCGGAagtag                                                                                                                     >  7:125382/1‑149 (MQ=255)
                               aaTACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAagtag                                                                                                                     <  3:108252/149‑1 (MQ=255)
                                 tACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑tt                                                                                                                 <  8:125382/149‑1 (MQ=255)
                                      ggCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACg                                                                                                             >  2:263404/1‑148 (MQ=255)
                                       gCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTa                                                                                                           <  2:60274/149‑1 (MQ=255)
                                            gAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGgcgc                                                                                                      <  4:266541/149‑1 (MQ=255)
                                            gAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGgc                                                                                                        <  8:15561/147‑1 (MQ=255)
                                              aCATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCg                                                                                                    <  6:166597/149‑1 (MQ=255)
                                              aCATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCg                                                                                                    >  2:23068/1‑149 (MQ=255)
                                                  catGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGc                                                                                                  >  3:288386/1‑147 (MQ=255)
                                                   atGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAg                                                                                               >  4:133768/1‑149 (MQ=255)
                                                    tGAAGACGACCACGCAGGAGGTGAAGTCACCCACCGGCGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGCTACCGCCAGGTAGAACCAGGTGATCCCCAGCAGCAGCGGAATCGAGCTTAAGTAG‑‑TTGACGTAGGCGCCGGCCAgg                                                                                              >  8:150588/1‑149 (MQ=255)
                                                     gAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGt                                                                                             >  5:100091/1‑149 (MQ=255)
                                                       aGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTc                                                                                           <  5:120083/149‑1 (MQ=255)
                                                            aCCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGAcag                                                                                      <  4:135072/149‑1 (MQ=255)
                                                             ccACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCAc                                                                                                                                                 <  6:244100/91‑1 (MQ=255)
                                                             ccACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCAc                                                                                                                                                 >  5:244100/1‑91 (MQ=255)
                                                              cACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGAcagca                                                                                    <  4:288386/149‑1 (MQ=255)
                                                              cACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGATACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGAcagca                                                                                    <  1:212385/149‑1 (MQ=255)
                                                                   aggaggTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTg                                                                               >  7:83052/1‑149 (MQ=255)
                                                                   aggaggTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTg                                                                               >  2:278801/1‑149 (MQ=255)
                                                                    ggaggTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGc                                                                              <  7:74538/149‑1 (MQ=255)
                                                                         tGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTgac                                                                         >  3:305048/1‑149 (MQ=255)
                                                                          gAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGt                                                                                                                                                       <  4:81311/72‑1 (MQ=255)
                                                                          gAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGt                                                                                                                                                       >  3:81311/1‑72 (MQ=255)
                                                                             ggCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTgacga                                                                       >  4:32969/1‑147 (MQ=255)
                                                                              gCACCCACCGGGGTATCTTCGCCGGTGATCCAGAGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGACGCCCAGGTTCGACAGCAGTTTGCTTGACGACAg                                                                    >  5:74014/1‑149 (MQ=255)
                                                                                acccaccGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGgc                                                                  <  3:133768/149‑1 (MQ=255)
                                                                                        gggTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGt                                                          <  1:23068/149‑1 (MQ=255)
                                                                                          gTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGc                                                        <  6:100091/149‑1 (MQ=255)
                                                                                               ttCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTCTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAg                                                     <  1:263404/149‑1 (MQ=255)
                                                                                                  gCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATg                                                 <  1:278801/148‑1 (MQ=255)
                                                                                                   ccGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGt                                               >  2:173665/1‑149 (MQ=255)
                                                                                                          tCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAg                                         >  6:184067/1‑148 (MQ=255)
                                                                                                           ccAGCGCAGCACGAACGGTACCGCCAGGTAGAACTAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGAGGTAGGCGCCGGCCAGGTTCGACAGCAGTTTg                                                                               >  1:105498/1‑109 (MQ=255)
                                                                                                           ccAGCGCAGCACGAACGGTACCGCCAGGTAGAACTAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTg                                                                               <  2:105498/109‑1 (MQ=255)
                                                                                                            cAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCa                                       <  8:83052/148‑1 (MQ=255)
                                                                                                                       gAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGGAGTTTGCTTGACGACAGGCGCATCAGTGcc                                                       <  1:30162/121‑1 (MQ=255)
                                                                                                                       gAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGcc                                                       >  2:30162/1‑121 (MQ=255)
                                                                                                                       gAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGa                           <  3:97981/149‑1 (MQ=255)
                                                                                                                            gTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCTGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACt                      <  6:132009/149‑1 (MQ=255)
                                                                                                                            gTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACt                      >  1:45784/1‑149 (MQ=255)
                                                                                                                             tACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGTTACCGACGATGACTc                     >  2:185735/1‑149 (MQ=255)
                                                                                                                                  ccAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACg                                     <  7:242504/128‑1 (MQ=255)
                                                                                                                                  ccAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACg                                     >  8:242504/1‑128 (MQ=255)
                                                                                                                                    aGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGAc              <  2:45784/149‑1 (MQ=255)
                                                                                                                                         gAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCTCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCa         <  4:45160/149‑1 (MQ=255)
                                                                                                                                             cAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGttt     >  6:11085/1‑149 (MQ=255)
                                                                                                                                               ggTGATCACCAGCAGCAGCGGATTCGATCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCATGATGGTGCGCATCACGATACCTCCGATGACTCCCATGTCCATCAGTTTGa   <  6:74014/149‑1 (MQ=255)
                                                                                                                                                gTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAg  <  1:173665/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GACTGCACACCGGCGCGGACGATTTCGCAGAAATACGCGGCCTCGAACATCATGAAGGCGACCACGCAGGAGGTGAAGGCACCCACCGGGGTATCTTCGCCGGTGATCCAGCGCAGCACGAACGGTACCGCCAGGTAGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAG‑‑TTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAG  >  NC_002947/1224319‑1224611

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: