Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,224,607 T→C 5.6% K24R (AAA→AGA) ‡ 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,6070TC 5.6% 211.5 / 2.2 71K24R (AAA→AGA) ‡gltKglutamate/aspartate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base C (2/2);  total (36/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 = 7.53e-01

CCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCAT  >  NC_002947/1224459‑1224746
                                                                                                                                                    |                                                                                                                                           
ccAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGtt                                                                                                                                             >  7:92309/1‑149 (MQ=255)
   ggTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGa                                                                                                                                          <  1:6537/149‑1 (MQ=255)
   cacattCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGa                                                                                                                                          <  6:346550/144‑1 (MQ=255)
    gTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAg                                                                                                                                         >  7:60872/1‑149 (MQ=255)
        tCACCAGCAGCAGCGGAATCGAGCGGAAGTATTTGACGTAGGCGCCGGCCAGGTTCGACAGCATTTTGGTTGACGACAGGCGCATCAGTGCCAGGATGTTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTTAGGGTc                                                                                                                                     <  8:343167/149‑1 (MQ=255)
          aCCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCAt                                                                                                                                   <  6:202295/149‑1 (MQ=255)
           ccAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCAt                                                                                                                                   >  5:124545/1‑148 (MQ=255)
             agcagcagcGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCCTCAAACTGATGGTCATggg                                                                                                                                <  5:83597/149‑3 (MQ=255)
                       gAATCGAGCGGAAGTACTTGACGTAGGCGCCGGCCAGGTTCGACAGCACGTTGCTTGACGACAGGCGCATCAGTGCCAGGAGGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGGCATGACCATGCCCTcc                                                                                                                      <  8:92309/149‑1 (MQ=255)
                          tCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATTCCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCag                                                                                                                   <  2:118343/149‑1 (MQ=255)
                            gAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGAAAGGCGCATCAGTGCAAGGATGGTGCCCAGAACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCagag                                                                                                                 >  8:62169/1‑149 (MQ=255)
                             aGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAgg                                                                                                                >  6:189946/1‑149 (MQ=255)
                                 gAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcag                                                                                                            <  5:189946/149‑1 (MQ=255)
                                 gAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcag                                                                                                            >  4:352823/1‑149 (MQ=255)
                                          gACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCgg                                                                                                   >  5:135444/1‑149 (MQ=255)
                                          gACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCg                                                                                                    <  3:226057/148‑1 (MQ=255)
                                           aCGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCgg                                                                                                   <  8:231172/148‑1 (MQ=255)
                                           aCGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGACGATGACTCCCATGACCATCAGTTTGTGGGTCATGACCATGCCCTCCCAGAGGGAAGGCAGGTCggg                                                                                                  >  8:246099/1‑149 (MQ=255)
                                             gTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCACGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGga                                                                                                >  4:187484/1‑149 (MQ=255)
                                              tAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGgat                                                                                               <  8:60872/149‑1 (MQ=255)
                                               aGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggca                                                                                                             >  3:350338/1‑134 (MQ=255)
                                               aGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggca                                                                                                             <  4:350338/134‑1 (MQ=255)
                                                 gcgcCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAgg                                                                                                                >  2:66196/1‑129 (MQ=255)
                                                 gcgcCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAgg                                                                                                                <  1:66196/129‑1 (MQ=255)
                                                     cGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCACGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCg                                                                                        <  3:187484/149‑1 (MQ=255)
                                                       ggcAGGTTCGACAGGAGTTTGCTTGACGACAGGCTCATCAGTGCCATGATGGTGCCAGCACGATACCGCCGATGTCTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTg                                                                                     <  6:205560/147‑1 (MQ=255)
                                                       gCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAgg                                                                                                                <  6:184763/123‑1 (MQ=255)
                                                       gCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAgg                                                                                                                >  5:184763/1‑123 (MQ=255)
                                                       gCCAGGGTCGACAGCAGGTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCt                                                                                      <  5:320087/149‑1 (MQ=255)
                                                                gACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccat                                                                             >  5:118015/1‑149 (MQ=255)
                                                                gACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCAGCGAAAtcatc                                                                             <  6:193273/149‑4 (MQ=255)
                                                                            cTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACtt                                                                 >  5:303374/1‑149 (MQ=255)
                                                                                       gcgcATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAt                                                                                 <  4:368324/122‑1 (MQ=255)
                                                                                         gcATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGa                                                     >  6:254439/1‑148 (MQ=255)
                                                                                           aTCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAAt                                                  >  7:12406/1‑149 (MQ=255)
                                                                                              aGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATca                                                <  6:303374/148‑1 (MQ=255)
                                                                                               gTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcagg                                              >  6:355766/1‑149 (MQ=255)
                                                                                               gTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcagg                                              <  6:124545/149‑1 (MQ=255)
                                                                                                tGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggc                                             >  2:193070/1‑149 (MQ=255)
                                                                                                 gCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggc                                             >  1:64004/1‑148 (MQ=255)
                                                                                                  ccAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggc                                             >  3:63649/1‑147 (MQ=255)
                                                                                                  ccAGGATGGTGCCCAGCACGATACCGCCGATGACTACCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGGTGAAATCCATTTCCATTTACTTGACCCCCACGGAAATcaggaca                                           >  3:9505/1‑149 (MQ=255)
                                                                                                          gTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACgg                                    <  2:64004/148‑1 (MQ=255)
                                                                                                           tGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcc                                            <  1:309033/139‑1 (MQ=255)
                                                                                                           tGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcc                                            >  2:309033/1‑139 (MQ=255)
                                                                                                                 gCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttc                            >  5:345705/1‑149 (MQ=255)
                                                                                                                    cGATACCGACGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGTGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttcttc                         <  7:246099/149‑1 (MQ=255)
                                                                                                                     gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttcttc                         >  5:249420/1‑148 (MQ=255)
                                                                                                                     gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                        <  5:355766/149‑1 (MQ=255)
                                                                                                                     gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                        >  8:53278/1‑149 (MQ=255)
                                                                                                                     gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                        >  1:361379/1‑149 (MQ=255)
                                                                                                                      aTACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTtctc                       <  6:135444/149‑1 (MQ=255)
                                                                                                                      aTACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTtctc                       <  3:352823/149‑1 (MQ=255)
                                                                                                                          cgccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGAc                    >  6:152338/1‑148 (MQ=255)
                                                                                                                            ccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCAt                 >  5:345711/1‑149 (MQ=255)
                                                                                                                                   tggtcATGACCCTCAAACTGATGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcagg                                              >  2:150223/5‑113 (MQ=38)
                                                                                                                                   tggtcATGACCCTCAAACTGATGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcagg                                              <  1:150223/109‑1 (MQ=38)
                                                                                                                                   cTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATACATTTCAATTTACTTGCACACCACGGAAATAAGGCAAGGCACGTCTCATTTCTTCTCGACCATGCGAAt           >  8:157946/1‑148 (MQ=255)
                                                                                                                                   cTCCCATGAACATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCAt           <  7:62169/148‑1 (MQ=255)
                                                                                                                                      tcATGACCCTCAAACTGATGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGag        >  1:183766/2‑148 (MQ=255)
                                                                                                                                      ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTtctc                       >  4:5964/1‑133 (MQ=255)
                                                                                                                                      ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTtctc                       <  3:5964/133‑1 (MQ=255)
                                                                                                                                      ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGagc       <  4:63649/149‑1 (MQ=255)
                                                                                                                                      ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGag        >  2:182562/1‑148 (MQ=255)
                                                                                                                                           aCCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCAt  <  5:254439/149‑1 (MQ=255)
                                                                                                                                            ccATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATc                                                 <  2:12620/101‑1 (MQ=255)
                                                                                                                                            ccATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATc                                                 >  1:12620/1‑101 (MQ=255)
                                                                                                                                             cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcca                                           <  4:1296/106‑1 (MQ=255)
                                                                                                                                             cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcca                                           >  3:1296/1‑106 (MQ=255)
                                                                                                                                             cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGa         >  6:34433/1‑140 (MQ=255)
                                                                                                                                             cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGa         <  5:34433/140‑1 (MQ=255)
                                                                                                                                             cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCAt  >  4:76872/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCAT  >  NC_002947/1224459‑1224746

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: