Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 362,372 G→C 5.9% D45E (GAC→GAG) ‡ PP_0301 ← betainyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947362,3720GC 5.9% 219.0 / 3.0 68D45E (GAC→GAG) ‡PP_0301betainyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base C (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GGCTGTGATACACGTGAAGACGCTTGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTGGTCTCCTTAA  >  NC_002947/362224‑362520
                                                                                                                                                       |                                                                                                                                                    
ggCTGTGATACACGTGAAGACGCTTGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCg                                                                                                                                                      <  8:37614/149‑1 (MQ=255)
ggCTGTGATACACGTGAAGACGCTTGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCg                                                                                                                                                      >  8:129697/1‑149 (MQ=255)
         acacGTGAAGACGCTTGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGcc                                                                                                                                              >  3:216417/1‑148 (MQ=255)
             gTGAAGACGCTTGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGc                                                                                                                                         >  4:233266/1‑149 (MQ=255)
                aaGACGCTTGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGt                                                                                                                                      <  3:233266/149‑1 (MQ=255)
                     gCTTGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAAATGAGATTAATATGCCCAACCAGCTTAAAAACCAATTAACCGATTTCCCCCCGGATAAGGGGGACCAGCCCGATGCggac                                                                                                                                     >  1:218787/1‑146 (MQ=255)
                       ttGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCAAt                                                                                                                                  >  6:38856/1‑146 (MQ=255)
                           gATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCCTGGACGCCGAtag                                                                                                                                         <  8:240671/135‑3 (MQ=255)
                           gATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCCTGGACGCCGAtag                                                                                                                                         >  7:240671/1‑133 (MQ=255)
                            aTCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCAt                                                                                                                          <  7:129697/149‑1 (MQ=255)
                               aaaGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCCTGGACGCCGAtagc                                                                                                                                        <  1:197145/132‑4 (MQ=255)
                               aaaGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCCTGGACGCCGAtagc                                                                                                                                        >  2:197145/1‑129 (MQ=255)
                                aaGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGt                                                                                                                      >  1:247325/1‑149 (MQ=255)
                                  gCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTgg                                                                                                                    >  3:212543/1‑149 (MQ=255)
                                     gAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCg                                                                                                                  >  3:126615/1‑148 (MQ=255)
                                     gAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGt                                                                                                                 <  5:15981/149‑1 (MQ=255)
                                       ggATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTa                                                                                                                <  4:33287/148‑1 (MQ=255)
                                        gATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGc                                                                                                              <  4:260963/149‑1 (MQ=255)
                                                    cACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATcag                                                                                                    <  4:98544/147‑1 (MQ=255)
                                                    cACCCACACCTCCGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTCAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATcag                                                                                                    <  4:94756/147‑1 (MQ=255)
                                                       ccACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCaga                                                                                                                              <  1:7271/118‑1 (MQ=255)
                                                       ccACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCaga                                                                                                                              >  2:7271/1‑118 (MQ=255)
                                                       ccACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATcagc                                                                                                   <  4:247225/145‑1 (MQ=255)
                                                       ccACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATcagc                                                                                                   >  3:247225/1‑145 (MQ=255)
                                                       ccACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAgg                                                                                                >  6:7874/1‑148 (MQ=255)
                                                       ccACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGCTGGCGTAGCTGAAGATCAGCAgg                                                                                                >  2:20665/1‑148 (MQ=255)
                                                       ccACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATAGGTTGCGTAGCTGAAGATCAGCAgg                                                                                                >  4:169572/1‑148 (MQ=255)
                                                            ccTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAg                                                                                          >  1:81873/1‑149 (MQ=255)
                                                             cTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCgg                                                                                                                                       >  6:244088/1‑103 (MQ=255)
                                                             cTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAgg                                                                                         >  5:189061/1‑149 (MQ=255)
                                                                ggTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCACCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCa                                                                                       <  5:232293/148‑1 (MQ=255)
                                                                    cccAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCACCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATcgcg                                                                                  >  2:21854/1‑149 (MQ=255)
                                                                     ccAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGTCCGATGCCGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATagcg                                                                                  >  2:50742/1‑148 (MQ=255)
                                                                     ccAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTa                                                                                              >  1:253827/1‑136 (MQ=255)
                                                                     ccAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTa                                                                                              <  2:253827/136‑1 (MQ=255)
                                                                       aGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAg                                                                               <  4:212543/149‑1 (MQ=255)
                                                                         cTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGt                                                                             <  2:24964/149‑1 (MQ=255)
                                                                                tCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTATGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGtt                                                                      >  1:202182/1‑149 (MQ=255)
                                                                                   tGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGt                                                                    >  6:4929/1‑148 (MQ=255)
                                                                                    gCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGa                                                                                                           <  3:48658/108‑1 (MQ=255)
                                                                                    gCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGa                                                                                                           >  4:48658/1‑108 (MQ=255)
                                                                                     cAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTaa                                                                  <  4:126615/148‑1 (MQ=255)
                                                                                       aaaTAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATc                                                                <  6:189061/149‑1 (MQ=255)
                                                                                       aaaTAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATc                                                                <  6:151740/149‑1 (MQ=255)
                                                                                          taatTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAg                                                                                          >  4:215506/1‑120 (MQ=255)
                                                                                          taatTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAg                                                                                          <  3:215506/120‑1 (MQ=255)
                                                                                              atTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCaga                                                                                                                              >  7:189209/1‑82 (MQ=255)
                                                                                              atTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCaga                                                                                                                              <  8:189209/82‑1 (MQ=255)
                                                                                              atTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGAcc                                                            <  2:218787/148‑1 (MQ=255)
                                                                                                               gCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCgtg                                          <  4:216417/149‑1 (MQ=255)
                                                                                                                cGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCgtgg                                         >  2:110743/1‑149 (MQ=255)
                                                                                                                   gAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGc                                      >  6:147045/1‑149 (MQ=255)
                                                                                                                            gTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGa                             >  3:4977/1‑149 (MQ=255)
                                                                                                                             ttACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGGGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGAt                            <  3:169572/149‑1 (MQ=255)
                                                                                                                                     ttGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCACGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCgg                    <  2:247325/149‑1 (MQ=255)
                                                                                                                                      tGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGTATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCgg                    <  2:81873/148‑1 (MQ=255)
                                                                                                                                      tGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCgg                    <  5:147045/148‑1 (MQ=255)
                                                                                                                                        ccggcgGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGTAGATCAGCAGGTTGAAGGGATCGCGCAGGTGGCCGTTGTAATCGTCGCTGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCa                 <  1:50742/146‑1 (MQ=255)
                                                                                                                                        cccGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCa                 >  5:215577/1‑149 (MQ=255)
                                                                                                                                         ccGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCgg                                    >  2:226182/1‑129 (MQ=255)
                                                                                                                                         ccGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCgg                                    <  1:226182/129‑1 (MQ=255)
                                                                                                                                         ccGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCAt                >  1:267342/1‑149 (MQ=255)
                                                                                                                                            cgGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCtt             <  5:7874/149‑1 (MQ=255)
                                                                                                                                             tgCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATATGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCGTTTGTAATCTACCCAGTCCTCCTGGACCGTGGTTCGGTAGGTGATCAGAGCGGGCATCTTg            <  1:20665/148‑1 (MQ=255)
                                                                                                                                                tATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTGGtc         >  3:133549/1‑149 (MQ=255)
                                                                                                                                                   cGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTg            <  6:205705/143‑1 (MQ=255)
                                                                                                                                                   cGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTg            >  5:205705/1‑143 (MQ=255)
                                                                                                                                                   cGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTGGTCTcc      >  1:31587/1‑149 (MQ=255)
                                                                                                                                                    ggCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTGGTCTcc      <  5:4929/148‑1 (MQ=255)
                                                                                                                                                    ggCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTGGTCTcc      <  1:110743/148‑1 (MQ=255)
                                                                                                                                                       gTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTGGTCTCCTTaa  <  2:202182/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GGCTGTGATACACGTGAAGACGCTTGCGATCAAAGCCGAGGATCTGTGTCTGCACCCACACCTCGGTGCCCAGCTTCACCTCGTGCAAATAATTGATATGCGCCTCCAGCGTGAACAGCGAGTTACCGCTTTGCCCGCGGCTATCGGCGTCCAGGCCGATGCGGTCCATCAGAGCATCGGTGGCGTAGCTGAAGATCAGCAGGTAGAAGGCATCGCGCAGGTGGCCGTTGTAATCGACCCAGTCCTCCTGGACCGTGGTGCGGTAGGTGATCAGAGCGGGCATCTTGGTCTCCTTAA  >  NC_002947/362224‑362520

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: