Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,224,601 A→C 5.4% M26R (ATG→AGG)  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,6010AC 5.4% 238.9 / 4.0 74M26R (ATG→AGG) gltKglutamate/aspartate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base C (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

GAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCAT  >  NC_002947/1224456‑1224746
                                                                                                                                                 |                                                                                                                                                 
gAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGCCCAt                                                                                                                                                  >  4:226057/1‑147 (MQ=255)
 aaCCAGTTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCa                                                                                                                                                <  7:163290/148‑1 (MQ=255)
   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                                                                                                                <  1:66196/129‑1 (MQ=255)
                                                    gcgcCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAgg                                                                                                                >  2:66196/1‑129 (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                                            >  2:309033/1‑139 (MQ=255)
                                                                                                              tGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcc                                            <  1:309033/139‑1 (MQ=255)
                                                                                                                    gCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttc                            >  5:345705/1‑149 (MQ=255)
                                                                                                                       cGATACCGACGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGTGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttcttc                         <  7:246099/149‑1 (MQ=255)
                                                                                                                        gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttcttc                         >  5:249420/1‑148 (MQ=255)
                                                                                                                        gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                        >  1:361379/1‑149 (MQ=255)
                                                                                                                        gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                        >  8:53278/1‑149 (MQ=255)
                                                                                                                        gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                        <  5:355766/149‑1 (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                       <  3:5964/133‑1 (MQ=255)
                                                                                                                                         ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTtctc                       >  4:5964/1‑133 (MQ=255)
                                                                                                                                         ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGagc       <  4:63649/149‑1 (MQ=255)
                                                                                                                                         ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGag        >  2:182562/1‑148 (MQ=255)
                                                                                                                                              aCCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCAt  <  5:254439/149‑1 (MQ=255)
                                                                                                                                               ccATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATc                                                 >  1:12620/1‑101 (MQ=255)
                                                                                                                                               ccATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATc                                                 <  2:12620/101‑1 (MQ=255)
                                                                                                                                                cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcca                                           <  4:1296/106‑1 (MQ=255)
                                                                                                                                                cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcca                                           >  3:1296/1‑106 (MQ=255)
                                                                                                                                                cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGa         <  5:34433/140‑1 (MQ=255)
                                                                                                                                                cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGa         >  6:34433/1‑140 (MQ=255)
                                                                                                                                                cATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCAt  >  4:76872/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCAT  >  NC_002947/1224456‑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: