Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,976,207 A→C 6.4% G188G (GGT→GGG)  flgJ ← peptidoglycan hydrolase FlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,976,2070AC 6.4% 171.2 / 4.3 63G188G (GGT→GGG) flgJpeptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base A (25/34);  new base C (2/2);  total (27/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 = 1.10e-01

GCCTGCTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATG  >  NC_002947/4976061‑4976354
                                                                                                                                                  |                                                                                                                                                   
gcCTGCTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGAcc                                                                                                                                                   <  6:242073/149‑1 (MQ=255)
gcCTGCTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGAcc                                                                                                                                                   <  2:116173/149‑1 (MQ=255)
     cTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCAc                                                                                                                                              <  7:305777/149‑1 (MQ=255)
       cGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACgg                                                                                                                                            <  7:250923/149‑1 (MQ=255)
            aTTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcagg                                                                                                                                       <  8:222681/149‑1 (MQ=255)
            aTTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcagg                                                                                                                                       <  1:338338/149‑1 (MQ=255)
             ttGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGggcagg                                                                                                                                       >  1:201263/1‑148 (MQ=255)
                   gcaTGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCAcggcg                                                                                                                                <  7:292894/149‑1 (MQ=255)
                   cattggtgtCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCAcggcg                                                                                                                                <  4:242153/146‑1 (MQ=255)
                      tggtggCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGgcggcg                                                                                                                             >  3:67081/1‑149 (MQ=255)
                       ggtggCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGt                                                                                                                            >  6:110265/1‑149 (MQ=255)
                         tggCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTa                                                                                                                          >  8:285867/1‑149 (MQ=255)
                            cgacgaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTg                                                                                                                       <  5:273551/149‑1 (MQ=255)
                             gacgaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTgg                                                                                                                      >  5:338458/1‑149 (MQ=255)
                             gacgaAGGCGGCGCTGGCGGGGAAAGCCTTGTTCGGCGCCAGCGGGGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTgg                                                                                                                      <  5:381196/149‑1 (MQ=255)
                              acgaAGGCGGGGCGGTCGGTCAACGCCTTGTTCGGCGCCTGCGGTGGCGGGGCCACGGCACGGCCGACAATGCGCCGGCCGTTGCCCGGCACGGCTAAGGCCTTGGCTACTGCCTCACCGTCCCGGGCAGGCACGGCGGCGTTATTGGc                                                                                                                     <  2:201263/149‑1 (MQ=255)
                                gaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCAc                                                                                                                                              <  6:275291/122‑1 (MQ=255)
                                gaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCAc                                                                                                                                              >  5:275291/1‑122 (MQ=255)
                                gaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCGGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAg                                                                                                                   <  7:334917/149‑1 (MQ=255)
                                   ggCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGttt                                                                                                                <  8:111510/149‑1 (MQ=255)
                                    gCGTCGCTGTCGGTGAAAGCCTTGTGCGTCTCCAGAGGTGGCTGGGCCAATGCACGGCGTACAATGCGCAGGCCGTTGTCCGGCACGGCTAAGGCGTTTTCTACTGCCTGACCGTCAGGGGAAGGCACGGCGGCGTTATTGGCAGttt                                                                                                                <  6:281155/148‑1 (MQ=255)
                                       tCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTccgc                                                                                                            >  2:71064/1‑149 (MQ=255)
                                          cTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGTAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGa                                                                                                         >  7:274397/1‑149 (MQ=255)
                                               ggTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCa                                                                                                    >  8:52650/1‑149 (MQ=255)
                                                 tGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCAc                                                                                                   >  1:228940/1‑148 (MQ=255)
                                                    aaGCCTTGTTCTGCGCCAGGGGGGGCTGGGCCACGGCACGGCCTACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATa                                                                                               <  8:187922/149‑1 (MQ=255)
                                                    aaGCCTTGTGCGGCCCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGCCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGAt                                                                                                <  7:52180/148‑1 (MQ=255)
                                                       ccTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGCCCGTGACGgtcagg                                                                                                                                       >  4:255044/1‑106 (MQ=255)
                                                           gTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGc                                                                                          <  2:270910/147‑1 (MQ=255)
                                                           gTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCa                                                                                        >  4:295397/1‑149 (MQ=255)
                                                              cGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCaa                                                                                       >  1:215279/1‑147 (MQ=255)
                                                                gcgcCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAgg                                                                                   >  3:334871/1‑149 (MQ=255)
                                                                 cgcCCGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGACGGCCTTGGCTACTCCCTGACCGTCCCGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGc                                                                                  <  4:60015/149‑1 (MQ=255)
                                                                  gcCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGAAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcg                                                                                 >  5:387698/1‑149 (MQ=255)
                                                                    cAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGTAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggt                                                                               <  8:274397/149‑1 (MQ=255)
                                                                     aGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtc                                                                              >  3:329297/1‑149 (MQ=255)
                                                                      gCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGcggtcg                                                                             <  7:285867/149‑1 (MQ=255)
                                                                          gggCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGa                                                                         <  5:147804/148‑1 (MQ=255)
                                                                              tGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGttt                                                                     >  5:295966/1‑149 (MQ=255)
                                                                                ggCCACGTCACGGCCGACAATGCTCAGTCCTTTGTCCGTCACGGCGAAGGCCTTGTCTACTGCCTGACCGTCACGGGCATGCACGGCGTCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGc                                                                   <  6:387698/149‑1 (MQ=255)
                                                                                  ccACGGCACGGCGGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCtt                                                                 <  8:197882/149‑1 (MQ=255)
                                                                                    acggcacggcCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTgg                                                               >  1:244588/1‑149 (MQ=255)
                                                                                           ggcCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGc                                                         >  5:344861/1‑148 (MQ=255)
                                                                                                 cAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGgcgc                                                  >  2:218922/1‑149 (MQ=255)
                                                                                                    tGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGcgccgc                                               <  6:338458/149‑1 (MQ=255)
                                                                                                       gcAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCcgcgcg                                            >  8:103581/1‑149 (MQ=255)
                                                                                                        cAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGt                                           <  6:295966/149‑1 (MQ=255)
                                                                                                        cAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGt                                           <  5:106263/149‑1 (MQ=255)
                                                                                                         aGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGtt                                          <  7:198441/149‑1 (MQ=255)
                                                                                                           gCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGa                                        <  5:110265/149‑1 (MQ=255)
                                                                                                              gTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTg                                     <  4:334871/149‑1 (MQ=255)
                                                                                                                 ggccGGCACGGCGACGGCCTTGGCTACTGCCGGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTgca                                   <  2:152277/146‑1 (MQ=255)
                                                                                                                     ggCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAAc                              >  4:175214/1‑149 (MQ=255)
                                                                                                                        aCGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAAcgtc                           <  4:67081/149‑1 (MQ=255)
                                                                                                                             gAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTgc                      <  1:71064/149‑1 (MQ=255)
                                                                                                                                   cTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAg                <  7:52650/149‑1 (MQ=255)
                                                                                                                                      ggCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCg              <  4:329297/148‑1 (MQ=255)
                                                                                                                                         tACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCCCGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGAc           >  2:262041/1‑148 (MQ=255)
                                                                                                                                             gCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCTGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGc      >  5:104170/1‑149 (MQ=255)
                                                                                                                                              ccTGCCCGTGACGGTCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCACCAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGc      >  5:393721/1‑148 (MQ=255)
                                                                                                                                                tGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCAt   >  4:144835/1‑149 (MQ=255)
                                                                                                                                                tGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCAt   >  4:43500/1‑149 (MQ=255)
                                                                                                                                                 gACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATg  <  3:295397/149‑1 (MQ=255)
                                                                                                                                                  aCCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATg  <  1:218922/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCTGCTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATG  >  NC_002947/4976061‑4976354

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