Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,116,556 C→A 5.0% A255A (GCC→GCA)  lpd → dihydrolipoyl dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,116,5560CA 5.0% 257.8 / 2.6 79A255A (GCC→GCA) lpddihydrolipoyl dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (35/40);  new base A (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TGGACCGCATCTGCCCAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAGCACCACCGTACCT  >  NC_002947/6116408‑6116701
                                                                                                                                                    |                                                                                                                                                 
tGGACCGCATCTGCCCAGGCACTGACACGGAAACGGCCAAGACACTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGCCTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTgcc                                                                                                                                                   >  6:101846/1‑149 (MQ=255)
  gACCGCATCTGCCCAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTgccgc                                                                                                                                                 >  7:287605/1‑149 (MQ=255)
      gCATCTGCCAAGGCACTGACACGGCAACGGCCAAGACCATGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCggg                                                                                                                                             <  8:287605/149‑1 (MQ=255)
       cATCTGCCCAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCgggg                                                                                                                                            <  7:228493/149‑1 (MQ=255)
        aTCTGCCCAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCggggg                                                                                                                                           >  6:228853/1‑149 (MQ=255)
            gCCCAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTAcc                                                                                                                                       >  2:276717/1‑149 (MQ=255)
              ccAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCAGCCGGGGGTACCg                                                                                                                                      >  3:264007/1‑148 (MQ=255)
              ccAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCAGCCGGGGGTACCg                                                                                                                                      >  3:259462/1‑148 (MQ=255)
               cAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCa                                                                                                                                    <  8:23562/149‑1 (MQ=255)
                      tGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGc                                                                                                                             >  4:263454/1‑149 (MQ=255)
                        acacGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCt                                                                                                                            <  3:270660/148‑1 (MQ=255)
                        acacGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTg                                                                                                                           >  5:1948/1‑149 (MQ=255)
                        acacGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTg                                                                                                                           >  3:136919/1‑149 (MQ=255)
                        acacGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTg                                                                                                                           >  6:70107/1‑149 (MQ=255)
                           cGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCAGCCGGGGGTACCGCAGAATCGCTGCaa                                                                                                                        <  4:264007/149‑1 (MQ=255)
                           cGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCAGCCGGGGGTACCGCAGAATCGCTGCaa                                                                                                                        <  4:259462/149‑1 (MQ=255)
                             gAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGc                                                                                                                      <  5:228853/149‑1 (MQ=255)
                             gAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGc                                                                                                                      >  8:47997/1‑149 (MQ=255)
                               aaCGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCg                                                                                                                    <  7:146702/149‑1 (MQ=255)
                               aaCGGCCAAGACCCTGCAGAAGGCGCTAGCCAAACAAAGCATGGTGTTCAAACTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATAGCTGCAAGcc                                                                                                                     >  7:50490/1‑148 (MQ=255)
                                    ccAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTAc                                                                                                               >  3:300243/1‑149 (MQ=255)
                                     cAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGTCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACg                                                                                                              >  7:176285/1‑149 (MQ=255)
                                     cAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACg                                                                                                              >  5:245715/1‑149 (MQ=255)
                                      aaGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACgt                                                                                                             <  2:122134/149‑1 (MQ=255)
                                        gACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACgtgt                                                                                                           <  2:145135/149‑1 (MQ=255)
                                          cccTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTtg                                                                                                         <  4:135337/149‑1 (MQ=255)
                                            cTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTtggt                                                                                                       >  3:164482/1‑149 (MQ=255)
                                             tGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTtggtg                                                                                                      >  8:132591/1‑149 (MQ=255)
                                             tGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTtggtg                                                                                                      >  4:122374/1‑149 (MQ=255)
                                               cagcagGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTtggtg                                                                                                      <  4:25592/147‑1 (MQ=255)
                                               cAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTtggtg                                                                                                      <  7:47997/147‑1 (MQ=255)
                                                aGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGc                                                                                                    >  5:8732/1‑148 (MQ=255)
                                                   aGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATc                                                                                                <  1:91928/149‑1 (MQ=255)
                                                             caagcaagGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGcc                                                                                      >  3:80376/1‑149 (MQ=255)
                                                              aagcaagGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGc                                                                                                                                     >  8:25270/1‑101 (MQ=255)
                                                              aagcaagGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGc                                                                                                                                     <  7:25270/101‑1 (MQ=255)
                                                              aagcaagGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCg                                                                                               <  7:123853/139‑1 (MQ=255)
                                                              aagcaagGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCg                                                                                               >  8:123853/1‑139 (MQ=255)
                                                              aagcaagGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGccc                                                                                     <  3:288226/149‑1 (MQ=255)
                                                                          gtgtTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAg                                                                          <  5:8886/148‑1 (MQ=255)
                                                                           tgtTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAggg                                                                        <  4:136919/149‑1 (MQ=255)
                                                                            gtTCAAGCTCGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGc                                                                       >  3:117360/1‑149 (MQ=255)
                                                                               cAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTaa                                                                    <  8:198877/149‑1 (MQ=255)
                                                                               cAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTAGCGGCCGTCGCCCTTATACAAAAGGGCTa                                                                     >  5:301410/1‑148 (MQ=255)
                                                                                aaGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTTTGTATAAGGGcgacg                                                                   <  6:8732/149‑5 (MQ=255)
                                                                                aaGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTaaa                                                                   <  6:1948/149‑1 (MQ=255)
                                                                                aaGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTaa                                                                    <  5:70107/148‑1 (MQ=255)
                                                                                      ggCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTgg                                                              <  4:300243/148‑1 (MQ=255)
                                                                                      ggCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTg                                                               <  7:132591/147‑1 (MQ=255)
                                                                                        cAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGaaa                                                           <  3:263454/149‑1 (MQ=255)
                                                                                             aaGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGTCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTg                                                      <  8:176285/149‑1 (MQ=255)
                                                                                             aaGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTg                                                      <  4:80376/149‑1 (MQ=255)
                                                                                             aaGTGACTCACGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTg                                                      >  8:169132/1‑149 (MQ=255)
                                                                                              aGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAAGGCTAAACCTGGAAAGCGTgg                                                     <  6:246757/149‑1 (MQ=255)
                                                                                                tGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGc                                                   >  2:73306/1‑149 (MQ=255)
                                                                                                  aCTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCt                                                 >  5:52296/1‑149 (MQ=255)
                                                                                                   cTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTgg                                                     <  1:245711/144‑1 (MQ=255)
                                                                                                   cTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTgg                                                     >  2:245711/1‑144 (MQ=255)
                                                                                                   cTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTg                                                >  2:252874/1‑149 (MQ=255)
                                                                                                        gCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCt                                                                      <  6:289689/122‑1 (MQ=255)
                                                                                                        gCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCt                                                                      >  5:289689/1‑122 (MQ=255)
                                                                                                         cTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGcc                                          <  6:52296/149‑1 (MQ=255)
                                                                                                               cAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAg                                    <  6:123366/149‑1 (MQ=255)
                                                                                                                aGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCTCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCGTTATACAAAAGGGCTAAACGTGGAAAGCGTGGGCCTGGAAGCCGACAAgc                                   <  7:169132/149‑1 (MQ=255)
                                                                                                                aGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAgc                                   <  4:164482/149‑1 (MQ=255)
                                                                                                                  gTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAgcgc                                 <  7:103695/149‑1 (MQ=255)
                                                                                                                     ccGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCgg                               >  1:109706/1‑148 (MQ=255)
                                                                                                                     ccGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCgg                               <  6:301410/148‑1 (MQ=255)
                                                                                                                     ccGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGc                              <  6:245715/149‑1 (MQ=255)
                                                                                                                        aTGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTgg                                                     >  4:273568/1‑123 (MQ=255)
                                                                                                                        aTGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTgg                                                     <  3:273568/123‑1 (MQ=255)
                                                                                                                                 gCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCaacaac                  <  3:122374/149‑1 (MQ=255)
                                                                                                                                     gACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAg               >  7:80525/1‑148 (MQ=255)
                                                                                                                                        gCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAGcac            >  4:149084/1‑148 (MQ=255)
                                                                                                                                         cTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTTTACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAGcacca          <  1:276717/149‑1 (MQ=255)
                                                                                                                                         cTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCt                                                 <  5:288128/110‑1 (MQ=255)
                                                                                                                                         cTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCt                                                 >  6:288128/1‑110 (MQ=255)
                                                                                                                                            gAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAGCACCACCg       >  2:315516/1‑149 (MQ=255)
                                                                                                                                            gAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAGCACCACCg       >  6:56472/1‑149 (MQ=255)
                                                                                                                                               ccTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAGCACCACCGTAc    >  6:19812/1‑149 (MQ=255)
                                                                                                                                                 tGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAGCACCACCGTACCt  >  1:65865/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGGACCGCATCTGCCCAGGCACTGACACGGAAACGGCCAAGACCCTGCAGAAGGCGCTCGCCAAGCAAGGCATGGTGTTCAAGCTGGGCAGTAAAGTGACTCAGGCTACCGCAGGTGCCGATGGCGTCAGCCTGACGCTGGAGCCTGCCGCCGGGGGTACCGCAGAATCGCTGCAAGCCGATTACGTGTTGGTGGCTATCGGCCGTCGCCCTTATACAAAAGGGCTAAACCTGGAAAGCGTGGGCCTGGAAGCCGACAAGCGCGGCATGCTCAACAACGAGCACCACCGTACCT  >  NC_002947/6116408‑6116701

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: ‑