Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 83,541 T→C 5.7% L94P (CTG→CCG)  hemF → coproporphyrinogen‑III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294783,5410TC 5.7% 251.7 / 4.0 87L94P (CTG→CCG) hemFcoproporphyrinogen‑III oxidase
Reads supporting (aligned to +/- strand):  ref base T (43/39);  new base C (2/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GCCGCACGCGGGTTATCGGTGAGGGCAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCAACGAAGAAGACTGCATT  >  NC_002947/83394‑83689
                                                                                                                                                   |                                                                                                                                                    
gCCGCACGCGGGTTATCGGTGAGGGCAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTg                                                                                                                                                     >  2:226641/1‑149 (MQ=255)
  cGCACGCGGGTTATCGGGGAGGGCAACGTGATCGAAAACCGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGt                                                                                                                                                   <  8:161051/149‑1 (MQ=255)
              aTCGGTGAGGGCAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGca                                                                                                                                       >  2:291096/1‑149 (MQ=255)
              aTCGGTGAGGGCAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTATCACCAGCGACACAcccaa                                                                                                                                        <  1:235322/148‑3 (MQ=255)
               tCGGTGAGGGCAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTATCACCAGCGACACAcccaagg                                                                                                                                      <  3:264651/149‑5 (MQ=255)
                  gTGAGGTCTACGTGATCGAAAAAGGCGGGGTCAACTTCTACCATGTGTTTGGTAGCGGCATGCCGCCGTCGGACAGTGCCCACCGTCCCTAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAAc                                                                                                                                   <  4:362313/149‑1 (MQ=255)
                  gTGAGGGCAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAAc                                                                                                                                   >  5:48769/1‑149 (MQ=255)
                         cAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTATCACCAGCGACACAcccaaggc                                                                                                                                     <  8:342131/140‑6 (MQ=255)
                         cAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTATCACCAGCGACACAcccaaggc                                                                                                                                     >  7:342131/1‑135 (MQ=255)
                          aaCGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCAtgt                                                                                                                           <  7:281295/149‑1 (MQ=255)
                             gTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGcc                                                                                                                        <  6:161555/149‑1 (MQ=255)
                             gTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGcc                                                                                                                        <  4:243592/149‑1 (MQ=255)
                               gATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCa                                                                                                                      >  8:36831/1‑149 (MQ=255)
                                  cGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACtt                                                                                                                   >  1:418778/1‑149 (MQ=255)
                                   gAAAAAGGCGGGGGCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTc                                                                                                                  <  7:336289/149‑1 (MQ=255)
                                    aaaaaGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTcc                                                                                                                 <  1:374926/149‑1 (MQ=255)
                                    aaaaaGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTcc                                                                                                                 >  2:162031/1‑149 (MQ=255)
                                    aaaaaGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTcc                                                                                                                 >  2:203856/1‑149 (MQ=255)
                                     aaaaGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTccc                                                                                                                <  8:44152/149‑1 (MQ=255)
                                     aaaaGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTcc                                                                                                                 >  8:27104/1‑148 (MQ=255)
                                            ggggTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCaa                                                                                                         >  5:403278/1‑149 (MQ=255)
                                              ggTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACg                                                                                                       >  5:293546/1‑149 (MQ=255)
                                                tCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGtt                                                                                                     <  1:162031/149‑1 (MQ=255)
                                                tCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGtt                                                                                                     <  8:142160/149‑1 (MQ=255)
                                                  aaCTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCg                                                                                                   >  1:417622/1‑149 (MQ=255)
                                                   aCTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGc                                                                                                  <  5:189819/149‑1 (MQ=255)
                                                     ttCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGctt                                                                                                >  6:122821/1‑149 (MQ=255)
                                                          ccATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGcttctt                                                                                             <  8:327357/147‑1 (MQ=255)
                                                          ccATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCttcttc                                                                                            <  5:122821/148‑1 (MQ=255)
                                                            aTGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATc                                                                                         <  8:93242/149‑1 (MQ=255)
                                                            aTGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATc                                                                                         >  1:14947/1‑149 (MQ=255)
                                                               tgtTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGcc                                                                                      <  1:85852/149‑1 (MQ=255)
                                                               tgtTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGcc                                                                                      >  5:421511/1‑149 (MQ=255)
                                                                gtTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGcc                                                                                      >  8:255763/1‑148 (MQ=255)
                                                                    ggTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGaaa                                                                                  <  2:418778/148‑1 (MQ=255)
                                                                       cacggcctgccagttcgggacggtgggcactggccgacggcggcaggccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGcttc                                                                                               >  8:259370‑M2/51‑132 (MQ=255)
                                                                       cacggcctgccagttcgggacggtgggcactggccgacggcggcaggccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGcttc                                                                                               <  7:259370‑M2/82‑1 (MQ=255)
                                                                        gCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCGCAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAaaggaa                                                                             <  7:114782/149‑1 (MQ=255)
                                                                        gCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAaagga                                                                              <  1:291096/148‑1 (MQ=255)
                                                                        acggcctgccagttcgggacggtgggccctggccgaccgcggcaggccgTGGCGTCCAAGACTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCAACTTCCCCAGCGAACGTTCGCTTCTTCAGCGCCGAAaagg                                                                               >  7:300855‑M2/50‑147 (MQ=255)
                                                                         cggcctgccagttcgggacggtgggcactggccgacggcggcaggccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGcttc                                                                                               <  2:394193‑M2/82‑1 (MQ=255)
                                                                         cggcctgccagttcgggacggtgggcactggccgacggcggcaggccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGcttc                                                                                               >  1:394193‑M2/49‑130 (MQ=255)
                                                                         cggcctgccagttcgggacggtgggcactggccgacggcggcaggccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGcttc                                                                                               >  1:389001‑M2/49‑130 (MQ=255)
                                                                         cggcctgccagttcgggacggtgggcactggccgacggcggcaggccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGcttc                                                                                               <  2:389001‑M2/82‑1 (MQ=255)
                                                                         CGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCGGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAaaggaag                                                                            <  6:84391‑M1/149‑1 (MQ=255)
                                                                             CTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCCCACAACCCGCATGTGCACACTTCCCACGCCAACGTTCGCTTCTTCAACGCCGAAAAGGAAGGTg                                                                         <  4:211036‑M1/148‑1 (MQ=255)
                                                                              TGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTg                                                                         <  6:385515‑M1/147‑1 (MQ=255)
                                                                                CCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAag                                                                      >  1:226641‑M1/1‑148 (MQ=255)
                                                                                 CGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAagag                                                                    <  1:51835‑M1/149‑1 (MQ=255)
                                                                                  gccgTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAgg                                                                   >  3:13121/1‑149 (MQ=255)
                                                                                  gccgTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAgg                                                                   >  4:161442/1‑149 (MQ=255)
                                                                                  gccgTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCCGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAgg                                                                   >  2:265424/1‑149 (MQ=255)
                                                                                   ccgTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGACGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGc                                                                  >  7:285742/1‑149 (MQ=255)
                                                                                        ggCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCCTGTGCCCACTTCCCACGCCAACGTTc                                                                                                    <  1:232173/110‑1 (MQ=255)
                                                                                        ggCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTc                                                                                                    >  2:232173/1‑110 (MQ=255)
                                                                                        ggCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGACTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTAt                                                             >  1:64719/1‑149 (MQ=255)
                                                                                              tGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATggtggt                                                       >  8:399191/1‑149 (MQ=255)
                                                                                               gCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGcc                                                                 >  2:50558/1‑138 (MQ=255)
                                                                                               gCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGcc                                                                 <  1:50558/138‑1 (MQ=255)
                                                                                                cccACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGc                                                                                       <  8:72797/115‑1 (MQ=255)
                                                                                                cccACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGc                                                                                       >  7:72797/1‑115 (MQ=255)
                                                                                                    ccGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTggcg                                                 <  6:403278/149‑1 (MQ=255)
                                                                                                     cGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTggcgg                                                <  2:51835/149‑1 (MQ=255)
                                                                                                        cccGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTggcggcgg                                             >  1:52059/1‑149 (MQ=255)
                                                                                                         ccGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTggcggcggc                                            >  2:33320/1‑149 (MQ=255)
                                                                                                          cGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCt                                           >  7:44854/1‑149 (MQ=255)
                                                                                                             aCTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTg                                        <  2:14947/149‑1 (MQ=255)
                                                                                                              cTGGCAGGGCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGa                                       <  8:285742/149‑1 (MQ=255)
                                                                                                                 gcaggccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGAcc                                     <  2:417622‑M2/140‑1 (MQ=255)
                                                                                                                   aggccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGa                                  <  3:328735‑M2/143‑1 (MQ=255)
                                                                                                                     gccgTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACg                                >  3:393722‑M2/5‑149 (MQ=255)
                                                                                                                         tGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCAt                            <  5:385515/149‑1 (MQ=255)
                                                                                                                             ttCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTa                        >  7:406330/1‑149 (MQ=255)
                                                                                                                             ttCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTa                        >  3:133163/1‑149 (MQ=255)
                                                                                                                               cGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCACCGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATg                      >  6:251071/1‑149 (MQ=255)
                                                                                                                                gAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATgg                     >  1:410354/1‑149 (MQ=255)
                                                                                                                                 aaGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATgg                     >  3:21965/1‑148 (MQ=255)
                                                                                                                                   gCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGcc                              <  8:146711/137‑1 (MQ=255)
                                                                                                                                   gCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGcc                              >  7:146711/1‑137 (MQ=255)
                                                                                                                                   gCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCaa                  <  7:399191/149‑1 (MQ=255)
                                                                                                                                         ggTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCAACgaag             >  3:95938/1‑148 (MQ=255)
                                                                                                                                             tgtCGCTGGTGCTACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAATGTGAAGAGGCCGTATGGTGGTTTGGCGGCGTCTTTGACCTTACGCCATACTATGGCAACGAagaaga         <  6:48769/148‑1 (MQ=255)
                                                                                                                                               tCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCAACGAAGAAGACt       <  6:421511/148‑1 (MQ=255)
                                                                                                                                               tCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCAACGAAGAAGACTg      >  5:64283/1‑149 (MQ=255)
                                                                                                                                               tCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTAGGGAAACGGAGAAGACTg      >  4:17157/1‑149 (MQ=255)
                                                                                                                                                 gCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCAACGAAGAAGACTGCa    >  7:407042/1‑149 (MQ=255)
                                                                                                                                                 gCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCAACGAAGAAGACTGCa    >  2:23558/1‑149 (MQ=255)
                                                                                                                                                   tGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCAACGAAGAAGACTGCAtt  <  2:410354/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCGCACGCGGGTTATCGGTGAGGGCAACGTGATCGAAAAAGGCGGGGTCAACTTCTCCCATGTGTTTGGTAGCGGCCTGCCGCCGTCGGCCAGTGCCCACCGTCCCGAACTGGCAGGCCGTGGCTTCGAAGCCTTGGGTGTGTCGCTGGTGATACACCCGCACAACCCGCATGTGCCCACTTCCCACGCCAACGTTCGCTTCTTCATCGCCGAAAAGGAAGGTGAAGAGGCCGTATGGTGGTTTGGCGGCGGCTTTGACCTGACGCCATACTATGGCAACGAAGAAGACTGCATT  >  NC_002947/83394‑83689

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: