Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,948,990 C→A 6.6% G119V (GGC→GTC)  fliN ← flagellar basal‑body C‑ring protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,948,9900CA 6.6% 231.1 / 4.0 73G119V (GGC→GTC) fliNflagellar basal‑body C‑ring protein
Reads supporting (aligned to +/- strand):  ref base C (30/38);  new base A (2/3);  total (32/41)
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.82e-01

GCGCATCGTGCCCTTCACTCAGCGCAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCGGGCCTTCGAGGCTGACATTTTC  >  NC_002947/4948842‑4949138
                                                                                                                                                       |                                                                                                                                                     
gcgcATCGTGCCCTTCACTCAGCGCAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAAGTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGc                                                                                                                                                       >  2:261046/1‑149 (MQ=255)
gcgcATCGTGCCCTTCACTCAGCGCAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGc                                                                                                                                                       >  4:288181/1‑149 (MQ=255)
       gTGCCCTTCACTCAGCGCAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAg                                                                                                                                                <  7:155734/149‑1 (MQ=255)
            cTTCACTCAGCGCAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGt                                                                                                                                           >  1:213200/1‑149 (MQ=255)
             ttCACTCAGCGCAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTc                                                                                                                                          >  1:315389/1‑149 (MQ=255)
                aCTCAGCGCAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGa                                                                                                                                        <  8:31486/148‑1 (MQ=255)
                       gcAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGAt                                                                                                                                 >  7:130781/1‑148 (MQ=255)
                           cttcttGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGCGCTCGATCAcc                                                                                                                             >  8:245217/1‑148 (MQ=255)
                            ttcttGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCAcc                                                                                                                             >  1:70283/1‑147 (MQ=255)
                              cttGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGa                                                                                                                           <  5:37724/147‑1 (MQ=255)
                                 gATACGTTCGCTGGGGCTGATCACGTCAGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGcc                                                                                                                        <  7:37322/147‑1 (MQ=255)
                                   tACTTTCGCTGGGGCTGATGACGTAGTTCAGGCGTATGCCGTACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGAGTGCCGTTGACCAGCAC‑‑‑GTGGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGGTGATCGAGctcga                                                                                                                     <  5:243754/148‑5 (MQ=255)
                                   tACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTg                                                                                                                    <  7:178971/148‑1 (MQ=255)
                                      gTTCGCTGGGGCTGATCCCGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGtt                                                                                                                 >  3:138109/1‑148 (MQ=255)
                                      gTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGtt                                                                                                                 >  7:21472/1‑148 (MQ=255)
                                      gTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTg                                                                                                                <  3:256532/149‑1 (MQ=255)
                                           cTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTg                                                                                                           <  4:31186/149‑1 (MQ=255)
                                              gggCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCAcc                                                                                                                             >  6:175357/1‑129 (MQ=255)
                                              gggCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCAcc                                                                                                                             <  5:175357/129‑1 (MQ=255)
                                              gggCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTgca                                                                                                         <  2:315389/148‑1 (MQ=255)
                                                gCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTgcag                                                                                                        >  5:106012/1‑147 (MQ=255)
                                                gCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTgcag                                                                                                        >  4:142822/1‑147 (MQ=255)
                                                 cTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGGTCGACCGCCTGGCCGGCGAGc                                                                                                                                      >  5:328244/1‑117 (MQ=255)
                                                 cTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGACCGCCTGGCCGGCGAGc                                                                                                                                      <  6:328244/117‑1 (MQ=255)
                                                    aTCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAgg                                                                                                    <  2:213200/147‑1 (MQ=255)
                                                    aTCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGtt                                                                                                  >  3:286403/1‑149 (MQ=255)
                                                    aTCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGtt                                                                                                  >  2:134132/1‑149 (MQ=255)
                                                     tCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGtt                                                                                                                 <  3:16314/133‑1 (MQ=255)
                                                     tCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGtt                                                                                                                 >  4:16314/1‑133 (MQ=255)
                                                         gTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGGTGAGCTGCAGCAGGTTGCGGa                                                                                             <  8:100714/149‑1 (MQ=255)
                                                            ggTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGt                                                                                          <  7:69763/149‑1 (MQ=255)
                                                            ggTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGt                                                                                          <  8:130781/149‑1 (MQ=255)
                                                              tCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCgg                                                                                                                                                    <  2:145158/90‑1 (MQ=255)
                                                              tCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTg                                                                                        <  1:134132/149‑1 (MQ=255)
                                                                 ggCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGAtc                                                                                     <  3:142822/149‑1 (MQ=255)
                                                                         ccGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCAGCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGtgct                                                                             >  3:302494/1‑149 (MQ=255)
                                                                          cGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGc                                                                                                                                                   >  6:1769/1‑79 (MQ=255)
                                                                          cGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGc                                                                                                                                                   <  5:1769/79‑1 (MQ=255)
                                                                          cGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTgca                                                                                                         >  3:97621/1‑120 (MQ=255)
                                                                          cGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTgca                                                                                                         <  4:97621/120‑1 (MQ=255)
                                                                          cGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGtgctg                                                                            <  3:288181/149‑1 (MQ=255)
                                                                               ttttCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCAc                                                                       >  8:152348/1‑149 (MQ=255)
                                                                                tttCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACt                                                                      >  4:4166/1‑149 (MQ=255)
                                                                                  tagTTGACCACGAGCACTTCGCCTTGGGCGATCAGCGTGCAGTTGACCAGCCCTCTGTCGTCCCGCTCGACGGCCAGGCGGTGGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCTGTGCTGCCCa                                                                        <  1:105741/146‑1 (MQ=255)
                                                                                    gTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCa                                                                  >  2:123209/1‑149 (MQ=255)
                                                                                         ccacgaccacTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGtt                                                                                         <  3:54415/121‑1 (MQ=255)
                                                                                         ccacgaccacTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGtt                                                                                         >  4:54415/1‑121 (MQ=255)
                                                                                         ccacgaccacTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGa                                                              >  6:91715/1‑148 (MQ=255)
                                                                                          cacgaccacTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGaaa                                                            >  6:77162/1‑149 (MQ=255)
                                                                                              accacTTCGCCATCGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGtt                                                        <  2:45625/149‑1 (MQ=255)
                                                                                               ccacTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGtt                                                        >  4:218468/1‑148 (MQ=255)
                                                                                                  cTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTCTGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCAc                                                     <  1:261046/149‑1 (MQ=255)
                                                                                                  cTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCAcc                                                    <  7:152348/149‑1 (MQ=255)
                                                                                                  cTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCAcc                                                    >  1:302480/1‑149 (MQ=255)
                                                                                                   ttCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCg                                                   >  3:75266/1‑149 (MQ=255)
                                                                                                      gCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGt                                                                                >  4:277312/1‑117 (MQ=255)
                                                                                                      gCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGt                                                                                <  3:277312/117‑1 (MQ=255)
                                                                                                      gCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGTTGCCCACTTCCATGGAAATGTTCACCGGGa                                                <  4:138109/149‑1 (MQ=255)
                                                                                                            ggCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTcc                                           <  2:302480/148‑1 (MQ=255)
                                                                                                              cGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAgg                                        <  5:91715/149‑1 (MQ=255)
                                                                                                              cGATCAGCGGGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGATGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAgg                                        <  2:234170/149‑1 (MQ=255)
                                                                                                                 tCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGAt                                      >  4:176053/1‑148 (MQ=255)
                                                                                                                 tCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATc                                     <  4:155128/149‑1 (MQ=255)
                                                                                                                         ccGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAg                             <  4:286403/149‑1 (MQ=255)
                                                                                                                         ccGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAg                             <  8:21472/149‑1 (MQ=255)
                                                                                                                           gTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGt                           >  5:113328/1‑149 (MQ=255)
                                                                                                                           gTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGt                           <  2:70283/149‑1 (MQ=255)
                                                                                                                                      c‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCAGCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCGATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCGGGCCTTCg                <  4:302494/149‑1 (MQ=255)
                                                                                                                                           tCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCGGGCCTTCGAg              <  2:214872/149‑1 (MQ=255)
                                                                                                                                                gagCTCGACCGCCTGGCCGGCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCGGGCCTTCGAGGCt           >  3:214862/3‑147 (MQ=255)
                                                                                                                                                gagCTCGACCGCCTGGCCGGCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCGGGCCTTCGAGGCt           <  4:214862/145‑1 (MQ=255)
                                                                                                                                                  gCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCggg                      <  4:44785/134‑1 (MQ=255)
                                                                                                                                                  gCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCggg                      >  3:44785/1‑134 (MQ=255)
                                                                                                                                                   cTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCGGGCCTTCGAGGCTGACAt      <  6:106012/149‑1 (MQ=255)
                                                                                                                                                    tCGCCGGCCAGGCGGTAGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATCTAGATCTCGGTGCTGCCCAATTCCATGTAAATGAACACCCGGATGACCAGGATCACCTCCAGTTTCGGGCCTTCGAGGATGACAtt     >  6:236203/1‑149 (MQ=255)
                                                                                                                                                       ccGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCGGGCCTTCGAGGCTGACATTTTc  <  1:123209/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCGCATCGTGCCCTTCACTCAGCGCAGCTTCTTGATACGTTCGCTGGGGCTGATCACGTCGGTCAGGCGGATGCCGAACTTTTCGTTGACCACGACCACTTCGCCATGGGCGATCAGCGTGCCGTTGACCAGCAC‑‑‑GTCGAGCGGCTCGCCGGCCAGGCGGTCGAGCTCGATCACCGAGCCCTGGTTGAGCTGCAGCAGGTTGCGGATGTTGATCTCGGTGCTGCCCACTTCCATGGAAATGTTCACCGGGATGTCCAGGATCACGTCCAGGTTCGGGCCTTCGAGGCTGACATTTTC  >  NC_002947/4948842‑4949138

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: