Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,681 A→T 7.0% intergenic (‑29/+28) PP_5073 ← / ← hemE hypothetical protein/uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,790,6810AT 7.0% 350.8 / 13.4 114intergenic (‑29/+28)PP_5073/hemEhypothetical protein/uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (43/63);  new base T (5/3);  total (48/66)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 116 total aligned reads displayed.

CCGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTGGCCA  >  NC_002947/5790532‑5790829
                                                                                                                                                     |                                                                                                                                                    
ccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGaa                                                                                                                                                      <  1:1102590/150‑1 (MQ=255)
    aaaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCg                                                                                                                                                  >  1:2097254/1‑150 (MQ=255)
     aaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                                 <  2:1698138/150‑1 (MQ=255)
     aaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                                 <  1:1401093/150‑1 (MQ=255)
     aaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                                 <  1:405531/150‑1 (MQ=255)
     aaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                                 <  2:735977/150‑1 (MQ=255)
       tGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATACTCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGt                                                                                                                                               <  1:2110847/150‑1 (MQ=255)
           ccAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTAc                                                                                                                                           <  1:2185235/150‑1 (MQ=255)
            cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACt                                                                                                                                          >  2:2035219/1‑150 (MQ=255)
            cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACt                                                                                                                                          >  2:1846883/1‑150 (MQ=255)
                 cTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACt                                                                                                                                          >  2:769895/1‑145 (MQ=255)
                 cTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACt                                                                                                                                          <  1:769895/145‑1 (MQ=255)
                 cTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCAGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgca                                                                                                                                     <  1:342109/150‑1 (MQ=255)
                 cTGAACAAGGGGGCGAGGAAGAGGTCGAGAGAACGGATTTTCATGATAATCGGGTCCTGTTTCGGGTTCGACGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGAGTTTACTgtgca                                                                                                                                     >  2:1332387/1‑150 (MQ=255)
                   gAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagt                                                                                                                                   <  1:2359242/150‑1 (MQ=255)
                   gAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTgtgcagt                                                                                                                                   <  2:2405142/150‑1 (MQ=255)
                    aaCAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGCTTACTgtgcagtg                                                                                                                                  >  1:2150869/1‑150 (MQ=255)
                         gggggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTTCGATGCCGGGtttt                                                                                                                                            <  1:2024950/135‑2 (MQ=255)
                         gggggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCa                                                                                                                             <  2:603185/150‑1 (MQ=255)
                           gggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                                 >  2:581758/1‑128 (MQ=255)
                           gggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCgg                                                                                                                                                 <  1:581758/128‑1 (MQ=255)
                           gggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCaat                                                                                                                           <  2:1971729/150‑2 (MQ=255)
                            ggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGt                                                                                                                          <  2:1641853/150‑1 (MQ=255)
                            ggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCAGGTTTACTGTGCAGTGCATCAAGt                                                                                                                          <  1:248026/150‑1 (MQ=255)
                               gAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCAc                                                                                                                       >  1:1429297/1‑150 (MQ=255)
                                  gATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtgg                                                                                                                    >  2:1488365/1‑150 (MQ=255)
                                    tGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtg                                                                                                                  <  2:801245/150‑1 (MQ=255)
                                           gagaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgc                                                                                                           >  1:2311001/1‑150 (MQ=255)
                                           gagaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgc                                                                                                           <  2:1367042/150‑1 (MQ=255)
                                           gagaAAACGGATTTTGATGTTAATCGGGTCCTGGTTCGGTTTCGATGGGACAAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgc                                                                                                           <  2:267783/150‑1 (MQ=255)
                                            agaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgcg                                                                                                          <  2:1429297/150‑1 (MQ=255)
                                            agaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgcg                                                                                                          <  2:2188446/150‑1 (MQ=255)
                                             gaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTgcgc                                                                                                         <  2:429916/150‑1 (MQ=255)
                                                 aCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCAt                                                                                                                               <  1:2186973/124‑1 (MQ=255)
                                                 aCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCAt                                                                                                                               >  2:2186973/1‑124 (MQ=255)
                                                 aCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGAc                                                                                                     >  2:189060/1‑150 (MQ=255)
                                                      ttttCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTc                                                                                                <  2:2150869/150‑1 (MQ=255)
                                                         tCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTg                                                                                             <  2:1371604/150‑1 (MQ=255)
                                                                  tCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAATTTATCGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGtt                                                                                    <  1:983038/150‑1 (MQ=255)
                                                                  tCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCCAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGtt                                                                                    <  1:189060/150‑1 (MQ=255)
                                                                  tCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGtt                                                                                    <  1:719554/150‑1 (MQ=255)
                                                                    gggTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTtga                                                                                  <  1:1453209/150‑1 (MQ=255)
                                                                      gTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACCTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTtgatg                                                                                >  2:1694386/1‑150 (MQ=255)
                                                                            tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgtg                                                                                                                                       <  2:1905347/89‑1 (MQ=255)
                                                                            tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACtgtg                                                                                                                                       >  1:1905347/1‑89 (MQ=255)
                                                                            tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCa                                                                          >  2:1058814/1‑150 (MQ=255)
                                                                               cggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCAccc                                                                       >  2:32308/1‑150 (MQ=255)
                                                                                 gtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGc                                                                     >  2:1723511/1‑150 (MQ=255)
                                                                                 gtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGc                                                                     >  2:2192770/1‑150 (MQ=255)
                                                                                  tttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGTGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCa                                                                    <  2:1164173/150‑1 (MQ=255)
                                                                                  tttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCa                                                                    >  1:949879/1‑150 (MQ=255)
                                                                                   ttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCAt                                                                   <  1:1058814/150‑1 (MQ=255)
                                                                                      gATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGtt                                                                >  2:1460521/1‑150 (MQ=255)
                                                                                            aCCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGt                                                          <  1:1460521/150‑1 (MQ=255)
                                                                                              cAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCa                                                        >  2:1109339/1‑150 (MQ=255)
                                                                                                 aTTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggt                                                                                                                   <  1:1573749/88‑1 (MQ=255)
                                                                                                 aTTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggt                                                                                                                   >  2:1573749/1‑88 (MQ=255)
                                                                                                       gacgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAATTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCa                                                                          >  2:1662339/1‑123 (MQ=255)
                                                                                                       gacgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCa                                                                          <  1:1662339/123‑1 (MQ=255)
                                                                                                         cgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGt                                             <  1:32308/150‑1 (MQ=255)
                                                                                                          gaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGt                                                          <  1:1442270/136‑1 (MQ=255)
                                                                                                            tCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTa                                                                                                                                            <  2:488875/52‑1 (MQ=255)
                                                                                                            tCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTa                                                                                                                                            >  1:488875/1‑52 (MQ=255)
                                                                                                               cTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGcc                                       >  2:838163/1‑150 (MQ=255)
                                                                                                                tGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGt                                                                                              >  2:2389857/1‑94 (MQ=255)
                                                                                                                tGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGt                                                                                              >  2:2389848/1‑94 (MQ=255)
                                                                                                                tGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGt                                                                                              <  1:2389848/94‑1 (MQ=255)
                                                                                                                tGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCg                                      >  1:1744912/1‑150 (MQ=255)
                                                                                                                cgAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGt                                                                                              <  1:2390186/93‑1 (MQ=255)
                                                                                                                cgAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGt                                                                                              >  2:2390186/2‑94 (MQ=255)
                                                                                                                 gAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGt                                     <  1:209423/150‑1 (MQ=255)
                                                                                                                  aaTTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTg                                    <  1:1469333/150‑1 (MQ=255)
                                                                                                                  aaTTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTg                                    <  1:1109339/150‑1 (MQ=255)
                                                                                                                   aTTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTgg                                   <  1:375258/150‑1 (MQ=255)
                                                                                                                    ttCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGt                                             >  2:2161024/1‑139 (MQ=255)
                                                                                                                       aaaCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGtt                                                                                                 <  2:1627745/84‑1 (MQ=255)
                                                                                                                         aCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAgg                             <  2:1859122/150‑1 (MQ=255)
                                                                                                                         aCTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAgg                             <  2:452432/150‑1 (MQ=255)
                                                                                                                              aggtgcAAAACGCCCGGCATCGAAGCCGGGGTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGGATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGaa                        <  2:18368/145‑1 (MQ=255)
                                                                                                                                 tACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAacac                     <  1:1694386/150‑1 (MQ=255)
                                                                                                                                  aCAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAAcaca                    <  1:1798848/150‑1 (MQ=255)
                                                                                                                                   cAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGAACGGCGTTGATGAATCCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACAt                   <  1:1488365/150‑1 (MQ=255)
                                                                                                                                    agtaaaCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGcc                                       >  1:632219/4‑129 (MQ=255)
                                                                                                                                    aaaaaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATg                  >  2:577916/1‑150 (MQ=255)
                                                                                                                                       aaaCCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGcc               >  2:1458097/1‑150 (MQ=255)
                                                                                                                                        aaCCCGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCa              >  1:2241762/1‑150 (MQ=255)
                                                                                                                                          cccGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGcc                                       >  1:809006/1‑123 (MQ=255)
                                                                                                                                          cccGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGcc                                       <  2:809006/123‑1 (MQ=255)
                                                                                                                                           ccGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTgg                                   >  1:2266748/1‑126 (MQ=255)
                                                                                                                                           ccGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTgg                                   <  2:2266748/126‑1 (MQ=255)
                                                                                                                                           ccGGCTTCGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCATGCTCTACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAggg           >  1:146389/1‑148 (MQ=255)
                                                                                                                                            cGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCAc                                                                         <  1:1444281/87‑1 (MQ=255)
                                                                                                                                            cGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCAc                                                                         >  2:1444281/1‑87 (MQ=255)
                                                                                                                                            cGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGc          >  2:860516/1‑150 (MQ=255)
                                                                                                                                            cGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGc          >  2:1048296/1‑150 (MQ=255)
                                                                                                                                            cGGCATCGAAGCCGGGTTTACAGTACAGTGCATCATGTCACTGGTGGTACGGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCGCCAGCATGTTCCGGGTAAACTGGCGGTGTGATGCCGGCGCCGAGGTGGAGGGCAGGGCCAGTGc          <  2:1207571/150‑1 (MQ=255)
                                                                                                                                             ggCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGcc         >  1:854043/1‑150 (MQ=255)
                                                                                                                                                aTCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGaa                                                                              <  2:1445326/78‑1 (MQ=255)
                                                                                                                                                 tCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTgg     <  1:838163/150‑1 (MQ=255)
                                                                                                                                                    aaGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTGGCCa  <  1:2035219/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCATCGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTGGCCA  >  NC_002947/5790532‑5790829

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