Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,574,791 A→G 7.1% F168S (TTC→TCC) ‡ motB ← flagellar motor rotation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,574,7910AG 7.1% 200.7 / 4.6 70F168S (TTC→TCC) ‡motBflagellar motor rotation protein
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base G (3/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

TTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGA  >  NC_002947/5574647‑5574931
                                                                                                                                                |                                                                                                                                            
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                                                                                                                               >  6:334789/1‑149 (MQ=255)
           gCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTtct                                                                                                                               >  2:308205/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                                                                                                                      <  2:186793/140‑1 (MQ=255)
                             cTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCa                                                                                                                      >  1:186793/1‑140 (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                                                                                <  7:9883/142‑1 (MQ=255)
                                                                 cGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAAc                                                                                >  8:9883/1‑142 (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                                                                                   <  5:114619/103‑1 (MQ=255)
                                                                                                     gATATCTTCGAAGTAGGGCTGCAGGTGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCg                                                                                   >  6:114619/1‑103 (MQ=255)
                                                                                                       tatCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCAcc                                   >  7:76932/1‑149 (MQ=255)
                                                                                                                    gggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCATCTGCGGGTTCTCTTCCACCTTGTTCTGCAGt                       >  2:172533/1‑148 (MQ=255)
                                                                                                                     ggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTc                     <  7:119267/149‑1 (MQ=255)
                                                                                                                     ggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTc                     <  3:201571/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                           >  8:352767/6‑131 (MQ=255)
                                                                                                                                 aaccgaCCGATGTTCGACATCGGTAGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGGACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTg                           <  7:352767/126‑1 (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)
                                                                                                                                                |                                                                                                                                            
TTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGA  >  NC_002947/5574647‑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: