Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,574,789 C→T 6.7% 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 6.7% 220.1 / 3.1 74D169N (GAC→AAC) motBflagellar motor rotation protein
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base T (3/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

CCCAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGA  >  NC_002947/5574642‑5574931
                                                                                                                                                   |                                                                                                                                              
cccAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTc                                                                                                                                                >  8:328429/1‑148 (MQ=255)
 ccAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGCTATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCg                                                                                                                                               >  4:359078/1‑148 (MQ=255)
 ccAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTc                                                                                                                                                <  3:156039/147‑1 (MQ=255)
 ccAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTc                                                                                                                                                <  3:150263/147‑1 (MQ=255)
 ccAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCg                                                                                                                                               >  7:29782/1‑148 (MQ=255)
 ccAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCg                                                                                                                                               <  7:134583/148‑1 (MQ=255)
 ccAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCg                                                                                                                                               >  1:345485/1‑148 (MQ=255)
     ttACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAt                                                                                                                                          <  8:214063/149‑1 (MQ=255)
     ttACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAt                                                                                                                                          >  8:331380/1‑149 (MQ=255)
       aCCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAtcg                                                                                                                                        >  4:201571/1‑149 (MQ=255)
        ccGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTTCGACAtcgg                                                                                                                                       <  5:377324/149‑1 (MQ=255)
           tATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAtcggtcg                                                                                                                                    >  5:109275/1‑149 (MQ=255)
           tATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAtcggtcg                                                                                                                                    >  3:133/1‑149 (MQ=255)
              tCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGtt                                                                                                                                 >  1:269617/1‑149 (MQ=255)
                gCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTtct                                                                                                                               >  2:308205/1‑149 (MQ=255)
                gCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTtct                                                                                                                               >  6:334789/1‑149 (MQ=255)
                  cAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCg                                                                                                                             <  3:122157/149‑1 (MQ=255)
                         ggCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCa                                                                                                                      <  3:131555/149‑1 (MQ=255)
                          gCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCAt                                                                                                                     <  2:201243/149‑1 (MQ=255)
                               cGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATc                                                                                                                 <  2:148205/148‑1 (MQ=255)
                                ggCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATc                                                                                                                        >  1:344487/1‑140 (MQ=255)
                                ggCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATc                                                                                                                        <  2:344487/140‑1 (MQ=255)
                                  gttGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGa                                                                                                             <  8:32548/148‑1 (MQ=255)
                                  cTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATcggtcggt                                                                                                                                  <  5:113657/128‑1 (MQ=255)
                                  cTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATcggtcggt                                                                                                                                  >  6:113657/1‑128 (MQ=255)
                                  cTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCa                                                                                                                      >  1:186793/1‑140 (MQ=255)
                                  cTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCa                                                                                                                      <  2:186793/140‑1 (MQ=255)
                                     ggCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGAt                                                                                                            <  5:334789/147‑1 (MQ=255)
                                     ggCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATAc                                                                                                          >  5:295433/1‑149 (MQ=255)
                                       cATCGCTATGCCCGATGAAGCAGACCTTGTTCGGCACCGCCTAGATGGTGTCAGCCATGGCCAGCAGGATAACTTCGAAGTAGGGCTTCAGGCGCGCGCTACCGATGTCGAACTTCGGTCGGTTCTCGTCATCCATGATCTGGATACg                                                                                                         >  2:311420/1‑148 (MQ=255)
                                         tCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGTGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAg                                                                                                      <  3:261803/149‑1 (MQ=255)
                                           ggTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGc                                                                                                    >  8:143364/1‑149 (MQ=255)
                                           ggTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGc                                                                                                    <  2:127227/149‑1 (MQ=255)
                                            gTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGcc                                                                                                   >  8:65525/1‑149 (MQ=255)
                                                 gCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCt                                                                                              <  2:345485/149‑1 (MQ=255)
                                                 gCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCt                                                                                              <  8:29782/149‑1 (MQ=255)
                                                   cGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                            <  7:331380/149‑1 (MQ=255)
                                                    gCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAg                                                                                           >  4:122019/1‑149 (MQ=255)
                                                           cTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGAtctc                                                                                    <  6:109275/149‑1 (MQ=255)
                                                                      cGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAAc                                                                                >  8:9883/1‑142 (MQ=255)
                                                                      cGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAAc                                                                                <  7:9883/142‑1 (MQ=255)
                                                                         cACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGt                                                                      <  7:65525/149‑1 (MQ=255)
                                                                          aCCGCCTTGATGGTGTCAGCCATGGCCAGCAGGAGATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTc                                                                     <  8:60933/149‑1 (MQ=255)
                                                                           ccGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTc                                                                                                <  1:207972/121‑1 (MQ=255)
                                                                           ccGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTc                                                                                                >  2:207972/1‑121 (MQ=255)
                                                                             gCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCtt                                                                  <  1:30748/149‑1 (MQ=255)
                                                                                  gATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACt                                                             >  7:294532/1‑149 (MQ=255)
                                                                                    tGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTGGAACAGGATCTGGTCCTTGAACTTc                                                           >  4:294108/1‑149 (MQ=255)
                                                                                      gtgtCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTg                                                         <  2:269617/149‑1 (MQ=255)
                                                                                         tCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCa                                                       <  4:133/148‑1 (MQ=255)
                                                                                         tCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAg                                                      >  8:115339/1‑149 (MQ=255)
                                                                                           aGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGAGCTCGAACAGGATCTCGTCCTTGAACTTATGCAGCt                                                    >  7:329214/1‑149 (MQ=255)
                                                                                            gCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTg                                                   >  5:134305/1‑149 (MQ=255)
                                                                                               aTGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTcc                                                                    <  6:190912/129‑1 (MQ=255)
                                                                                                        aggttctCTTCGAAGTATGGCTGCAGGCGCGCGCTACCGATGTCGAACATGGGTCGGTTCTCGCCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTc                                       <  8:329214/143‑1 (MQ=255)
                                                                                                          gATATCTTCGAAGTAGGGCTGCAGGTGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCg                                                                                   >  6:114619/1‑103 (MQ=255)
                                                                                                          gATATCTTCGAAGTAGGGCTGCAGGTGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCg                                                                                   <  5:114619/103‑1 (MQ=255)
                                                                                                            tatCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCAcc                                   >  7:76932/1‑149 (MQ=255)
                                                                                                                         gggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCATCTGCGGGTTCTCTTCCACCTTGTTCTGCAGt                       >  2:172533/1‑148 (MQ=255)
                                                                                                                          ggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTc                     <  3:201571/149‑1 (MQ=255)
                                                                                                                          ggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTc                     <  7:119267/149‑1 (MQ=255)
                                                                                                                           gCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTcc                    >  7:320051/1‑149 (MQ=255)
                                                                                                                               cAGGCGCGCGCTACCGATGTCGAACCTCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTg                  <  3:241358/147‑1 (MQ=255)
                                                                                                                               cAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgca                <  3:122019/149‑1 (MQ=255)
                                                                                                                               cAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgc                 <  3:359078/148‑1 (MQ=255)
                                                                                                                               cAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTAGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGGGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTg                  <  3:294108/147‑1 (MQ=255)
                                                                                                                                aGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgcag               <  6:134305/149‑1 (MQ=255)
                                                                                                                                    gcgcgcTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAAc           <  7:328429/149‑1 (MQ=255)
                                                                                                                                      aaccgaCCGATGTTCGACATCGGTAGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGGACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTg                           <  7:352767/126‑1 (MQ=255)
                                                                                                                                      aaccgaCCGATGTTCGACATCGGTAGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGGACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTg                           >  8:352767/6‑131 (MQ=255)
                                                                                                                                          tACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTc      >  2:135794/1‑148 (MQ=255)
                                                                                                                                          gaCCGATGTTCGACATCGGTCGGTTCACGGCATCCATGATCTGAATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCAGGTCCTTGAACTTCAGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTc      >  1:269125/2‑148 (MQ=255)
                                                                                                                                           aCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTc      >  6:329149/1‑147 (MQ=255)
                                                                                                                                            ccGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCa    <  1:308205/148‑1 (MQ=255)
                                                                                                                                             cGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGCACTTCTGCAGCTGCGGGTTCTCTTCCCCCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAgc  >  8:166844/1‑148 (MQ=255)
                                                                                                                                             cGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCAc                                    <  8:26385/115‑1 (MQ=255)
                                                                                                                                             cGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCAc                                    >  7:26385/1‑115 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCCAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGA  >  NC_002947/5574642‑5574931

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: