Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,038,228 G→A 7.2% A389A (GCC→GCT PP_0897 ← class 1 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,038,2280GA 7.2% 244.0 / 5.9 84A389A (GCC→GCTPP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base A (3/3);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

AACGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATCT  >  NC_002947/1038081‑1038373
                                                                                                                                                        |                                                                                                                                                 
aaCGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCgg                                                                                                                                                  >  6:139707/1‑149 (MQ=255)
      cATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGc                                                                                                                                            <  4:210177/149‑1 (MQ=255)
         cAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTa                                                                                                                                         >  6:168833/1‑149 (MQ=255)
          aGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTa                                                                                                                                         <  6:251726/148‑1 (MQ=255)
           ggTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTCCAAATCATGCCCAGCAGGCG‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACATGTACATTCGGGAAGCAGTGGTCGGGCCGGCTGGGCCTAc                                                                                                                                        >  8:46390/1‑148 (MQ=255)
               cacGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTAGTCCCAGCCGGCCCGACCACt                                                                                                                                    <  8:146816/148‑1 (MQ=255)
               cacGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTAGGCCCAGCCGGCCCGACCACt                                                                                                                                    <  2:178732/148‑1 (MQ=255)
                acGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTc                                                                                                                                  >  4:279676/1‑149 (MQ=255)
                acGGCCTTGTTGTCCTTGATCGCCTCGATAGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTc                                                                                                                                  >  1:21156/1‑149 (MQ=255)
                 cGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTc                                                                                                                                  >  8:259411/1‑148 (MQ=255)
                 cGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTc                                                                                                                                  <  7:259411/148‑1 (MQ=255)
                 cGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCg                                                                                                                                 >  5:360918/1‑149 (MQ=255)
                  ggCCTTGTTGTCCTTGATCGCCTCGATAGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGt                                                                                                                                <  2:21156/149‑1 (MQ=255)
                     cTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTAGGCCCAGCCGGCCCGACCACtgctaccc                                                                                                                             <  7:69975/149‑8 (MQ=38)
                          tgtTCTTTATCGCCTCGATGGCGGTGGGGCCAAGCTCGGACTTGTCAATAATTCCTTCCTTTTCTTTTTTTTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTca                                                                                                                              <  1:311262/148‑1 (MQ=255)
                           gtCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCg                                                                                                                       >  5:155996/1‑149 (MQ=255)
                            tCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCgg                                                                                                                      >  1:9252/1‑149 (MQ=255)
                             ccTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCggg                                                                                                                     >  8:42124/1‑149 (MQ=255)
                               ttGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGAC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTc                                                                                                                   >  3:307181/1‑149 (MQ=255)
                                 gATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTc                                                                                                                   <  5:214924/147‑1 (MQ=255)
                                   tCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGAcc                                                                                                               >  1:37625/1‑149 (MQ=255)
                                      ccTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCggg                                                                                                            >  5:111557/1‑149 (MQ=255)
                                      ccTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCggg                                                                                                            >  3:2456/1‑149 (MQ=255)
                                      ccTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCAACCACCGGGTCGACCggg                                                                                                            >  3:39956/1‑149 (MQ=255)
                                       cTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGc                                                                                                           >  4:312417/1‑149 (MQ=255)
                                        tCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGcc                                                                                                          >  4:206685/1‑149 (MQ=255)
                                              gCGGTCGGGCCACGCTCGGAGTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCTTACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACgta                                                                                                    <  4:296207/149‑1 (MQ=255)
                                               cGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACgtag                                                                                                   <  6:270965/149‑1 (MQ=255)
                                                 gTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACgtagta                                                                                                 <  7:324441/149‑1 (MQ=255)
                                                      gCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATg                                                                                             >  2:238521/1‑148 (MQ=255)
                                                           gCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCgg                                                                                       <  5:168833/149‑1 (MQ=255)
                                                                gACTTGCCAATCATTCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGAGCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                  <  3:351428/149‑1 (MQ=255)
                                                                   ttGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                               >  4:80534/1‑149 (MQ=255)
                                                                     gCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAgg                                                                              <  8:270722/148‑1 (MQ=255)
                                                                     gCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAgg                                                                              >  8:135716/1‑148 (MQ=255)
                                                                     gCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                               <  3:279676/147‑1 (MQ=255)
                                                                     gCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                               <  5:139707/147‑1 (MQ=255)
                                                                        aaTCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTa                                                                          <  2:9252/149‑1 (MQ=255)
                                                                         aTCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGaa                                                                                           >  5:102420/1‑131 (MQ=255)
                                                                         aTCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGaa                                                                                           <  6:102420/131‑1 (MQ=255)
                                                                         aTCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTAc                                                                         <  6:208675/149‑1 (MQ=255)
                                                                            aTGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCg                                                                       <  3:80534/148‑1 (MQ=255)
                                                                                ccAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCGTACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGc                                                                                                           >  5:209580/1‑108 (MQ=255)
                                                                                ccAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCGTACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGc                                                                                                           <  6:209580/108‑1 (MQ=255)
                                                                                ccAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCaa                                                                   <  7:135716/148‑1 (MQ=255)
                                                                                     agGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                             <  6:360918/149‑1 (MQ=255)
                                                                                     agGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                             <  3:312417/149‑1 (MQ=255)
                                                                                     agGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                             <  6:179089/149‑1 (MQ=255)
                                                                                      ggCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCg                                                            <  7:42124/149‑1 (MQ=255)
                                                                                      ggCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCg                                                            <  6:132487/149‑1 (MQ=255)
                                                                                         c‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTc                                                                                                                               >  7:6104/1‑79 (MQ=255)
                                                                                         c‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTc                                                                                                                               <  8:6104/79‑1 (MQ=255)
                                                                                                   tGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGtt                                                    >  5:25256/1‑149 (MQ=255)
                                                                                                     cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGtt                                                    <  6:111557/147‑1 (MQ=255)
                                                                                                     cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGa                                                  <  2:37625/149‑1 (MQ=255)
                                                                                                      tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                               <  2:349613/119‑1 (MQ=255)
                                                                                                      tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                               >  1:349613/1‑119 (MQ=255)
                                                                                                      tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGTCAACTCTTCGCCACGGTTGAg                                                 >  8:253040/1‑149 (MQ=255)
                                                                                                           gATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                               >  3:196124/1‑114 (MQ=255)
                                                                                                           gATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                               <  4:196124/114‑1 (MQ=255)
                                                                                                           gATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCAtt                                             >  1:219482/1‑148 (MQ=255)
                                                                                                               tGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAAc                                        >  8:33576/1‑149 (MQ=255)
                                                                                                                ggCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAAc                                        >  5:249493/1‑148 (MQ=255)
                                                                                                                 gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCa                                          <  4:147174/145‑1 (MQ=255)
                                                                                                                   gAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATg                                    >  2:77166/1‑149 (MQ=255)
                                                                                                                   gAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATg                                    >  3:154618/1‑149 (MQ=255)
                                                                                                                   gAGTGAACTTGTCCATACGGGTAGCAGTCGACGGGACGGCTGGGCATACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTACTAGATGAAGCCGCCTTTCAGCTCTACCGGCAACTCTTCGCCACGGATGAGCATTTCAACCATg                                    >  5:299918/1‑149 (MQ=255)
                                                                                                                       gAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCt                                >  2:271922/1‑149 (MQ=255)
                                                                                                                        aaCTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCtt                               <  1:238521/149‑1 (MQ=255)
                                                                                                                          cTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTtgt                             <  6:155996/149‑1 (MQ=255)
                                                                                                                           ttGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTAcc                                                                        >  8:130525/1‑105 (MQ=255)
                                                                                                                           ttGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTAcc                                                                        <  7:130525/105‑1 (MQ=255)
                                                                                                                           ttGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTtgtg                            >  1:130001/1‑149 (MQ=255)
                                                                                                                              tCCATAAGTGTAGCAGTGGTCGGTCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTgcg                          <  7:319541/148‑1 (MQ=255)
                                                                                                                                cATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGcgcg                        <  1:271922/148‑1 (MQ=255)
                                                                                                                                      tgacgaaGTGGTAGGCCCAGCCGGCCCGACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCAc                 >  7:214499/7‑149 (MQ=255)
                                                                                                                                      tgacgaaGTGGTAGGCCCAGCCGGCCCGACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCAc                 >  8:42703/7‑149 (MQ=255)
                                                                                                                                       gTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACg                <  4:2456/149‑1 (MQ=255)
                                                                                                                                           cAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCg            <  4:154618/149‑1 (MQ=255)
                                                                                                                                           aaGTGGTAGGCCCAGCCGGCCCGACCCCTTCGTCCCCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGGTGAGCATTTCAACCATGCGCTTGTGCGCGGCTTCACGGCCg            >  7:247704/2‑149 (MQ=255)
                                                                                                                                              tGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCggac         >  4:360115/1‑147 (MQ=255)
                                                                                                                                                gTCGGGCCTGCTGTGACTACCACTTCGTCACCCACCGGGTCGACCGGTCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCa        <  4:39956/148‑1 (MQ=255)
                                                                                                                                                gTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCg            <  1:29160/144‑1 (MQ=255)
                                                                                                                                                gTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCg            >  2:29160/1‑144 (MQ=255)
                                                                                                                                                     gCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATCt  >  2:84776/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
AACGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCC‑‑‑‑‑AGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCTACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATCT  >  NC_002947/1038081‑1038373

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: