Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,975,249 T→A 6.9% T119S (ACT→TCT)  flgK ← flagellar hook‑associated protein FlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,975,2490TA 6.9% 193.3 / 3.0 67T119S (ACT→TCT) flgKflagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base A (3/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGTTGTCCAG  >  NC_002947/4975101‑4975393
                                                                                                                                                      |                                                                                                                                                 
tGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGt                                                                                                                                                   >  2:248063/1‑149 (MQ=255)
  gccggccATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGTGTCAGCTGGGTGCTCACCTCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTc                                                                                                                                                  <  4:17033/148‑1 (MQ=255)
     ggccATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGATAGTGAACGGTTGCTCAGGGTTTGCCCATGGTACAGCAGCACATGGCGCGAGCCAACGACACCCGCATACGCCGCCGCAGACTg                                                                                                                                                >  1:126086/1‑149 (MQ=255)
             tgtCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCgag                                                                                                                                      <  2:53051/149‑1 (MQ=255)
               tCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAATCGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGagaa                                                                                                                                    <  1:288985/149‑1 (MQ=255)
                 gAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGagaaga                                                                                                                                  <  5:117231/149‑1 (MQ=255)
                 gAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGagaag                                                                                                                                   >  4:100355/1‑148 (MQ=255)
                  aGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAgaagaa                                                                                                                                 <  3:100355/149‑1 (MQ=255)
                       ggAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCg                                                                                                                             <  1:248063/148‑1 (MQ=255)
                         aGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGt                                                                                                                           >  6:310842/1‑148 (MQ=255)
                          gTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTaa                                                                                                                         <  6:22323/149‑1 (MQ=255)
                              aTCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGAc                                                                                                                    <  8:247567/149‑1 (MQ=255)
                                 gTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGAcc                                                                                                                   >  1:354040/1‑147 (MQ=255)
                                 gTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACACGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGAcc                                                                                                                   >  6:195925/1‑147 (MQ=255)
                                  tGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCgg                                                                                                                            >  6:76335/1‑138 (MQ=255)
                                  tGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCgg                                                                                                                            <  5:76335/138‑1 (MQ=255)
                                     cGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTg                                                                                                             <  7:357058/149‑1 (MQ=255)
                                           gctgGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCg                                                                                                       >  2:70538/1‑149 (MQ=255)
                                               ggTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGgtgt                                                                                                   >  1:189953/1‑149 (MQ=255)
                                               ggTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGgt                                                                                                     <  3:339748/147‑1 (MQ=255)
                                                        ggTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCCGAGCTGATGCCGGTGTCACGGTCg                                                                                           <  1:3897/149‑1 (MQ=255)
                                                         gTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGc                                                                                         >  8:376006/1‑149 (MQ=255)
                                                         gTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGc                                                                                         <  5:301525/149‑1 (MQ=255)
                                                           gCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCTCCGCCGCAGTCTGCAg                                                                                                                                             <  8:29914/96‑1 (MQ=255)
                                                           gCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAg                                                                                                                                             >  7:29914/1‑96 (MQ=255)
                                                            cTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGTcag                                                                                      <  2:16859/149‑1 (MQ=255)
                                                                ccGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCTCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGt                                                                                  <  2:189953/149‑1 (MQ=255)
                                                                 cgggTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCTCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGtt                                                                                 <  5:314155/146‑1 (MQ=255)
                                                                 cgcgTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGtt                                                                                 >  4:285491/1‑149 (MQ=255)
                                                                    gTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGt                                                                              <  1:70538/149‑1 (MQ=255)
                                                                          aaaCGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGAc                                                                          >  2:175034/1‑149 (MQ=255)
                                                                           aaCGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACg                                                                         >  8:260410/1‑149 (MQ=255)
                                                                                    tCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGtgatc                                                                >  7:298681/1‑149 (MQ=255)
                                                                                      aGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGtgatct                                                               >  6:96573/1‑148 (MQ=255)
                                                                                          tttGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATcc                                                          >  8:44419/1‑149 (MQ=255)
                                                                                             gcgcCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTg                                                        <  5:310842/148‑1 (MQ=255)
                                                                                              cgcCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTgg                                                       >  1:355402/1‑148 (MQ=255)
                                                                                               gcCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGaa                                                     <  2:355402/149‑1 (MQ=255)
                                                                                               gcCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGaa                                                     <  4:186847/149‑1 (MQ=255)
                                                                                                 cTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTg                                                                      >  5:11009/1‑130 (MQ=255)
                                                                                                 cTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTg                                                                      <  6:11009/130‑1 (MQ=255)
                                                                                                 cTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACg                                                   >  5:116978/1‑149 (MQ=255)
                                                                                                    ggTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCa                                                <  8:298681/149‑1 (MQ=255)
                                                                                                      tCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCAcc                                              >  6:311256/1‑149 (MQ=255)
                                                                                                         gcagcagcTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGagaaga                                                                                                                                  >  7:75482/1‑63 (MQ=255)
                                                                                                         gcagcagcTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGagaaga                                                                                                                                  <  8:75482/63‑1 (MQ=255)
                                                                                                         gcagcagcTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAg                                           >  7:100994/1‑149 (MQ=255)
                                                                                                          cagcagcTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGt                                          <  7:376006/149‑1 (MQ=255)
                                                                                                                   gcgcGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGtgatc                                                                >  7:116656/1‑118 (MQ=255)
                                                                                                                   gcgcGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGtgatc                                                                <  8:116656/118‑1 (MQ=255)
                                                                                                                     gcgAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAAc                               <  6:116978/149‑1 (MQ=255)
                                                                                                                       gAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAAc                               >  4:160289/1‑147 (MQ=255)
                                                                                                                          gCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGt                          <  1:175034/149‑1 (MQ=255)
                                                                                                                            aaCGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTgg                        >  4:367611/1‑149 (MQ=255)
                                                                                                                              cGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGcc                      <  5:311256/149‑1 (MQ=255)
                                                                                                                               gTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCg                     >  4:122992/1‑149 (MQ=255)
                                                                                                                                     cGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCa               <  5:96573/149‑1 (MQ=255)
                                                                                                                                        cTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCt            <  2:354040/149‑1 (MQ=255)
                                                                                                                                         ttCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTg           >  5:337885/1‑149 (MQ=255)
                                                                                                                                            ggcgctgcAGACTGCGGCGGCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTgg                        <  2:395867/131‑1 (MQ=255)
                                                                                                                                            ggcgctgcAGACTGCGGCGGCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTgg                        >  1:395867/3‑133 (MQ=255)
                                                                                                                                            gccgccgcATTCTGCAGGGCCGAGAAGAAGGCGGTAAGGACCGATCTGATTCCGTTTTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGTAACGTCAACGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGtt        <  7:255601/149‑1 (MQ=255)
                                                                                                                                            gccgccgcAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGtt        >  5:132244/1‑149 (MQ=255)
                                                                                                                                             gcgctgcAGACTGCGGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGTTg       >  6:332469/2‑149 (MQ=255)
                                                                                                                                             ccgccgcAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGTTg       >  3:219495/1‑149 (MQ=255)
                                                                                                                                              cgccgcAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGTTGt      <  3:367611/149‑1 (MQ=255)
                                                                                                                                                  gcATTCTGCAGCGCCGAGTAGAAGGCGTTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGTTGTCCAg  <  5:195925/149‑1 (MQ=255)
                                                                                                                                                   cAGACTGCGGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGTTGTCCa   <  5:332469/147‑1 (MQ=255)
                                                                                                                                                      tCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTctgcct                                    >  4:273646/1‑111 (MQ=255)
                                                                                                                                                      tCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTctgcct                                    <  3:273646/111‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
TGGCCGGCCATGGTGTCGAGCTGGGAGTTGATCGTGGCGTTCTGCTGGGTCAGCTGGGTGCTCACCGCGTTGAAACGGTTGCTCAGGGTTTGCGCCTGGGTCAGCAGCAGCTGGCGCGAGGCAACGTCACCCGGCTTCGCCGCCGCAGTCTGCAGCGCCGAGAAGAAGGCGGTAAGGACCGAGCTGATGCCGGTGTCACGGTCGGCCAGCAGTTTGTCGACGCTGGTGATCTGATCCTGGAACGTCACCGAGTCTGCCTGCAACGAGGTGGCCGTCTGCAGCTGGTTGTCCAG  >  NC_002947/4975101‑4975393

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: