Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,224,632 A→G 7.4% W16R (TGG→CGG) ‡ 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,6320AG 7.4% 190.2 / 4.8 67W16R (TGG→CGG) ‡gltKglutamate/aspartate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (34/28);  new base G (2/3);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

GAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAG  >  NC_002947/1224487‑1224768
                                                                                                                                                  |                                                                                                                                        
gAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCagag                                                                                                                                       <  4:263050/149‑1 (MQ=255)
  gCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAggg                                                                                                                                     <  7:240096/149‑1 (MQ=255)
     gAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcag                                                                                                                                  >  8:74539/1‑149 (MQ=255)
     gAAGTAGTTGACGTAGGCGCCGGCCAGCTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATCGTGCCCAGCACGATACCGCCGATGACTCCCATGACCAACAGTTTGAGGGTCATGACAATGCCCTCCCAGAGggcag                                                                                                                                  >  6:94134/1‑149 (MQ=255)
      aaGTAGTTGACGTAGGCGCCGGCCAGAGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATTCCCTCCCAGAGggcag                                                                                                                                  >  1:3403/1‑149 (MQ=255)
        gtagtTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggcatggt                                                                                                                               <  5:49892/149‑5 (MQ=255)
           gtTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcag                                                                                                                                  <  1:187720/143‑1 (MQ=255)
           gtTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcag                                                                                                                                  >  2:187720/1‑143 (MQ=255)
                      cgcCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGAtt                                                                                                                 >  1:138674/1‑149 (MQ=255)
                           gCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCgggg                                                                                                                       >  6:90040/1‑138 (MQ=255)
                           gCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCgggg                                                                                                                       <  5:90040/138‑1 (MQ=255)
                            ccAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggcatg                                                                                                                                 >  4:264656/1‑125 (MQ=255)
                            ccAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggcatg                                                                                                                                 <  3:264656/127‑3 (MQ=255)
                                 gTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAt                                                                                                       >  2:100561/1‑149 (MQ=255)
                                     gACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccat                                                                                                   <  7:74539/149‑1 (MQ=255)
                                        agcagTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccatttc                                                                                                <  3:211549/149‑1 (MQ=255)
                                                    gacgacAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGcc                                                                                    >  2:312322/1‑149 (MQ=255)
                                                    gacgacAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACAATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGcc                                                                                    >  2:306770/1‑149 (MQ=255)
                                                     acgacAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGccc                                                                                   <  1:312322/149‑1 (MQ=255)
                                                        acAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGccccc                                                                                 >  5:81575/1‑148 (MQ=255)
                                                          aGGCGCATCAGTGCCATGATGGTGCCCAGCACGATACCGCCGATGACTCACATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCAc                                                                              <  2:70736/149‑1 (MQ=255)
                                                          aGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCAc                                                                              >  4:285577/1‑149 (MQ=255)
                                                                   aGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATca                                                                      <  1:100561/148‑1 (MQ=255)
                                                                    gTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcagg                                                                    >  1:206512/1‑149 (MQ=255)
                                                                     tGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggc                                                                   >  8:276900/1‑149 (MQ=255)
                                                                        cAGGATGGTGCCCAGCACGTTAGCGCGGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCGCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggccag                                                                <  5:94134/149‑1 (MQ=255)
                                                                        cAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtc                                                                                                      >  7:192503/1‑111 (MQ=255)
                                                                        cAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtc                                                                                                      <  8:192503/111‑1 (MQ=255)
                                                                            aTGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCAc                                                            <  6:81575/149‑1 (MQ=255)
                                                                                tGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGc                                                         >  3:189047/1‑148 (MQ=255)
                                                                                tGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATACATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGc                                                         >  5:180811/1‑148 (MQ=255)
                                                                                  cccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGAc                                                      >  3:205485/1‑149 (MQ=255)
                                                                                  cccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGAc                                                      >  6:198516/1‑149 (MQ=255)
                                                                                  cccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGCAAATCAGGCCAGGCACGGCGa                                                       >  4:239911/1‑148 (MQ=255)
                                                                                   ccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                                                     >  4:135422/1‑149 (MQ=255)
                                                                                          gATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                                              <  5:198516/149‑1 (MQ=255)
                                                                                            tACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCg                                            >  2:53903/1‑149 (MQ=255)
                                                                                              ccgccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGa                                           <  2:206512/148‑1 (MQ=255)
                                                                                              ccgccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGAc                                          <  1:53903/149‑1 (MQ=255)
                                                                                               cgccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGCTGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGAc                                          >  5:240177/1‑148 (MQ=255)
                                                                                               cgccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGAcc                                         <  1:306770/149‑1 (MQ=255)
                                                                                               cgccgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGAcc                                         <  4:205485/149‑1 (MQ=255)
                                                                                                  cgATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATg                                      <  3:239911/149‑1 (MQ=255)
                                                                                                   gATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttc                                                  <  4:288956/136‑1 (MQ=255)
                                                                                                   gATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttc                                                  >  3:288956/1‑136 (MQ=255)
                                                                                                    aTGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATgcg                                    <  3:135422/149‑1 (MQ=255)
                                                                                                       aCTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCAt                                 <  7:276900/149‑1 (MQ=255)
                                                                                                        cTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCAt                                 >  3:198388/1‑148 (MQ=255)
                                                                                                         tCCCATGCCCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGACCCCCACGGAAATCAGGCCAGGCACCGCGCCTTTCTTATCGACCATGCGCAt                                 >  7:235517/1‑147 (MQ=255)
                                                                                                           tcATGACCCTCAAACTGATGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGag                              >  2:76201/2‑148 (MQ=255)
                                                                                                           ccATGACCCTGAGGTTGAGTGTCATGACCATGCCCTCCCAGAGGGGAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTTCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCTTGCGCATGagc                             <  8:235517/149‑1 (MQ=255)
                                                                                                               gACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAgcagca                         >  7:83616/1‑149 (MQ=255)
                                                                                                                 ccATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGCCCATGCGCATGAGCAGCATc                       <  8:83616/149‑1 (MQ=255)
                                                                                                                       gTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGccc                 <  4:189047/149‑1 (MQ=255)
                                                                                                                        tttGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                                                     <  6:63637/112‑1 (MQ=255)
                                                                                                                        tttGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                                                     >  5:63637/1‑112 (MQ=255)
                                                                                                                          tGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCAt               <  4:198388/148‑1 (MQ=255)
                                                                                                                          tGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATg              >  2:9462/1‑149 (MQ=255)
                                                                                                                             gggTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGtt            >  7:240134/1‑148 (MQ=255)
                                                                                                                                       cctGCCCTCTGGGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                                                     >  6:200547/3‑97 (MQ=255)
                                                                                                                                       cctGCCCTCTGGGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                                                     <  5:200547/95‑1 (MQ=255)
                                                                                                                                         tGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTAATTGCCCCCCACGGAAATCAGGCCAGGAACGGCGACTTTCTTCTCGACAATGAGCATGAGAAGCATCAGGCCCATGTTCAGGGTGAAg  >  5:246296/1‑146 (MQ=255)
                                                                                                                                          gCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                                              >  1:71795/1‑101 (MQ=255)
                                                                                                                                          gCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTtcttct                                              <  2:71795/101‑1 (MQ=255)
                                                                                                                                          gCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCa                      <  5:172693/125‑1 (MQ=255)
                                                                                                                                          gCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCa                      >  6:172693/1‑125 (MQ=255)
                                                                                                                                                ccAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGc                                                         >  7:254436/1‑84 (MQ=255)
                                                                                                                                                ccAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGc                                                         <  8:254436/84‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
GAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAGCATCAGGCCCATGTTCAGGGTGAAG  >  NC_002947/1224487‑1224768

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: