Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,099,027 A→G 7.2% intergenic (+71/+27) tpm → / ← htpX thiopurine S‑methyltransferase/protease HtpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,099,0270AG 7.2% 861.6 / 35.6 278intergenic (+71/+27)tpm/htpXthiopurine S‑methyltransferase/protease HtpX
Reads supporting (aligned to +/- strand):  ref base A (142/116);  new base G (8/12);  total (150/128)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 281 total aligned reads displayed.

TGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCATGT  >  NC_002947/2098878‑2099175
                                                                                                                                                       |                                                                                                                                                    
tgtgATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGaa                                                                                                                                                      >  2:3303053/1‑150 (MQ=255)
      tGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGcc                                                                                                                                                <  2:1443118/150‑1 (MQ=255)
         ggACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCtt                                                                                                                                             >  2:6339144/1‑150 (MQ=255)
          gACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACTTCGGTCGCCCttt                                                                                                                                            <  1:3534016/150‑1 (MQ=255)
          gACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCttt                                                                                                                                            <  1:3742881/150‑1 (MQ=255)
           aCAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACTTCGGTCGCCCtttt                                                                                                                                           <  2:5636004/150‑2 (MQ=255)
           aCAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACTTCGGTCGCCCtttt                                                                                                                                           <  1:5619646/150‑2 (MQ=255)
           aCAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTg                                                                                                                                           >  2:4701908/1‑150 (MQ=255)
                ncGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGttttt                                                                                                                                      >  1:6380734/2‑150 (MQ=255)
                ccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGttttt                                                                                                                                      >  2:4532813/1‑150 (MQ=255)
                ccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGttttt                                                                                                                                      <  1:6318735/150‑1 (MQ=255)
                 cGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTa                                                                                                                                     <  2:1751903/150‑1 (MQ=255)
                  gAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                    >  1:1052080/1‑150 (MQ=255)
                   aattttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATc                                                                                                                                   <  1:1848975/147‑1 (MQ=255)
                   aaGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATc                                                                                                                                   <  1:4145295/150‑1 (MQ=255)
                     gTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAg                                                                                                                                 <  1:3167938/150‑1 (MQ=255)
                      tttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAgg                                                                                                                                <  1:6000432/150‑1 (MQ=255)
                       ttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAACAAAGGGCGACTTCGGTCGCCCTttttttattcccg                                                                                                                                 <  2:6619913/150‑8 (MQ=37)
                        tgtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGt                                                                                                                               >  1:1065272/1‑149 (MQ=255)
                        tgtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGt                                                                                                                               <  2:1065422/149‑1 (MQ=255)
                        tgtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCGTTGTTTTTATCAGGTg                                                                                                                              >  1:3784941/1‑150 (MQ=255)
                         gtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTgg                                                                                                                             <  2:3197042/150‑1 (MQ=255)
                                  ggggCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCcgc                                                                                                                    <  2:6197669/150‑1 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AATGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                              >  1:2628191/1‑150 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                              <  2:1565807/150‑1 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                              >  2:822941/1‑150 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                              >  1:6093666/1‑150 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                              <  2:3564650/150‑1 (MQ=255)
                                        tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGAc                                                                                                              >  1:3611458/1‑150 (MQ=255)
                                          gAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGc                                                                                                            <  1:6665869/150‑1 (MQ=255)
                                          gAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCa                                                                                                           >  2:2766717/1‑150 (MQ=255)
                                            gTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTAACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAg                                                                                                          <  2:6786816/150‑1 (MQ=255)
                                             tttGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGttttt                                                                                                                                      >  2:623281/1‑121 (MQ=255)
                                             tttGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGttttt                                                                                                                                      <  1:623281/121‑1 (MQ=255)
                                             tttGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAgg                                                                                                         >  2:4862602/1‑150 (MQ=255)
                                             tttGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAgg                                                                                                         <  2:6483912/150‑1 (MQ=255)
                                                      gCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AATGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCCGCCGCGCGGACGCTGGGCCTCGATg                                                                                                >  1:5703483/1‑150 (MQ=255)
                                                       cGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGc                                                                                               <  2:4921695/150‑1 (MQ=255)
                                                         ggTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCgg                                                                                             >  1:5292223/1‑150 (MQ=255)
                                                          gtgtATCAGGTCAGCCGATAAGCAGGTATGGCAGCCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGt                                                                                            >  1:4231361/1‑150 (MQ=255)
                                                               tCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTcc                                                                                       >  1:1455125/1‑150 (MQ=255)
                                                                 aGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAg                                                                                     <  2:37115/150‑1 (MQ=255)
                                                                 aGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAg                                                                                     >  2:5240256/1‑150 (MQ=255)
                                                                  ggTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGc                                                                                    <  1:3397486/150‑1 (MQ=255)
                                                                    tCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCgg                                                                                  >  2:5434535/1‑150 (MQ=255)
                                                                     cAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                                 <  2:1699128/150‑1 (MQ=255)
                                                                       gCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgg                                                                               <  1:5154553/150‑1 (MQ=255)
                                                                         cGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTTCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggt                                                                             >  2:6160101/1‑150 (MQ=255)
                                                                           aTAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGgtgggtgg                                                                           <  2:6031694/150‑1 (MQ=255)
                                                                             aaGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCt                                                                         >  1:1707534/1‑150 (MQ=255)
                                                                               gCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtca                                                                       >  2:4199577/1‑150 (MQ=255)
                                                                               gCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtca                                                                       >  1:5244100/1‑150 (MQ=255)
                                                                                cAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcag                                                                      >  1:4587934/1‑149 (MQ=255)
                                                                                cAGGTAGGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcat                                                                      <  1:943975/150‑1 (MQ=255)
                                                                                 aGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatc                                                                     <  2:569155/150‑1 (MQ=255)
                                                                                  ggTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatca                                                                    <  2:2472890/150‑1 (MQ=255)
                                                                                  ggTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatca                                                                    <  1:4919394/150‑1 (MQ=255)
                                                                                            gCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTTGCTGATGAGCCAcccacc                                                          <  2:753204/150‑4 (MQ=255)
                                                                                             cTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGcc                                                         >  2:5368161/1‑150 (MQ=255)
                                                                                              tGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCa                                                        >  2:5088833/1‑150 (MQ=255)
                                                                                              tGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCa                                                        >  1:3480996/1‑150 (MQ=255)
                                                                                              tGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCa                                                        >  1:67703/1‑150 (MQ=255)
                                                                                                  ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                    >  2:1957534/1‑150 (MQ=255)
                                                                                                  ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                    >  1:1403796/1‑150 (MQ=255)
                                                                                                  ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                    >  1:4650696/1‑150 (MQ=255)
                                                                                                  ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGcc                                                    >  2:3843913/1‑150 (MQ=255)
                                                                                                       cccACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AATGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCt                                               >  1:4487272/1‑150 (MQ=255)
                                                                                                        ccACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                              >  1:5359709/1‑150 (MQ=255)
                                                                                                        ccACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                              >  1:530741/1‑150 (MQ=255)
                                                                                                        ccACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                              >  2:4465333/1‑150 (MQ=255)
                                                                                                         cacaGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTg                                             >  2:2122783/1‑150 (MQ=255)
                                                                                                          acaGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGa                                            >  2:4378390/1‑150 (MQ=255)
                                                                                                          acaGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGa                                            >  2:6609204/1‑150 (MQ=255)
                                                                                                           caGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAg                                           <  1:3957879/150‑1 (MQ=255)
                                                                                                            aGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAgg                                          >  2:2378603/1‑150 (MQ=255)
                                                                                                              aaGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgcc                                        <  1:1496679/150‑1 (MQ=255)
                                                                                                                ggCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGgccgc                                      >  1:3761207/1‑150 (MQ=255)
                                                                                                                    gCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCAtt                                  <  1:990853/150‑1 (MQ=255)
                                                                                                                      gTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGa                                <  2:931959/150‑1 (MQ=255)
                                                                                                                       tAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGAt                               <  2:6410342/150‑1 (MQ=255)
                                                                                                                        aaGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATg                              >  1:3381071/1‑150 (MQ=255)
                                                                                                                            ccGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGa                          >  1:1179525/1‑150 (MQ=255)
                                                                                                                             cGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGaa                         <  2:5292223/150‑1 (MQ=255)
                                                                                                                              gggAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAg                        <  2:5808447/150‑1 (MQ=255)
                                                                                                                              gggAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAg                        <  1:6354466/150‑1 (MQ=255)
                                                                                                                               ggAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAgg                       >  2:6747330/1‑150 (MQ=255)
                                                                                                                                ataaaaac‑‑aaaGGGCGACTTCGGTCGCCCTTTTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGcc                     >  1:4414697/5‑150 (MQ=38)
                                                                                                                                    aaac‑‑aaaGGGCGACTTCGGTCGCCCTTTTTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAAccc                                                             >  1:6179002/1‑107 (MQ=38)
                                                                                                                                    aaac‑‑aaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTc                  >  2:1274242/1‑150 (MQ=255)
                                                                                                                                      aa‑‑aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt                >  1:1143235/1‑150 (MQ=255)
                                                                                                                                      aa‑‑aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt                <  1:2283584/150‑1 (MQ=255)
                                                                                                                                      aa‑‑aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCAt                >  1:5148070/1‑150 (MQ=255)
                                                                                                                                       a‑‑aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg               <  1:2766717/150‑1 (MQ=255)
                                                                                                                                          aaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCGGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATgg              >  1:810486/1‑150 (MQ=255)
                                                                                                                                           aaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGgt             <  1:5383859/150‑1 (MQ=255)
                                                                                                                                                 gACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCgg       >  1:238548/1‑150 (MQ=255)
                                                                                                                                                 gACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCgg       >  2:3643115/1‑150 (MQ=255)
                                                                                                                                                   ccGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCa     >  2:1671426/1‑150 (MQ=255)
                                                                                                                                                    cGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCAt    <  1:2424491/150‑1 (MQ=255)
                                                                                                                                                     gAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCAtg   >  1:689091/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
TGTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCATGT  >  NC_002947/2098878‑2099175

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