Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 388,319 C→G 16.3% R175R (CGC→CGG tauD → taurine dioxygenase, 2‑oxoglutarate‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913388,3190CG16.3% 167.9 / 2.4 73R175R (CGC→CGGtauDtaurine dioxygenase, 2‑oxoglutarate‑dependent
Reads supporting (aligned to +/- strand):  ref base C (44/12);  new base G (2/11);  total (47/26)
Fisher's exact test for biased strand distribution p-value = 3.78e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCC  >  NC_000913/388098‑388554
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
gtgACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGCGGTAAGCCCTGGGGCAGCGGGCTTACCGGCCATGAGGCGGCCTATTTTTCCTTCCGCCCGCGTCCGGGTGGGGTGCGGGGGGAGCATGGATTTCCGTAAACCTTTCCCGAATAACAATAACGCAAAACCCAGGAGGGACATCCACCCATGCGGGCGGCGGGCGGCGaaaa                                                                                                                                                                                                                                  >  2:38600/1‑239 (MQ=255)
     aTTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGGTACGCCCTTGGCCCGCCGTATTTTGGCCTTTGAGGCGGTCTCTCTTTCCTTCCCCCCGCCGCTGAGTGGGCTGGCGGCGGGGGCTGGTTTTCCTTAATCCATTCCGGGCTACCACATACCCAAAACCCACGAGGAACAACAACCCTGGCCGGGGGCGGTGGCGaaaaaaacc                                                                                                                                                                                                                             >  2:372987/1‑234 (MQ=255)
                                                     aCCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGGAATGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGGGTGGGCTGGGGGCGGGGCCTTATTTCCCTAAAACGGTTCCGGAATACAAATACCGCCAAAACCAGGAGGAACATCCACCCTGGGCAGCGGCGGGAGAGCAAAACCCGCCGCTGGCACACCCGGCGGGGCGCAcacgccccggcagagggaaaaacg                                                                                                                                                                           >  2:367770/1‑216 (MQ=255)
                                                                    ttgttctctctGGGCCGGGGGTTTTGCGCCTTATGGGGCGCTTTCTGTCCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGGGGAGCATGATTCCCGTAAATCGTCCCCGGAATACAAATCCCGAAAAACCGAGGAGAAACACCAACGCGGGCGGGAGGGGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGtgaatg                                                                                                                                                         <  2:1299756/232‑1 (MQ=255)
                                                                      agcgcctctGGACCAGGTGTTGTGAGGGCGTTGGGGCGCTCTCTTTTCCCTTCGCCCGGCTGTGGGGGGGGTCGGGCGGGGATCATGTTTCCCGTAAATCGTTCCGGAAAAACAAAGACCGAAAACCCGGGGAGGAACAACAACCGTGGCGAGAGGCGGTCGAGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGtgaatga                                                                                                                                                        <  2:1219549/233‑1 (MQ=255)
                                                                              tGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACAACCGGTGGTGCGAACGCATCCGGTGACCGGTAAACCGGCGCTGTTTGTGAATGAAGGCTTtat                                                                                                                                               >  1:256765/1‑238 (MQ=255)
                                                                                   cgggggTTTTGGGGCCTTGGGGGCCTTCTTTGTTCCCCTCCGCAAGTGGCTGGGGGGGCTGGGGGGGGAGCTTGATTCCCGAAAATCGTCCCCGGAAAACAAATCCCGAAAGACCGAGGAGGAACATAAACGCTGGGGCGGGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTAcgc                                                                                                                                          <  2:353209/235‑1 (MQ=255)
                                                                                     gCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCt                                                                                                                                                                                                                    <  1:810553/166‑1 (MQ=255)
                                                                                     gCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCt                                                                                                                                                                                                                    >  2:810553/1‑166 (MQ=255)
                                                                                      cGGTATTGCGGCCTATGAGGCGCTCTTTGTTCCCTCCCGCCGGCTGCTGAGTGGGCTGCGGGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATCCCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGt                                                                                                                                                                                         <  1:681528/192‑1 (MQ=255)
                                                                                      cGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGACAAACCCGCCGTTGCTAACTCCGGGGGTGCGAAAGCATCCGGGGAGCGGTAAAAAGGCGCGGTTTGTGAATGGAGGCTTtaacaggggaa                                                                                                                                       >  2:1279787/1‑230 (MQ=255)
                                                                                      cGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGATCGGTAAACAGGCCGTGTTTGTGAATGAAGGCTTTACAAcgttct                                                                                                                                       >  1:1219549/1‑235 (MQ=255)
                                                                                          aTTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCCAACGCATCCGGTGAGGGGTAACAAAGGCCGGGTTTTGGATGAAGGCTTtactacccacattgg                                                                                                                                   >  2:561400/1‑230 (MQ=255)
                                                                                          aTTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGCTGTGCGAACGCATCCGGTGAGCGGTAAACAGGGGCTGTTTTGGAATGAAGGCGTtaatatgtgcaatgt                                                                                                                                   >  1:453391/1‑229 (MQ=255)
                                                                                            tGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTtgt                                                                                                                                                              >  1:124894/1‑210 (MQ=255)
                                                                                                   tatactgagcTATCTGTTCCCTCCCGCCACCTGATGACTGGGCTGCGTGCGGAGCATGATTTCCGGACAGCGTTCCCGGAATAAAAATACCGCAAAACCGAGTAGGAACATCAACGCTGGCGCGAGGCTGTCCCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCAAACGCACCGGGTGGGGGGGAAGAAGGCGACTGTTGTGCATGAGGGCTTTAACACGCGAATGGGTGaagggat                                                                                                                          >  2:139904/9‑233 (MQ=255)
                                                                                                   tATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAg                                                                                                                          >  1:1311401/1‑239 (MQ=255)
                                                                                                   tATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAg                                                                                                                          >  1:1159711/1‑239 (MQ=255)
                                                                                                     tGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGACAAACCAGCCGTTGCTACATCCGGTGGGGCGAACGCAACCGGTGAGGGGTAAACAGGCGCTGTTTTTGTATGAAGGCTTTACACGGCGAAttgttgttggtggcg                                                                                                                        >  2:609194/1‑233 (MQ=255)
                                                                                                                 ttttCCCTCCCGCCGGGGGGTGGGTGGGGTGGGGGCGGAGCATGTTTACCGAAAATCGTACCCGGAACACAAATCCCGAAAAACCGGGGAGGAACACCAACGGGGGCGGGGGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCg                                                                                                           <  2:569888/237‑1 (MQ=255)
                                                                                                                   ttCCCTTCCGCCAGCTGCTGAGTGGGCTGCGGGCGGAGCATGATTTCCGTAAATCGTTCCCGAAATACAAATCCCGCAAAACCGAGGAGGAACATCAACGCGGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGaa                                                                                                         <  1:375749/239‑1 (MQ=255)
                                                                                                                        ttCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGGCtg                                                                                                    >  1:548005/1‑238 (MQ=255)
                                                                                                                         tCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAAAGAGCGAAGCCTTg                                                                                                   >  1:732532/1‑239 (MQ=255)
                                                                                                                         tCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCCCGAAAAACCCGCCGTTGCTACATCCGGTGGGCCGAACGCGTGCGGGGAGCGGTAAACAGGCGCCGTTTTGTGAGGAAGGCTGTCCTACGCGAATTGTTGCTTGGAGCGACAAagagcacgaggacttt                                                                                                   >  2:39049/1‑227 (MQ=255)
                                                                                                                           cGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGCAAAACCCGCCGTTGCTACATCCGGGGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGGTTGTGACTGAAGGCTTTACTCCGCGCATTGTTGATGTGGGCGAGAagagaggcgagcccgtgat                                                                                                >  2:870272/1‑226 (MQ=255)
                                                                                                                           cGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTTAGCGGGAAAAGGAGCGAAGCCTTGt                                                                                                  >  1:569888/1‑238 (MQ=255)
                                                                                                                           cGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGGGAGCGAAGCCTTGtt                                                                                                 >  1:123999/1‑239 (MQ=255)
                                                                                                                           cGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTAAGCGACTTGTTTATGTGGGCGGGAAAGAGAGCGAAGCATtttt                                                                                                 >  2:999807/1‑239 (MQ=255)
                                                                                                                                tcttgtgggTGGGGTGCGGGGGGGACAGGTTTTCCGTAAATCGGTCCCGGAATACAAATACCGAAAAACGGAGGAGGAACATCAACGCTGGCGCGGGGCGGGCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGtt                                                                                            <  2:1159711/231‑1 (MQ=255)
                                                                                                                                gctgctGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAACTCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGGAGGCTTTACTACGCGAATTGTTGAGGTGAGCGGGAAAGAGGGCGAAGCCATGTTAAGtt                                                                                            >  2:915643/1‑239 (MQ=255)
                                                                                                                                gctgctGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAggt                                                                                            >  1:1062139/1‑237 (MQ=255)
                                                                                                                                              tgcgtgcgGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATCCCGCAAAACCGAGGCGGAACATCAACGCTGGCGCGAGGCGGCCGCGCAAAACCCGCCGTTCCTACATCCGGTGGTGGGACCGTATCCGGTGAGCGGGAATCCGGCCCCGTTTCTGAATGAACGTTTTACTACGCGCATTTCTGCTGTGTGCGTGCAGCAGGGGGAAGACTTGTTAAATTTTTTTTTTGtccca                                                                              >  2:271188/1‑237 (MQ=255)
                                                                                                                                              tgcgtgcgGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGGGCGAACGCATCCGGGGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCat                                                                              >  2:1008279/1‑239 (MQ=255)
                                                                                                                                                cgtgcgGAGCCTGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGTTCGCGAAAAACCCGCCGTTGCTACATACGGGGATGCGAACGCATCCGGTGCGCGTGTAACAGGCGCCGTTGGTGAATCAAGGCTTTACTATGGTAATTGTTGAGGTGGGGGGGGAAGGAAGGGAGGCGATCTTAAGTTTTTTGTTTGACAatac                                                                            >  2:1242190/1‑238 (MQ=255)
                                                                                                                                                 gtgcgGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGGGAATGAAGGCTTTACTACGCGAATTGTTGATGTGGGGGGGAAAGAGAGGGAAGCCTTGTTAAGTTTTTTTTTTGCCCAAATc                                                                           >  2:315456/1‑239 (MQ=255)
                                                                                                                                                 gtgcgGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGCGAGCGGTAAACAGGCGCTGCTTGTGAATGCAGGCTTTACTACCGGACATCGTGATGTGAGCGGGAAGGAAGAGGAAGCCTTGTTAAGTTTTTTTTTTGCCCaaaaa                                                                           >  2:420968/1‑235 (MQ=255)
                                                                                                                                                          aTGATTCCCGAAATTCGTTCCGGGAATACAAATCCCACAAACCGGAGGAGAAACATCAACGCGGGGGCGAGGCGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAAcc                                                                   <  1:315456/238‑1 (MQ=255)
                                                                                                                                                          aTGATTCCCGAAAATCGTCCCCGGAAACCAAATCCCGAAAACCCGAGGAGGACCATAAACGCTGGCGCGAGGCGGCGGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAAcc                                                                   <  1:322152/238‑1 (MQ=255)
                                                                                                                                                                 tCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCCTCCGGTGAGCGGTAAACCGACGCTGTTTGTGAGTGAAGGCGTTACTACGCGAATTGTTGATGTGAGCGAGAAGGAGAGCGAGGCCTTGTTGAGTTTTTTGTGTGCCCAAATCGCCAAACCGGAGttg                                                            >  2:1058005/1‑238 (MQ=255)
                                                                                                                                                                    gaaaaTGGTTCCGGGAATACAAATCCCGAAAACCCGGGGAGGAACATAACCGCGGGGGGGGGGGGGTCGCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAg                                                         <  1:561400/237‑1 (MQ=255)
                                                                                                                                                                    gTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTTAGCGGGAAAGAGACGGAAGCCTGGTTAATTTTTTTTTTTGCCCATATAACCAAACAGGGGTTTCAg                                                         >  2:84022/1‑239 (MQ=255)
                                                                                                                                                                             tCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTgcgctg                                                 >  1:943716/1‑238 (MQ=255)
                                                                                                                                                                             tCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGGgcgctg                                                 >  1:1354955/1‑238 (MQ=255)
                                                                                                                                                                               ccGGAAAAAAAATCCCACAAAACCGGGGGGAAAAACAAAGCGGGGCGGAGGGGGGCCCGAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTgcgctggcg                                              <  2:1354955/239‑1 (MQ=255)
                                                                                                                                                                                    acaCAAAGCCCGAAAAACGGGGGAGAAACCACCACTGTGGCGAGGGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGcgctggcgctgg                                          <  1:1206123/236‑1 (MQ=255)
                                                                                                                                                                                    aTACAAATACCCAAAACGCGGGGAGAACCAGAAACGCTGGCGCGAGGGGGTCGAAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCgctggcgctggc                                         <  1:925972/239‑1 (MQ=255)
                                                                                                                                                                                    aTAAAAATACCGCAACACGGGGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGcgctggcgctgg                                          <  1:356670/238‑1 (MQ=255)
                                                                                                                                                                                            tcccccAAAACACAGGAGGCACATCAAGGGTGGCGGGAGGCGGTCGGAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTCGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCaaa                                  <  1:1058005/234‑1 (MQ=255)
                                                                                                                                                                                              ccGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTAc                                                                                                                                                                                                                  <  1:1095068/65‑1 (MQ=255)
                                                                                                                                                                                              ccGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTAc                                                                                                                                                                                                                  >  2:1095068/1‑65 (MQ=255)
                                                                                                                                                                                              ccGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGGGCTGGCAACAAAATg                                >  1:431563/1‑239 (MQ=255)
                                                                                                                                                                                              ccGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGGGGCCTGGGCCTGGCAACCaaaga                                >  2:661588/1‑237 (MQ=255)
                                                                                                                                                                                                 cAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGGGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCCCTGGCCACCAAACGata                             >  1:614804/1‑239 (MQ=255)
                                                                                                                                                                                                 cAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGCAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGCCGCTGGCAACCAAATGata                             >  1:217565/1‑239 (MQ=255)
                                                                                                                                                                                                 cAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCCAACCGGAGTTTCAGGTGCGTTGGCGCTGGCAACCacacaaaa                             >  1:1037688/1‑232 (MQ=255)
                                                                                                                                                                                                 cAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTAATACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGGGTTTCAGATGCGCTGGCGCTGGCAACCAATTGaaa                             >  1:15332/1‑237 (MQ=255)
                                                                                                                                                                                                    aaCCGGGGGGGAAAACAAACGGTGGGGGGGGGGGGGCGCCGAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATAtt                           <  2:943716/239‑1 (MQ=255)
                                                                                                                                                                                                     cccGGGGGGGAACATCCACGGTGGGGGGGGGGGGGGGGGAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGGGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTg                          <  2:123999/238‑1 (MQ=255)
                                                                                                                                                                                                      ccGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGGGCTGGCAACCAAATGATATTGCg                        >  1:774518/1‑239 (MQ=255)
                                                                                                                                                                                                               gaACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGCGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGCGTTTCAGGTGCGCTGGCGCTGGCAACCAAATAAAATTTGGATTTGGGa                >  2:538144/1‑239 (MQ=255)
                                                                                                                                                                                                               gaACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGa                >  1:946928/1‑239 (MQ=255)
                                                                                                                                                                                                                 aCATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATa              >  1:87625/1‑239 (MQ=255)
                                                                                                                                                                                                                       cctgtGGCGTGGGGGGGTCGCAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:372987/235‑1 (MQ=255)
                                                                                                                                                                                                                       aggggggCGGGAGGGGGTCGCAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:999807/234‑1 (MQ=255)
                                                                                                                                                                                                                       accggggCGGGAGGCGGTGTCAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTTAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:735352/234‑1 (MQ=255)
                                                                                                                                                                                                                       aCGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGGTTTACTACGAGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGGTTTTTGTTTGCCCATATCCCCAAACCGGAGTTTCAGGTACGCTGGCGCTGGCAACCAAATGATAATCCGATTAGGGAGAACcgcg        >  2:244623/1‑239 (MQ=255)
                                                                                                                                                                                                                       aCGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGGGACGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGGCGGCTGGCGGTGGCACCCAAAAGATATTGCGATTTGGGAACACcgcg        >  2:343179/1‑239 (MQ=255)
                                                                                                                                                                                                                       aCGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCCCCAAACCGGAGTTTCAGGTGCGCTGGGGCTGGCCAACCAATTCTATTGCGGTTTggggacaccccg        >  2:375869/1‑230 (MQ=255)
                                                                                                                                                                                                                         gcgggCGCGAGGCGGACGCGAACACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:38600/233‑1 (MQ=255)
                                                                                                                                                                                                                         gCTGGGGGGAGGCGGTCGAAAAAACCCCGCCGTTGCTACATCCGGTGGGGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:112170/237‑1 (MQ=255)
                                                                                                                                                                                                                          tgggCGGGGGGGGGTCGGAAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:169983/234‑1 (MQ=255)
                                                                                                                                                                                                                          cTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCCCCAAACCGGAGTTTCAGTGTCGCTGGGGCTGGCAACCAAATAAAATTGCGATTTGGGTAAAcccggga      >  2:288436/1‑233 (MQ=255)
                                                                                                                                                                                                                           tGGCGCGAGGGGGTCGCAAAAAACCCGCCGTTGCTACATCCGGGTGTGCGAACGCAGCCGGGGAGGGGTAAACTGGCGCTGTTTGTGAAGGAAGGCTTTACTACGCGAATTGTTGATGTGATCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:673876/235‑1 (MQ=255)
                                                                                                                                                                                                                             gcgcGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGCATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTTCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGAACCCAAATGATATTGCGATTTTGGAAAACCGCGGGGccc  >  2:1017803/1‑239 (MQ=255)
                                                                                                                                                                                                                             gcgcGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGAt                                                                                                                               >  1:1203118/1‑115 (MQ=255)
                                                                                                                                                                                                                             gcgcGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGCAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGGGTTTCAGGTGGGCTGGCGGCGGGAACCAAATGGTATTGCGAATTGGGATAAACGCGTGAcca  >  2:1218870/1‑238 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
GTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCATGATTTCCGTAAATCGTTCCCGGAATACAAATACCGCAAAACCGAGGAGGAACATCAACGCTGGCGCGAGGCGGTCGCGAAAAACCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAAACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAGAGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCC  >  NC_000913/388098‑388554

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