Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,950,621 C→G 7.1% P318P (CCG→CCC flhB ← flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,950,6210CG 7.1% 973.6 / 46.2 296P318P (CCG→CCCflhBflagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (130/143);  new base G (13/8);  total (145/151)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

TCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGGTTGGT  >  NC_002947/4950472‑4950770
                                                                                                                                                     |                                                                                                                                                     
   ttCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGc                                                                                                                                                    >  1:1117359/1‑150 (MQ=255)
   ttCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGc                                                                                                                                                    >  1:4932610/1‑150 (MQ=255)
       gTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGGTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGCtt                                                                                                                                                >  2:1111357/1‑150 (MQ=255)
        tGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTc                                                                                                                                               <  1:4856464/150‑1 (MQ=255)
         ggCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGTTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtcta                                                                                                                                              <  1:5662785/150‑4 (MQ=255)
         ggCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtcta                                                                                                                                              <  1:3565361/150‑4 (MQ=255)
          gCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGa                                                                                                                                             <  2:4028225/150‑1 (MQ=255)
          gCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGa                                                                                                                                             <  2:40505/150‑1 (MQ=255)
          gCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGa                                                                                                                                             <  2:717764/150‑1 (MQ=255)
           ctctGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAg                                                                                                                                            >  2:1405558/1‑150 (MQ=255)
            tctGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtctatta                                                                                                                                           <  1:2156848/150‑7 (MQ=255)
            tctGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAgg                                                                                                                                           >  2:2859541/1‑150 (MQ=255)
                  cgcTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTgg                                                                                                                                     <  1:784991/150‑1 (MQ=255)
                   gcTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGa                                                                                                                                    <  1:1714537/150‑1 (MQ=255)
                     ttgcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtg                                                                                                                                  >  2:5693873/1‑150 (MQ=255)
                       gcccttgccGGCGCGGTATTGGCGAATCTGGAACACCTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGGTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                                <  1:5956746/150‑1 (MQ=255)
                       gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                                <  2:1532736/150‑1 (MQ=255)
                       gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                                <  2:2581448/150‑1 (MQ=255)
                       gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                                <  2:6602401/150‑1 (MQ=255)
                         ccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgttgt                                                                                                                              >  2:3896771/1‑150 (MQ=255)
                           ttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGc                                                                                                                            >  1:6787399/1‑150 (MQ=255)
                               cGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCa                                                                                                                        >  2:4432672/1‑150 (MQ=255)
                                ggCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCAc                                                                                                                       >  2:6777759/1‑150 (MQ=255)
                                     ggTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCaa                                                                                                                  >  1:6350811/1‑150 (MQ=255)
                                      gTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAAt                                                                                                                 <  2:6286910/150‑1 (MQ=255)
                                       tATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCCGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAAtt                                                                                                                >  2:6026153/1‑150 (MQ=255)
                                          tGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCa                                                                                                             <  1:4407432/150‑1 (MQ=255)
                                           ggCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCAc                                                                                                            <  2:2145301/150‑1 (MQ=255)
                                            gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGAGCTCCACGCCAATATCACg                                                                                                           >  2:6207670/1‑150 (MQ=255)
                                              gAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGaa                                                                                                         >  2:5505893/1‑150 (MQ=255)
                                                aTCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATc                                                                                                       <  2:558891/150‑1 (MQ=255)
                                                   tGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTc                                                                                                    <  1:6070598/150‑1 (MQ=255)
                                                   tGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTc                                                                                                    >  1:3739738/1‑150 (MQ=255)
                                                    ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCa                                                                                                   >  1:4328326/1‑150 (MQ=255)
                                                    ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCa                                                                                                   >  1:4011607/1‑150 (MQ=255)
                                                    ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCa                                                                                                   >  1:2766481/1‑150 (MQ=255)
                                                    ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCa                                                                                                   >  2:32/1‑150 (MQ=255)
                                                     gAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCaa                                                                                                  <  2:3050456/150‑1 (MQ=255)
                                                      aaCACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAAc                                                                                                 <  1:1255961/150‑1 (MQ=255)
                                                       acacATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACg                                                                                                >  1:2404457/1‑150 (MQ=255)
                                                        cacaGACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGc                                                                                               <  1:4037965/150‑1 (MQ=255)
                                                         acaTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGcc                                                                                              >  1:5571755/1‑150 (MQ=255)
                                                          caTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCa                                                                                             <  2:4592822/150‑1 (MQ=255)
                                                            tACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCGGCGGGGATAACGTGATAAGGTACGGGGGAGTAATATATCGCCCGCGCCAGCGCCGCCGATTCGAGGATGTGGAGGTTTTGCTCCACGCCAATTACACGAAGATACAACGCAATg                                                                                           >  2:2264019/1‑150 (MQ=255)
                                                               gCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGaaa                                                                                        >  2:277545/1‑150 (MQ=255)
                                                                ccAGCACCTGAGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAATTTCGGTGGAGTAATAGATCGACCGCGCCAGCGCCGGCGATACGCGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAAt                                                                                       >  2:345334/1‑150 (MQ=255)
                                                                     aCCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTa                                                                                  <  1:2037227/150‑1 (MQ=255)
                                                                     aCCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTa                                                                                  <  2:1735528/150‑1 (MQ=255)
                                                                            ccACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTgg                                                                           >  2:3714720/1‑150 (MQ=255)
                                                                              aCCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGcc                                                                         >  1:4964016/1‑150 (MQ=255)
                                                                                  ccAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCcagca                                                                     <  1:5288049/150‑1 (MQ=255)
                                                                                   cAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGTAATCGCTACCCTTGGCcagcag                                                                    >  1:2183412/1‑150 (MQ=255)
                                                                                     ggTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAgcagca                                                                  <  2:5781209/150‑1 (MQ=255)
                                                                                      gTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCaa                                                                 >  1:3814782/1‑150 (MQ=255)
                                                                                        aCAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACg                                                               >  2:1903069/1‑150 (MQ=255)
                                                                                              ccgccgGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAt                                                                                                                                                 <  2:5043271/62‑1 (MQ=25)
                                                                                                 ccgGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCAc                                                      >  2:6652661/1‑150 (MQ=255)
                                                                                                    ggATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGc                                                   <  1:1816105/150‑1 (MQ=255)
                                                                                                    ggATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGc                                                   >  1:2721137/1‑150 (MQ=255)
                                                                                                     gATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGcc                                                  >  1:5856568/1‑150 (MQ=255)
                                                                                                         tcCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTtc                                              >  1:912292/1‑150 (MQ=255)
                                                                                                         tcCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTtc                                              >  1:2965213/1‑150 (MQ=255)
                                                                                                            tGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTcc                                           <  2:46871/150‑1 (MQ=255)
                                                                                                             gCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCg                                          <  2:3976692/150‑1 (MQ=255)
                                                                                                              cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATc                                                                                                       <  2:3597999/88‑1 (MQ=255)
                                                                                                              cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                                         >  1:6173332/1‑150 (MQ=255)
                                                                                                              cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                                         >  2:2479278/1‑150 (MQ=255)
                                                                                                              cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                                         <  1:4528450/150‑1 (MQ=255)
                                                                                                               tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                        >  2:5126861/1‑150 (MQ=255)
                                                                                                               tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                        <  1:6443951/150‑1 (MQ=255)
                                                                                                               tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                        <  1:4695554/150‑1 (MQ=255)
                                                                                                                cAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGt                                       <  1:3913943/150‑1 (MQ=255)
                                                                                                                 aaGTTCGGTGGAGTAATAGATAGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCATCAAAGCCATGAAATAGCAACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTc                                      >  2:4295145/1‑150 (MQ=255)
                                                                                                                       ggTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTg                                <  1:546868/150‑1 (MQ=255)
                                                                                                                        gTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGc                               >  1:1600808/1‑150 (MQ=255)
                                                                                                                        gTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGc                               <  2:6451317/150‑1 (MQ=255)
                                                                                                                           gAGTAATAGATCGCCCGCGCCAGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                            <  1:6140042/150‑1 (MQ=255)
                                                                                                                           gAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                            <  2:931613/150‑1 (MQ=255)
                                                                                                                           gAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                            >  2:5557446/1‑150 (MQ=255)
                                                                                                                           gAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                            <  2:3579651/150‑1 (MQ=255)
                                                                                                                              tAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGcc                         >  2:4865926/1‑150 (MQ=255)
                                                                                                                              tAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGcc                         >  1:3448524/1‑150 (MQ=255)
                                                                                                                              cctcgaaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGcc                         >  1:1468607/7‑150 (MQ=255)
                                                                                                                               ctcgaaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa                        >  2:3351478/6‑150 (MQ=255)
                                                                                                                               aaTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAACTCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa                        <  1:5726373/150‑1 (MQ=255)
                                                                                                                                aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCAc                       >  1:6332224/1‑150 (MQ=255)
                                                                                                                                 tAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCAcc                      <  2:1857695/150‑1 (MQ=255)
                                                                                                                                   aaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGc                    >  1:516102/2‑150 (MQ=255)
                                                                                                                                   aaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCCCCTTTCTCCGGGTCGTATTGCAAAGCCACCGc                    >  1:2561022/2‑150 (MQ=255)
                                                                                                                                    aTCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTc                                                                                                              <  1:2142455/59‑1 (MQ=255)
                                                                                                                                       gCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATaa                <  1:1405558/150‑1 (MQ=255)
                                                                                                                                       gCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCATGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTCTTGCAAAGCCACCGCATaa                >  1:5386759/1‑150 (MQ=255)
                                                                                                                                        cccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACgg                                                              <  1:3222132/103‑1 (MQ=255)
                                                                                                                                         ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATg              >  1:1044948/1‑150 (MQ=255)
                                                                                                                                         ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATg              >  1:3474344/1‑150 (MQ=255)
                                                                                                                                           gcgcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCg            >  1:5055499/1‑150 (MQ=255)
                                                                                                                                             gcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTc          <  2:4960427/150‑1 (MQ=255)
                                                                                                                                              ccAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCg         <  2:2952906/150‑1 (MQ=255)
                                                                                                                                               cAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCgg        >  1:1530005/1‑150 (MQ=255)
                                                                                                                                                 gcgcCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGGt      <  1:246283/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGGTTGGT  >  NC_002947/4950472‑4950770

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: