Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,950,189 Δ1 bp 10.4% coding (155/969 nt) fliM ← flagellar biosynthesis switching/energizing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,950,1890G.10.4% 238.3 / 21.9 77coding (155/969 nt)fliMflagellar biosynthesis switching/energizing protein
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base . (4/4);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

GTACAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTCC  >  NC_002947/4950053‑4950332
                                                                                                                                        |                                                                                                                                                
gTACAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                     <  6:244068/149‑2 (MQ=255)
  aCAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATg                                                                                                                                             <  3:300115/140‑1 (MQ=255)
  aCAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATg                                                                                                                                             >  4:300115/1‑140 (MQ=255)
    aGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCAc                                                                                                                                 <  7:235957/149‑1 (MQ=255)
    aGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccc                                                                                                                                      <  4:209177/144‑1 (MQ=255)
    aGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccc                                                                                                                                      >  3:209177/1‑144 (MQ=255)
     gCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACg                                                                                                                                >  4:197242/1‑149 (MQ=255)
        aGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                     >  1:71672/1‑140 (MQ=255)
        aGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                     <  2:71672/141‑2 (MQ=255)
          tGCACGTATTCGTCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACg                                                                                                                           <  7:230690/149‑1 (MQ=255)
            cACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGt                                                                                                                         >  7:152286/1‑149 (MQ=255)
            cACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGt                                                                                                                         >  4:222657/1‑149 (MQ=255)
             aCGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTc                                                                                                                        >  4:67185/1‑149 (MQ=255)
             aCGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTc                                                                                                                        >  3:15773/1‑149 (MQ=255)
               gTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCt                                                                                                                      <  3:222657/149‑1 (MQ=255)
                 aTTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                    >  4:52020/1‑149 (MQ=255)
                 aTTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTAGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                    >  5:29158/1‑149 (MQ=255)
                  ttCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCa                                                                                                                                               >  1:107448/1‑122 (MQ=255)
                  ttCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCa                                                                                                                                               <  2:107448/122‑1 (MQ=255)
                  ttCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                    <  3:67185/148‑1 (MQ=255)
                     gCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACgg                                                                                                                          >  4:211691/1‑139 (MQ=255)
                     gCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACgg                                                                                                                          <  3:211691/139‑1 (MQ=255)
                     gCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTgg                                                                                                                >  1:225341/1‑149 (MQ=255)
                      ccGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGt                                                                                                               >  4:226514/1‑149 (MQ=255)
                      cagAACTTCATGACCTGCACGCCACCGACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGATGATACGGGTATAGCGGGCGAAGCGCCCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGt                                                                                                               <  1:250812/147‑1 (MQ=255)
                        gAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCa                                                                                                             <  3:197242/149‑1 (MQ=255)
                          aCTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAgg                                                                                                           <  3:226514/149‑1 (MQ=255)
                           cTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                    <  4:241280/139‑1 (MQ=255)
                           cTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGa                                                                                                                    >  3:241280/1‑139 (MQ=255)
                               aTGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCg                                                                                                        >  6:46229/1‑147 (MQ=255)
                                tGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACCAGTCAggaccgta                                                                                                     <  6:176014/149‑7 (MQ=255)
                                    cTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCttt                                                                                                 >  5:38945/1‑149 (MQ=255)
                                    cTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCATTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                     >  3:289149/1‑112 (MQ=255)
                                    cTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCATTGATCATTTCCAGGGTC‑GCATGCCGAccct                                                                                                                                     <  4:289149/113‑2 (MQ=255)
                                     tGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCtttt                                                                                                >  3:57463/1‑149 (MQ=255)
                                      gCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCtttt                                                                                                >  6:3805/1‑148 (MQ=255)
                                       cACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCtttt                                                                                                <  5:112289/147‑1 (MQ=255)
                                       cACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGa                                                                                              <  2:225341/149‑1 (MQ=255)
                                          gCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCt                                                                                                                      >  4:102076/1‑122 (MQ=255)
                                          gCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCt                                                                                                                      <  3:102076/122‑1 (MQ=255)
                                           ccacccacGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGGCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGGCAGGTCGTAGCTTTTGATACt                                                                                          <  3:220287/149‑1 (MQ=255)
                                            cacccacGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTg                                                                                         >  5:256324/1‑149 (MQ=255)
                                                    gCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTc                                                                                 >  7:57265/1‑149 (MQ=255)
                                                     ccACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCg                                                                                >  3:51537/1‑149 (MQ=255)
                                                            agcggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGc                                                                           <  6:240944/147‑1 (MQ=255)
                                                            agcggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGc                                                                           >  5:240944/1‑147 (MQ=255)
                                                              cggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGCCCCCGCACGAGACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCAAGGCTCGGATGCActc                                                                       >  5:292143/1‑149 (MQ=255)
                                                              cggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTg                                                                                         >  6:228995/1‑131 (MQ=255)
                                                              cggagcggCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTg                                                                                         <  5:228995/131‑1 (MQ=255)
                                                                          gcagGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCAc                                                            >  7:194872/1‑148 (MQ=255)
                                                                           cagGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCa                                                          <  2:196728/149‑1 (MQ=255)
                                                                                 gAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCcg                                                     <  8:57265/148‑1 (MQ=255)
                                                                                     aTGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATTCTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCgtcgtc                                                <  4:57463/149‑1 (MQ=255)
                                                                                      tGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGtcgtcg                                               <  6:38945/149‑1 (MQ=255)
                                                                                       gCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGa                                              >  8:84351/1‑149 (MQ=255)
                                                                                               gggTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGc                                      <  7:84351/149‑1 (MQ=255)
                                                                                                      gCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAggg                               <  4:15773/149‑1 (MQ=255)
                                                                                                       cGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAggg                               >  6:166073/1‑148 (MQ=255)
                                                                                                       cGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAgg                                <  8:247993/147‑1 (MQ=255)
                                                                                                       cGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAgg                                >  7:247993/1‑147 (MQ=255)
                                                                                                         ggCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCAt                            <  5:3805/149‑1 (MQ=255)
                                                                                                            gAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATg                                                                            >  8:128440/1‑98 (MQ=255)
                                                                                                            gAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATg                                                                            <  7:128440/98‑1 (MQ=255)
                                                                                                              aGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGAtt                       <  3:52020/149‑1 (MQ=255)
                                                                                                                      ttGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTg               <  4:51537/149‑1 (MQ=255)
                                                                                                                        gATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTgg              <  5:46229/148‑1 (MQ=255)
                                                                                                                        gATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTgg              >  4:109252/1‑148 (MQ=255)
                                                                                                                        gATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTgg              <  6:256324/148‑1 (MQ=255)
                                                                                                                          tCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGa            >  5:20373/1‑148 (MQ=255)
                                                                                                                          tCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGa            >  1:209210/1‑148 (MQ=255)
                                                                                                                          tCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGAc           <  8:194872/149‑1 (MQ=255)
                                                                                                                          tCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGCCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGa            >  2:264987/1‑148 (MQ=255)
                                                                                                                              ttCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGAcag         >  5:282796/1‑147 (MQ=255)
                                                                                                                                ccATGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCAc                                                            <  8:210627/94‑1 (MQ=255)
                                                                                                                                  ggggTC‑GCATGCCGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTc   >  4:63014/2‑149 (MQ=255)
                                                                                                                                  aGGTTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGCCGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTc   <  6:292143/149‑1 (MQ=255)
                                                                                                                                   gggTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCGACAGGGCATCGATTTCATCCTGGGACAGCAGGTcc  <  8:152286/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                
GTACAGCGAGTGCACGTATTCGCCGAACTTCATGACCTGCACGCCACCCACGGCCACGTCAGCGGAGCGGCGCAGCAGGTTGAACATGCTGATACGGGTATAGCGGGCGAAGCGCTCGTTGATCATTTCCAGGGTCGGCATG‑CGACCCCGCACGATACGGTCCTGACTGGTCAGGTCGTAGCTTTTGATACTGCCAGGCTCGGATGCACTCTCGGTCTGCACCAGCCCGTCGTCGACGCCATGCAACAGGGCATCGATTTCATCCTGGGACAGCAGGTCC  >  NC_002947/4950053‑4950332

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: ‑