Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,976,212 C→G 6.4% D187H (GAC→CAC)  flgJ ← peptidoglycan hydrolase FlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,976,2120CG 6.4% 188.5 / 3.7 63D187H (GAC→CAC) flgJpeptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base C (26/33);  new base G (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

CTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCAT  >  NC_002947/4976066‑4976358
                                                                                                                                                  |                                                                                                                                                  
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                                                                                                                                                  >  5:275291/1‑122 (MQ=255)
                           gaAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCAc                                                                                                                                                  <  6:275291/122‑1 (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)
                                                                                                                                                gTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGGAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCAt  >  8:179851/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTCGGCCATTGGCAGCATGGTGGCGACGAAGGCGTCGCTGTCGGTGAAAGCCTTGTTCGGCGCCAGCGGTGGCTGGGCCACGGCACGGCCGACAATGCGCAGGCCGTTGTCCGGCACGGCGAAGGCCTTGGCTACTGCCTGACCGTCACGGGCAGGCACGGCGGCGTTATTGGCAGTTTCCGCCGACGGCACGATACCGGCCAACAGGCGGTCGGTGAGTTTGCTTGGCAGGGCCAGGCGCCGCGCGTTGAGTGCAGCAACGTCGTTGCGCGTAGCGGACGAGGCCATGGCAT  >  NC_002947/4976066‑4976358

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