Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,137,952 C→T 5.5% G194G (GGC→GGT)  flgJ → flagellar rod assembly protein and murein hydrolase; flagellum‑specific muramidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,137,9520CT 5.5% 233.1 / 2.4 72G194G (GGC→GGT) flgJflagellar rod assembly protein and murein hydrolase; flagellum‑specific muramidase
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base T (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

TGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGC  >  NC_000913/1137804‑1138095
                                                                                                                                                    |                                                                                                                                               
tGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGc                                                                                                                                                 >  4:149992/1‑149 (MQ=255)
         aCAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGcgagccgagc                                                                                                                                        >  3:338295/1‑149 (MQ=255)
           aGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCta                                                                                                                                      >  6:252814/1‑149 (MQ=255)
             tAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCtat                                                                                                                                     >  8:4435/1‑148 (MQ=255)
              aaaGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATaa                                                                                                                                   >  6:188353/1‑149 (MQ=255)
               aaGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATaa                                                                                                                                   <  8:331214/148‑1 (MQ=255)
                    ttCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGtt                                                                                                                             >  8:284770/1‑149 (MQ=255)
                     tCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGGAAGCCAGCAAAGCGGTGTGCCACATCATTTGCTCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCGGttt                                                                                                                            >  1:185129/1‑149 (MQ=255)
                      ccTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGttt                                                                                                                            <  1:201850/148‑1 (MQ=255)
                          gcgcAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGTGAGCCGAGCTATAACCTGTTTGgtgt                                                                                                                       >  1:4646/1‑149 (MQ=255)
                          gcgcAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGgtgt                                                                                                                       >  8:233338/1‑149 (MQ=255)
                                     gCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGTGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCctc                                                                                                              <  2:4646/147‑1 (MQ=255)
                                     gCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTg                                                                                                            >  7:150665/1‑149 (MQ=255)
                                         cccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCGCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCCACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCGATAACCTGGTAGGTGTCAAAGCCTCTGGCaa                                                                                                        >  1:270902/1‑149 (MQ=255)
                                          ccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAAc                                                                                                       >  2:46015/1‑149 (MQ=255)
                                           cgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACTGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACt                                                                                                      <  6:311182/149‑1 (MQ=255)
                                             cccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAAGCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTg                                                                                                     >  8:11050/1‑148 (MQ=255)
                                                aaCTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCATCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGa                                                                                                   >  3:67021/1‑147 (MQ=255)
                                                   tGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGc                                                                                                          >  6:143785/1‑137 (MQ=255)
                                                   tGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAggg                                                                                              >  5:242206/1‑149 (MQ=255)
                                                   tGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAggg                                                                                              >  8:383572/1‑149 (MQ=255)
                                                    ggCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCaa                                                                                                        >  1:179286/1‑138 (MQ=255)
                                                    ggCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCaa                                                                                                        <  2:179286/138‑1 (MQ=255)
                                                     gCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAAc                                                                                                       <  7:218515/138‑1 (MQ=255)
                                                     gCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAggg                                                                                              <  4:96688/147‑1 (MQ=255)
                                                     gCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGc                                                                                             <  7:233338/148‑1 (MQ=255)
                                                     gCAAGCCAGCAAAGCGGAGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAAc                                                                                                       >  8:218515/1‑138 (MQ=255)
                                                      cAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCa                                                                                           <  7:323581/149‑1 (MQ=255)
                                                       aaGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAg                                                                                          <  2:244807/149‑1 (MQ=255)
                                                          ccAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTa                                                                                       >  7:86391/1‑149 (MQ=255)
                                                            aGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACt                                                                                     <  5:188353/149‑1 (MQ=255)
                                                                 aGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGgt                                                                                                                         >  3:112691/1‑108 (MQ=255)
                                                                 aGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGgt                                                                                                                         <  4:112691/108‑1 (MQ=255)
                                                                 aGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAAt                                                                                <  8:150665/149‑1 (MQ=255)
                                                                  gCGGTGTGCCACATCGTTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATc                                                                               <  7:352843/149‑1 (MQ=255)
                                                                      tgtgCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATcacca                                                                           <  5:252814/149‑1 (MQ=255)
                                                                          ccACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGTGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGAc                                                                       >  7:36968/1‑149 (MQ=255)
                                                                           cacaTCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACt                                                                      >  5:868/1‑149 (MQ=255)
                                                                            acaTCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACt                                                                      >  3:280846/1‑148 (MQ=255)
                                                                            acaTCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACt                                                                      >  2:119546/1‑148 (MQ=255)
                                                                             catcatTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTg                                                                                                                           <  3:275784/94‑1 (MQ=255)
                                                                             catcatTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTg                                                                                                                           >  4:275784/1‑94 (MQ=255)
                                                                               tcatTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATACCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACGGaa                                                                   >  6:280746/1‑148 (MQ=255)
                                                                                 atTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTgg                                                                                                    <  8:166078/113‑1 (MQ=255)
                                                                                 atTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTgg                                                                                                    >  7:166078/1‑113 (MQ=255)
                                                                                 atTTGATCCTCGCGCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGCGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAATGGGCCAGTTACTGAAATCACCCCGACTGAatat                                                                >  3:57387/1‑149 (MQ=255)
                                                                                      aTCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGaaaa                                                           <  2:82790/149‑1 (MQ=255)
                                                                                          tCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGgc                                                       >  4:306057/1‑149 (MQ=255)
                                                                                            ggtCAGGCGGCACTGTAAGCTGGTGGGGGGAAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATCTGAAAACGgc                                                       <  2:270902/145‑1 (MQ=255)
                                                                                            gCTCAGGCGGCAATGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGGGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAACCGgc                                                       <  1:270898/147‑1 (MQ=255)
                                                                                             cTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAACCGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGgcg                                                      >  6:363159/1‑147 (MQ=255)
                                                                                             cTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGgcgaa                                                    <  3:65187/149‑1 (MQ=255)
                                                                                                aGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGgcgaagc                                                  <  6:868/148‑1 (MQ=255)
                                                                                                 ggcggcACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGgcgaagcga                                                <  5:179266/149‑1 (MQ=255)
                                                                                                  gcggcACTGGAATCTGGTTGGGGGCAACGGAAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTGGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGGTACTGAAATCACCACGACTGAATATGAAAACGgcgaagcgaa                                               <  1:197621/149‑1 (MQ=255)
                                                                                                        cTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATAACCACGACTGAATATGAAAACGGCGCAGCGAAGAAAg                                          >  7:168033/1‑148 (MQ=255)
                                                                                                         tGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGt                                         >  1:46559/1‑148 (MQ=255)
                                                                                                                  gTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGa                                 <  7:383572/147‑1 (MQ=255)
                                                                                                                      ggggCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTc                           >  2:254932/1‑149 (MQ=255)
                                                                                                                           aaCGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTc                      <  7:123330/149‑1 (MQ=255)
                                                                                                                             cGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTa                    <  7:284770/149‑1 (MQ=255)
                                                                                                                             cGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTa                    >  4:163947/1‑149 (MQ=255)
                                                                                                                                cAAATCCGCCGCGAAAACGGTGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAg                 <  8:36968/149‑1 (MQ=255)
                                                                                                                                 aaaTCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGc                >  7:182741/1‑149 (MQ=255)
                                                                                                                                   aTCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTc              <  3:163947/149‑1 (MQ=255)
                                                                                                                                       gccgcGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTAt          <  4:338295/149‑1 (MQ=255)
                                                                                                                                         cgcgAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACGACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCt        <  4:67021/149‑1 (MQ=255)
                                                                                                                                           cgAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTgg      <  1:46015/149‑1 (MQ=255)
                                                                                                                                             aaaaCGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCt        <  8:187464/145‑1 (MQ=255)
                                                                                                                                             aaaaCGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCt        >  7:187464/1‑145 (MQ=255)
                                                                                                                                              aaaCGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTgg      >  2:354616/1‑146 (MQ=255)
                                                                                                                                              aaaCGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTgg      <  1:354616/146‑1 (MQ=255)
                                                                                                                                              aaaCGGCGAGCAGAGCTATAACCTGTTTGGTGTCAAAGCCTCTTGCAACTGGAAAGGGCAAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTATACAGCTAGTATCTGGaac   >  3:191774/1‑148 (MQ=255)
                                                                                                                                               aaCGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGc  <  8:86391/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGC  >  NC_000913/1137804‑1138095

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