Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 591,079 G→A 6.6% R109H (CGC→CAC)  PP_0501 → NAD‑dependent epimerase/dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947591,0790GA 6.6% 217.3 / 4.1 75R109H (CGC→CAC) PP_0501NAD‑dependent epimerase/dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base A (3/2);  total (38/37)
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.34e-01

GATGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAG  >  NC_002947/590931‑591218
                                                                                                                                                    |                                                                                                                                           
gATGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATgcg                                                                                                                                             >  8:31999/1‑149 (MQ=255)
gATGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATgcg                                                                                                                                             >  2:336346/1‑149 (MQ=255)
   gCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGt                                                                                                                                          >  3:106376/1‑149 (MQ=255)
     aGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTg                                                                                                                                         >  4:21502/1‑148 (MQ=255)
            gTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGt                                                                                                                                 >  7:230446/1‑149 (MQ=255)
            gTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGt                                                                                                                                 >  4:225163/1‑149 (MQ=255)
            gTCACTCAGGCAGCGGCTGGCTGCCGGGCTGTGGTGCATCTGGCGGCGGCAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGGTTGCGAAGCCATGCACGCGCATGGctt                                                                                                                                 <  4:35114/149‑3 (MQ=255)
             tCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGTGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTg                                                                                                                                <  2:160919/149‑1 (MQ=255)
               tcTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCAt                                                                                                                                      <  6:242229/140‑1 (MQ=255)
               aCTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCACGCATCGGTCGAAGCCCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATgagatcgg                                                                                                                              >  5:242229/1‑142 (MQ=255)
               aCTCAGGCTGCGGCTGGCTGCCGGGCTGTGGTGCATCTGGCGGCGTTAGCGTCGTTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTgc                                                                                                                               <  8:350435/148‑1 (MQ=255)
                cTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTACATGGCGTgcgc                                                                                                                             >  4:203365/1‑149 (MQ=255)
                 tCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGcc                                                                                                                            >  6:70422/1‑149 (MQ=255)
                      cTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTg                                                                                                                       >  7:272627/1‑149 (MQ=255)
                         cGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGctc                                                                                                                    <  8:209530/149‑1 (MQ=255)
                             tggctgCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCg                                                                                                                >  7:23719/1‑149 (MQ=255)
                               gctgCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGc                                                                                                               <  8:23719/148‑1 (MQ=255)
                                 tgCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGCCCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGCAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTc                                                                                                            >  4:43251/1‑149 (MQ=255)
                                 tgCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTAACCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTc                                                                                                            >  6:227349/1‑149 (MQ=255)
                                 tgCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTAAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTc                                                                                                            >  2:339958/1‑149 (MQ=255)
                                     gggCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagc                                                                                                        <  5:64203/147‑1 (MQ=255)
                                      tGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagcg                                                                                                       <  6:159301/149‑1 (MQ=255)
                                           tggtgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagcgcagcg                                                                                                  >  6:215327/1‑149 (MQ=255)
                                           tggtgCATCTGGCGGCGGGAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCATACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagcgcagcg                                                                                                  <  5:324072/149‑1 (MQ=255)
                                             gtgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgt                                                                                                >  6:177133/1‑149 (MQ=255)
                                             gtgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgt                                                                                                >  8:144180/1‑149 (MQ=255)
                                             gtgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgt                                                                                                >  8:179614/1‑149 (MQ=255)
                                             gtgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGCTGCTCTTCGCCTCCAGCGCAGCGgt                                                                                                >  4:196242/1‑149 (MQ=255)
                                               gCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgtg                                                                                               >  4:317127/1‑148 (MQ=255)
                                                    tGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACgg                                                                                          >  5:333541/1‑148 (MQ=255)
                                                      gcggcgGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTcc                                                                                                           <  2:365225/129‑1 (MQ=255)
                                                      gcggcgGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTACATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACg                                                                                           <  3:203365/145‑1 (MQ=255)
                                                                cGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATGGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGaa                                                                             <  3:81957/149‑1 (MQ=255)
                                                                cGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGaa                                                                             <  1:336346/149‑1 (MQ=255)
                                                                 gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAg                                                                            >  5:4697/1‑149 (MQ=255)
                                                                 gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAg                                                                            >  6:260108/1‑149 (MQ=255)
                                                                 gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCACCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAg                                                                            >  7:140760/1‑149 (MQ=255)
                                                                    ggTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAggg                                                                          <  5:260108/148‑1 (MQ=255)
                                                                         aGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATc                                                                     <  8:272627/148‑1 (MQ=255)
                                                                           gCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGca                                                                                        <  7:113086/127‑1 (MQ=255)
                                                                           gCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGca                                                                                        >  8:113086/1‑127 (MQ=255)
                                                                             aTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGAttt                                                                >  2:359632/1‑149 (MQ=255)
                                                                              tCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTc                                                               <  3:196242/149‑1 (MQ=255)
                                                                                ggTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTcc                                                              >  7:152373/1‑148 (MQ=255)
                                                                                ggTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTcc                                                              <  8:230446/148‑1 (MQ=255)
                                                                                 gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCtt                                                                                                                  <  5:103582/95‑1 (MQ=255)
                                                                                 gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCtt                                                                                                                  >  6:103582/1‑95 (MQ=255)
                                                                                 gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGa                                                            <  7:31999/149‑1 (MQ=255)
                                                                                  tCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                           <  6:4697/149‑1 (MQ=255)
                                                                                  tCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGGGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGGGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                           <  8:141047/149‑1 (MQ=255)
                                                                                       gACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGacac                                                      <  7:271392/149‑1 (MQ=255)
                                                                                              tGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAg                                               <  5:227349/149‑1 (MQ=255)
                                                                                                  tACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGgcgc                                           <  5:215327/149‑1 (MQ=255)
                                                                                                    cccACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGcc                                          <  4:106376/148‑1 (MQ=255)
                                                                                                     ccaccaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGt                                        >  8:120271/1‑149 (MQ=255)
                                                                                                     ccaccaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGt                                        <  5:177133/149‑1 (MQ=255)
                                                                                                       accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGc                                                    <  6:138077/135‑1 (MQ=255)
                                                                                                       accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGc                                                    >  5:138077/1‑135 (MQ=255)
                                                                                                         caGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGAc                                    <  7:144180/149‑1 (MQ=255)
                                                                                                         caGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGAc                                    <  3:317127/149‑1 (MQ=255)
                                                                                                                  ttCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGc                           >  1:126563/1‑149 (MQ=255)
                                                                                                                    gaTCGGTACCCTCAACGTTTGCGATGCCATGCGCGTGCTTGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAgtg                                                                          <  6:325414/99‑3 (MQ=255)
                                                                                                                    cATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTCCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAggg                                                                          >  5:325414/1‑100 (MQ=255)
                                                                                                                      tCGGTACCCTCAACGTTTGCGAAGCAATGAGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                       <  5:70422/149‑1 (MQ=255)
                                                                                                                         gTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAgg                                                                           <  7:343751/94‑1 (MQ=255)
                                                                                                                         gTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAgg                                                                           >  8:343751/1‑94 (MQ=255)
                                                                                                                           aCCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGAc                    >  5:94966/1‑147 (MQ=255)
                                                                                                                               tCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTg              >  4:60234/1‑149 (MQ=255)
                                                                                                                                 aaCGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTg              <  6:333541/147‑1 (MQ=255)
                                                                                                                                 aaCGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTg              <  3:21502/147‑1 (MQ=255)
                                                                                                                                   cGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCa          <  7:120271/149‑1 (MQ=255)
                                                                                                                                         gcacgCCATGCACGCGCATGGCTTGCGCCGTGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCCATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAg     >  2:242226/5‑148 (MQ=255)
                                                                                                                                         gcacgCCATGCACGCGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                           <  2:109483/90‑1 (MQ=255)
                                                                                                                                         gcacgCCATGCACGCGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                           >  1:109483/5‑94 (MQ=255)
                                                                                                                                          gAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAg     >  1:342254/1‑147 (MQ=255)
                                                                                                                                          gAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCa   >  8:59515/1‑149 (MQ=255)
                                                                                                                                           acgCCATGCACGCGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAg  >  3:13278/3‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GATGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAG  >  NC_002947/590931‑591218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑