Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,948,640 T→G 6.2% D162A (GAC→GCC)  PP_2582 ← heme oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,948,6400TG 6.2% 241.7 / 7.6 90D162A (GAC→GCC) PP_2582heme oxygenase
Reads supporting (aligned to +/- strand):  ref base T (47/37);  new base G (2/4);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CTGTGCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCACTGCAAT  >  NC_002947/2948498‑2948787
                                                                                                                                              |                                                                                                                                                   
cTGTGCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCt                                                                                                                                              <  2:838586/150‑1 (MQ=255)
 tgtgCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTg                                                                                                                                             <  2:115283/150‑1 (MQ=255)
  gtgCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGGACAGCCTgg                                                                                                                                            <  2:934338/150‑1 (MQ=255)
       gcaCCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCt                                                                                                                                              <  1:370994/143‑1 (MQ=255)
       gcaCCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCt                                                                                                                                              >  2:370994/1‑143 (MQ=255)
             gCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTcg                                                                                                                                 >  2:608742/1‑150 (MQ=255)
             gCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTcg                                                                                                                                 >  1:1212167/1‑150 (MQ=255)
             gCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGGACAGCCTggtgctgg                                                                                                                                      >  1:1073048/1‑139 (MQ=255)
             gCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGGACAGCCTggtgctgg                                                                                                                                      <  2:1073048/145‑7 (MQ=255)
                    gTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAAt                                                                                                                          >  1:181744/1‑150 (MQ=255)
                       ccATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCa                                                                                                                       >  1:635010/1‑150 (MQ=255)
                        cATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAg                                                                                                                      >  1:212527/1‑150 (MQ=255)
                              tcTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGGCCGGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCg                                                                                                                <  1:423720/149‑1 (MQ=255)
                               cTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCgg                                                                                                               >  2:559918/1‑150 (MQ=255)
                                tCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGc                                                                                                              >  1:892053/1‑150 (MQ=255)
                                  ggTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCACTGCAGCGCTCggaga                                                                                                            >  2:423720/1‑147 (MQ=255)
                                   gTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCAt                                                                                                           >  2:39857/1‑150 (MQ=255)
                                    tCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATc                                                                                                          >  2:234276/1‑150 (MQ=255)
                                       aCTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCgg                                                                                                       <  2:498973/150‑1 (MQ=255)
                                       aCTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCgg                                                                                                       <  1:944983/150‑1 (MQ=255)
                                        ctgcctgcAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCggg                                                                                                      >  2:478599/1‑150 (MQ=255)
                                        ctgcctgcAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCggg                                                                                                      >  1:1211636/1‑150 (MQ=255)
                                             tgcAAATACCTCACGCAGATGCACCCGATAACTGTGGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTGCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCACTGCAGCGCTCGGGCAACCGGGTGAgg                                                                                                 >  2:57071/1‑150 (MQ=255)
                                             tgcAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAgg                                                                                                 >  1:230536/1‑150 (MQ=255)
                                             tgcAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAgg                                                                                                 >  1:795132/1‑150 (MQ=255)
                                                aaaTACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGc                                                                                              <  1:233104/150‑1 (MQ=255)
                                                   tACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAg                                                                                           <  1:39857/150‑1 (MQ=255)
                                                      cTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTg                                                                                        >  1:281719/1‑150 (MQ=255)
                                                       tCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTgg                                                                                       >  1:829880/1‑150 (MQ=255)
                                                       tCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGATGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTgg                                                                                       >  1:254956/1‑150 (MQ=255)
                                                          cGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCgg                                                                                    <  2:1212167/150‑1 (MQ=255)
                                                                   aCCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTg                                                                           >  1:1096496/1‑150 (MQ=255)
                                                                    cccGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTcgc                                                                                                                                >  1:994112/1‑96 (MQ=255)
                                                                    cccGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTcgc                                                                                                                                <  2:994112/96‑1 (MQ=255)
                                                                       gATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCAc                                                                       >  1:695825/1‑150 (MQ=255)
                                                                             tGTCGAACGCTGCACTCGCCCCCGCTCCCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGTTCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCa                                                                 <  2:529274/150‑1 (MQ=255)
                                                                                cGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGcc                                                              <  2:829880/150‑1 (MQ=255)
                                                                                 gAACGCTGTAATCGCCCCCGCTACCACCTCTGCTCCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGcc                                                              <  2:960100/149‑1 (MQ=255)
                                                                                 gAACGCTGTAATCGCCCCCGCTACCACCTCTGCTCCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGcc                                                              >  1:960100/1‑149 (MQ=255)
                                                                                 gAACGCTGCAATCGCAACCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGccc                                                             >  2:955339/1‑150 (MQ=255)
                                                                                  aaCGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCa                                                            >  1:1110862/1‑150 (MQ=255)
                                                                                  aaCGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCa                                                            >  1:1182044/1‑150 (MQ=255)
                                                                                   aCGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCaa                                                           <  2:230536/150‑1 (MQ=255)
                                                                                       tGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCgg                                                                                                       <  2:45919/102‑1 (MQ=255)
                                                                                       tGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCgg                                                                                                       >  1:45919/1‑102 (MQ=255)
                                                                                       tGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACg                                                       >  2:1131820/1‑150 (MQ=255)
                                                                                       tGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACg                                                       >  1:1000202/1‑150 (MQ=255)
                                                                                        gCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGt                                                      >  2:716694/1‑150 (MQ=255)
                                                                                        gCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGATGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGt                                                      <  2:254956/150‑1 (MQ=255)
                                                                                             cGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCgt                                                 <  2:1182044/150‑1 (MQ=255)
                                                                                               cccccGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCCAGACGTTGCgtgt                                               >  2:1174425/1‑150 (MQ=255)
                                                                                                 cccGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCGAGGTTGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGc                                             <  2:181744/150‑1 (MQ=255)
                                                                                                 cccGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGc                                             >  2:890202/1‑150 (MQ=255)
                                                                                                   cGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCTGCAGACGGCCATCCGGCCGATGTGCCAGGGGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCtt                                           <  1:478599/150‑1 (MQ=255)
                                                                                                      tACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTa                                        <  1:890202/150‑1 (MQ=255)
                                                                                                       accaccTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTTCTTGTag                                       <  1:234276/150‑1 (MQ=255)
                                                                                                           ccTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGaa                                   >  2:836988/1‑150 (MQ=255)
                                                                                                           ccTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGaa                                   >  2:292872/1‑150 (MQ=255)
                                                                                                               tGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTaga                                      <  1:180757/143‑1 (MQ=255)
                                                                                                               tGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTaga                                      >  2:180757/1‑143 (MQ=255)
                                                                                                                gCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAAcgcc                              <  2:57231/150‑1 (MQ=255)
                                                                                                                      tcctGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAAcgccgc                            <  1:530109/146‑1 (MQ=255)
                                                                                                                      tcctGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAAcgccgc                            >  2:530109/1‑146 (MQ=255)
                                                                                                                         tGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAgg                     <  1:1131820/150‑1 (MQ=255)
                                                                                                                           cTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGtt                   <  1:955339/150‑1 (MQ=255)
                                                                                                                            tCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCa                                                                        >  2:1261419/1‑96 (MQ=255)
                                                                                                                            tCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCa                                                                        <  1:1261419/96‑1 (MQ=255)
                                                                                                                             cccccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTagag                                     >  1:672364/1‑130 (MQ=255)
                                                                                                                             cccccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTagag                                     <  2:672364/130‑1 (MQ=255)
                                                                                                                              ccccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGa                >  2:129940/1‑150 (MQ=255)
                                                                                                                               cccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGaa               >  2:1010613/1‑150 (MQ=255)
                                                                                                                                ccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAAc              <  1:448255/150‑1 (MQ=255)
                                                                                                                                ccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAAc              <  1:608742/150‑1 (MQ=255)
                                                                                                                                ccAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAAc              >  2:448255/1‑150 (MQ=255)
                                                                                                                                 cAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAAcc             >  1:192930/1‑150 (MQ=255)
                                                                                                                                  gttgccAGGCTGGACAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAAc                                 <  1:476185/125‑1 (MQ=255)
                                                                                                                                  gttgccAGGCTGGACAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAAc                                 >  2:476185/5‑129 (MQ=255)
                                                                                                                                   gCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCtt           <  2:695825/150‑1 (MQ=255)
                                                                                                                                     aCCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCa         <  1:1174425/150‑1 (MQ=255)
                                                                                                                                     aCCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCa         <  2:635010/150‑1 (MQ=255)
                                                                                                                                      ccAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCAc        <  1:1010613/150‑1 (MQ=255)
                                                                                                                                       cAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGa                                                         >  2:1000981/1‑100 (MQ=255)
                                                                                                                                       cAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGa                                                         <  1:1000981/100‑1 (MQ=255)
                                                                                                                                       cAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGccc                          <  2:1030820/131‑1 (MQ=255)
                                                                                                                                       cAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGccc                          >  1:1030820/1‑131 (MQ=255)
                                                                                                                                        aGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGaa               >  2:667297/1‑141 (MQ=255)
                                                                                                                                        aGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGaa               <  1:667297/141‑1 (MQ=255)
                                                                                                                                        aGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCACTg      >  1:1048145/1‑150 (MQ=255)
                                                                                                                                           cTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGCGGCAAAAAGTAACGCGGGCCAGTTTCACCCTCCAGCCCAAGACGTTGAGTGTGCTTGCTTATAAACCCCGCCCCCAGGTTCGAACCTTCACTGCaa   <  1:292872/150‑1 (MQ=255)
                                                                                                                                            tGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCACTGCAAt  >  2:1057726/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CTGTGCAGCACCCGCTAAAGGTGCCATTGCGCTCGGTCAACTGCCTGCAAATACCTCACGCAGATGCACCCGATAACTGTCGAACGCTGCAATCGCCCCCGCTACCACCTCTGCTTCCTCCTGCTCCCCCAGCACCAGGCTGTCCAGCCTGGCAACAAATTCGCGCCAATGCAGCGCTCGGCCATCCGGGTGAGGGGCCAGGTGGCGGGCCCCATTGTCACCCTCCAGCCCAAGACGTTGCGTGTGCTTGTAGAGAAACGCCGCCCCCAGGTTCGAACCTTCACTGCAAT  >  NC_002947/2948498‑2948787

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