Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,574,789 C→T 8.1% D169N (GAC→AAC)  motB ← flagellar motor rotation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,574,7890CT 8.1% 215.3 / 6.0 74D169N (GAC→AAC) motBflagellar motor rotation protein
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base T (2/4);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTC  >  NC_002947/5574650‑5574927
                                                                                                                                           |                                                                                                                                          
ccGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAtcg                                                                                                                                    <  3:132156/148‑1 (MQ=255)
 cGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAtcggt                                                                                                                                  >  5:147821/1‑149 (MQ=255)
  gTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTTCGACAtcggta                                                                                                                                 <  5:301077/149‑2 (MQ=255)
  gTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTTCGACAtcggta                                                                                                                                 <  7:242584/149‑2 (MQ=255)
       cGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTTCGACAtcggtagcgc                                                                                                                             <  7:321857/148‑6 (MQ=255)
        gCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTAGGTtct                                                                                                                           >  7:245493/1‑149 (MQ=255)
               ccGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCAt                                                                                                                     <  2:128652/148‑1 (MQ=255)
                cGGCATTCGGCTTGGCATCGGTATGGCCGCTTATGCTGACCTTGTTCGTCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATcc                                                                                                                   <  8:245493/149‑1 (MQ=255)
                cGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATc                                                                                                                    >  4:66887/1‑148 (MQ=255)
                        ggCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTg                                                                                                           >  8:103940/1‑149 (MQ=255)
                              gCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACg                                                                                                     >  2:3112/1‑149 (MQ=255)
                                aTCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCAt                                                                                                                 >  7:195991/1‑135 (MQ=255)
                                aTCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCAt                                                                                                                 <  8:195991/135‑1 (MQ=255)
                                   ggTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGc                                                                                                <  6:147821/149‑1 (MQ=255)
                                           cGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                        <  1:313973/149‑1 (MQ=255)
                                           cGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                        <  1:3112/149‑1 (MQ=255)
                                                 tGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGAtc                                                                                  <  7:103940/149‑1 (MQ=255)
                                                    tGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGAtctc                                                                                >  4:319498/1‑148 (MQ=255)
                                                     gACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCg                                                                               >  7:208592/1‑148 (MQ=255)
                                                     gACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCg                                                                               >  6:193644/1‑148 (MQ=255)
                                                         ttGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAg                                                                          >  7:167791/1‑149 (MQ=255)
                                                                gCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTgg                                                                   <  3:66887/149‑1 (MQ=255)
                                                                gCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTgg                                                                   >  7:1235/1‑149 (MQ=255)
                                                                    cGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGaa                                                                             >  6:332871/1‑135 (MQ=255)
                                                                      ccTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCtt                                                              >  5:207734/1‑148 (MQ=255)
                                                                          gATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACt                                                         >  1:325589/1‑149 (MQ=255)
                                                                           aTGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGAtc                                                                                  <  5:175202/123‑1 (MQ=255)
                                                                           aTGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGAtc                                                                                  >  6:175202/1‑123 (MQ=255)
                                                                            tGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCt                                                                     <  5:60475/135‑1 (MQ=255)
                                                                            tGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCt                                                                     >  6:60475/1‑135 (MQ=255)
                                                                                gtCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                        >  8:81454/1‑112 (MQ=255)
                                                                                gtCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                        <  7:81454/112‑1 (MQ=255)
                                                                                 tCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGTCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCtt                                                              <  2:81771/137‑1 (MQ=255)
                                                                                 tCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGTCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCtt                                                              <  2:76862/137‑1 (MQ=255)
                                                                                 tCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGTCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCtt                                                              >  1:81771/1‑137 (MQ=255)
                                                                                 tCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGTCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCtt                                                              >  1:76862/1‑137 (MQ=255)
                                                                                 tCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCa                                                   <  8:208592/148‑1 (MQ=255)
                                                                                     ccATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGc                                              >  5:259053/1‑149 (MQ=255)
                                                                                     ccATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGc                                              <  2:325589/149‑1 (MQ=255)
                                                                                     ccATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGc                                              >  2:168546/1‑149 (MQ=255)
                                                                                       aTGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACa                                                                           <  4:246034/118‑1 (MQ=255)
                                                                                       aTGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACa                                                                           >  3:246034/1‑118 (MQ=255)
                                                                                        tGGCCAGTAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                        <  7:296338/104‑1 (MQ=255)
                                                                                        tGGCCAGTAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                        >  8:296338/1‑104 (MQ=255)
                                                                                                  gATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATGCGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCa                                 >  5:48370/1‑149 (MQ=255)
                                                                                                  gATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCa                                 <  8:1235/149‑1 (MQ=255)
                                                                                                        ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCgg                                                                                                                        <  1:189963/56‑1 (MQ=255)
                                                                                                        ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCgg                                                                                                                        >  2:189963/1‑56 (MQ=255)
                                                                                                        ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGtt                                         <  3:82960/135‑1 (MQ=255)
                                                                                                        ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGtt                                         >  4:82960/1‑135 (MQ=255)
                                                                                                        ttCGAAGTAGGGCTGCAGGCGCGCGCTAACGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGt                           >  2:248974/1‑149 (MQ=255)
                                                                                                         tCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                          >  8:9670/1‑149 (MQ=255)
                                                                                                          cGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                          <  6:259053/148‑1 (MQ=255)
                                                                                                           gAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTc                         <  1:168546/148‑1 (MQ=255)
                                                                                                           gAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTc                         <  7:9670/148‑1 (MQ=255)
                                                                                                            aaGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTg                       <  8:167791/149‑1 (MQ=255)
                                                                                                              gTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCa                     <  5:193644/149‑1 (MQ=255)
                                                                                                                aGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGt                   >  6:70374/1‑149 (MQ=255)
                                                                                                                 gggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGt                   >  4:251779/1‑148 (MQ=255)
                                                                                                                  ggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCTGTTc                 <  1:248974/149‑1 (MQ=255)
                                                                                                                   gCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCggg                                           >  1:65632/1‑122 (MQ=255)
                                                                                                                   gCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCggg                                           <  2:65632/122‑1 (MQ=255)
                                                                                                                       cAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgca            <  5:70374/149‑1 (MQ=255)
                                                                                                                       cAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgc             <  3:251779/148‑1 (MQ=255)
                                                                                                                         ggCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgcag           >  5:86074/1‑148 (MQ=255)
                                                                                                                            gcgcgcTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAAc       <  6:207734/149‑1 (MQ=255)
                                                                                                                              gcgcTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGc                                                    <  3:26875/102‑1 (MQ=255)
                                                                                                                              gcgcTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGc                                                    >  4:26875/1‑102 (MQ=255)
                                                                                                                                gcTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCt   >  2:317210/1‑149 (MQ=255)
                                                                                                                                gcTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCt   >  6:192702/1‑149 (MQ=255)
                                                                                                                                 cgaCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTc  >  3:330779/3‑149 (MQ=255)
                                                                                                                                  tACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTc  >  6:45260/1‑148 (MQ=255)
                                                                                                                                  gaCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgca            >  5:72800/2‑138 (MQ=255)
                                                                                                                                  gaCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgca            <  6:72800/137‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTC  >  NC_002947/5574650‑5574927

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: