Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,478,782 G→C 5.5% I33M (ATC→ATG)  PP_4816 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,478,7820GC 5.5% 235.1 / 3.6 73I33M (ATC→ATG) PP_4816hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base C (2/2);  total (35/38)
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.65e-01

TGTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGACC  >  NC_002947/5478637‑5478926
                                                                                                                                                 |                                                                                                                                                 
tGTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATc                                                                                                                                                <  5:20285/149‑1 (MQ=255)
   ccTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGc                                                                                                                                             >  6:108451/1‑149 (MQ=255)
       ggCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGccc                                                                                                                                           <  6:127950/147‑1 (MQ=255)
        gCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCgg                                                                                                                                         <  7:218213/148‑1 (MQ=255)
        gCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCgg                                                                                                                                         <  7:166271/148‑1 (MQ=255)
         cATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATCATCGGcctgttg                                                                                                                                       <  2:17344/149‑6 (MQ=255)
          atatACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGCTCGGCCCGGGt                                                                                                                                       >  6:121453/1‑148 (MQ=255)
              aCCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTaaaaa                                                                                                                                  >  1:167492/1‑149 (MQ=255)
                cATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACg                                                                                                                                >  6:146757/1‑149 (MQ=255)
                cATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACg                                                                                                                                >  2:102218/1‑149 (MQ=255)
                 aTCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGTAACAACAGGCCGATGATCGGCCCGGGTAAAAACgg                                                                                                                               <  8:197325/149‑1 (MQ=255)
                    aCCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCa                                                                                                                             <  8:257871/148‑1 (MQ=255)
                     ccccGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAc                                                                                                                            >  7:225823/1‑148 (MQ=255)
                      cccGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCa                                                                                                                             <  2:108204/146‑1 (MQ=255)
                      cccGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCa                                                                                                                             >  1:108204/1‑146 (MQ=255)
                      cccGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑g                                                                                                                          >  8:25888/1‑148 (MQ=255)
                       ccGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑Gaa                                                                                                                        >  5:223679/1‑149 (MQ=255)
                          aCTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAg                                                                                                                     <  5:146757/149‑1 (MQ=255)
                           cTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAgg                                                                                                                    <  4:298083/149‑1 (MQ=255)
                            tGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATCATCGGcctg                                                                                                                                          <  5:261025/127‑3 (MQ=255)
                            tGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATCATCGGcctg                                                                                                                                          >  6:261025/1‑125 (MQ=255)
                              ccGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGtga                                                                                                                 <  7:265114/149‑1 (MQ=255)
                              ccGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGtga                                                                                                                 <  5:121453/149‑1 (MQ=255)
                                 gcggcACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGtgattg                                                                                                              <  6:262988/149‑1 (MQ=255)
                                    gcACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATc                                                                                                           <  6:159427/149‑1 (MQ=255)
                                            agcagcGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTccc                                                                                                   <  8:74071/149‑1 (MQ=255)
                                              cagcGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCca                                                                                                 <  5:108451/149‑1 (MQ=255)
                                                gcGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAACACGGCAC‑GAACAGGTGATTGATCGCCGTCCCcag                                                                                                >  4:209095/1‑148 (MQ=255)
                                                gcGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCcagc                                                                                               <  7:25888/149‑1 (MQ=255)
                                                 cGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGc                                                                                                                              <  8:273459/118‑1 (MQ=255)
                                                 cGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGc                                                                                                                              >  7:273459/1‑118 (MQ=255)
                                                    gagaTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCt                                                                                           <  1:102218/149‑1 (MQ=255)
                                                    gagaTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCt                                                                                           >  2:44805/1‑149 (MQ=255)
                                                      gaTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGGAACAGGTGATTGATCGCCGTCCCCAGCAGCTg                                                                                          <  1:44805/149‑1 (MQ=255)
                                                       aTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCca                                                                                                 <  6:57985/140‑1 (MQ=255)
                                                       aTAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCca                                                                                                 >  5:57985/1‑140 (MQ=255)
                                                        tAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGaa                                                                                       <  6:223679/149‑1 (MQ=255)
                                                            gcagcaAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAg                                                                                    >  1:191709/1‑148 (MQ=255)
                                                            gcagcaAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAg                                                                                    >  6:17305/1‑148 (MQ=255)
                                                              agcaAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCCC‑GAACAGGTGATTGATCGCCGTCCCCAGAAGCTGGAACAGTAc                                                                                 >  5:254784/1‑149 (MQ=255)
                                                                caAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTAcca                                                                               <  2:6109/149‑1 (MQ=255)
                                                                    ctgctTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTAccagcca                                                                           <  8:225823/149‑1 (MQ=255)
                                                                     tgctTGCGGACTCGTTCAGCGTATTGCCGACCTCACCTCGGGCCATGAGGAAGAGCAGCAGCAACAACATGGCGATGATCGGCCCGTGTAATAACGGCAC‑GAACAGGTGTTTGTTCGCCGTCGCCAGGAGCTGGAACAGTAccagccag                                                                          <  7:2971/149‑1 (MQ=255)
                                                                       ctTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGt                                                                        <  7:263170/149‑1 (MQ=255)
                                                                       ctTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGt                                                                        <  5:17305/149‑1 (MQ=255)
                                                                               cTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAACAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAACCGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTAGAACAGTACCAGCCAGGTCAcaccac                                                                >  1:234486/1‑149 (MQ=255)
                                                                                tCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTAccagccag                                                                          >  6:330122/1‑138 (MQ=255)
                                                                                  gTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAg                                                                                                                     <  5:158401/93‑1 (MQ=255)
                                                                                  gTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAg                                                                                                                     >  6:158401/1‑93 (MQ=255)
                                                                                  gTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACg                                                               <  2:167492/147‑1 (MQ=255)
                                                                                   ttCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCaa                                                            >  2:134885/1‑149 (MQ=255)
                                                                                     cAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCaa                                                            <  4:53675/147‑1 (MQ=255)
                                                                                     cAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCaa                                                            >  3:53675/1‑147 (MQ=255)
                                                                                       gCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAg                                                         >  7:304317/1‑148 (MQ=255)
                                                                                           cTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCAt                                                      >  3:181913/1‑147 (MQ=255)
                                                                                             tGCCGACCTCCCCGCCGCCCAGGCGGAACAGCAGCAGCACCACCACCCCGATGATCGGCCCGCGTAAAACCGGCAC‑GAACCGGTGACTGAGCGCCGTCCCCCGCAGCTGGAACACTACCAGCCAGGTCAAACCACGCAACACCATAAGc                                                  <  4:44580/149‑1 (MQ=255)
                                                                                                    cTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGaa                                                                                       >  4:238529/1‑105 (MQ=255)
                                                                                                      cACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATaa                                                    <  6:149699/138‑1 (MQ=255)
                                                                                                      cACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATaa                                                    >  5:149699/1‑138 (MQ=255)
                                                                                                      cACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTc                                         >  3:192195/1‑149 (MQ=255)
                                                                                                         cgcgGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCggg                                      <  2:27220/149‑1 (MQ=255)
                                                                                                                     ggAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACa                                                          >  5:206692/1‑117 (MQ=255)
                                                                                                                     ggAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACa                                                          <  6:206692/117‑1 (MQ=255)
                                                                                                                      gAAGAGCAGCAGCAACAACAGGCCGATGAGCGGCCCGCGTAAAAACGGCAA‑GAACAGGTGATTGATCGCCGTCCACAGCAGCTGGAACAGTACCATCCAGGTCAACACACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCAt                          >  3:82/1‑148 (MQ=255)
                                                                                                                          agcagcagcAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATaa                     <  3:50997/149‑1 (MQ=255)
                                                                                                                              gcagcaACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCAc                 >  4:295558/1‑149 (MQ=255)
                                                                                                                                  caacaaCAGGCCGATGATCGGCCCTGGTATAAACGGCAG‑GAACAGGTGTGTGATCGCCGTCCCCAGAAGCTGGAACAGTACCAGCCAGGTCAACCCACGCAACAGCATAAGCATCTCCCTCGGGCATGGACGCCATTATAAGCACGCTc             <  6:254784/149‑1 (MQ=255)
                                                                                                                                    acaacaGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCgg           >  5:219147/1‑149 (MQ=255)
                                                                                                                                      cccgggCCGATCATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCt         >  5:263969/5‑149 (MQ=255)
                                                                                                                                         aGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTa      <  1:123036/149‑1 (MQ=255)
                                                                                                                                           gCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGa    >  4:158850/1‑149 (MQ=255)
                                                                                                                                             cGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGAcc  >  5:151437/1‑149 (MQ=255)
                                                                                                                                             cGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGAcc  >  2:20385/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGTCCTTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCAC‑GAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGGCTGTAGACC  >  NC_002947/5478637‑5478926

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