Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,817,948 G→T 5.4% E63D (GAG→GAT dinD → DNA damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,9480GT 5.4% 1787.1 / 65.4 507E63D (GAG→GATdinDDNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base G (274/205);  new base T (14/13);  total (288/219)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-03

Only 100 of 533 total aligned reads displayed.

ACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAGC  >  NC_000913/3817682‑3818220
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                
                                         ctACTTTCTAAGCAATTGATATCACATGGAGTGGGCAATGAACGAATATCATCAACCTATTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGCACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTAGACCTGTATTAACGCGCGCCGACGAAGCCTGTGGGAACCGTGGTCACACAAATGATGACCAGTTCGAGGAGCTCCACGATACGGTCAAAATTGGCTCAAATGCAAAAAGAGCATTaaaa                                                                                                                                                                                                                                  >  4:243229/2‑274 (MQ=255)
                                                   tGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCAAAAAGAGCATTAAAACACATCGt                                                                                                                                                                                                                          >  2:514855/1‑272 (MQ=255)
                                                          gATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAtt                                                                                                                                                                                                                                      >  2:202328/1‑253 (MQ=255)
                                                              tCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaaa                                                                                                                                                                                                                                  <  4:202327/253‑1 (MQ=255)
                                                                aactgggatGGGAAAAGAACCAACTGCATNGACCTTTTGAAGAGATAAAACTGATNAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGTGTTATTCAGAGTAACGTCACTTTATACCNGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGGGGTCACACAAATGATGACCATTACGGGGAGATCATCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                                                                                                                                                          <  3:87359/267‑1 (MQ=255)
                                                                aaaTGGAGTTGGCAAATAACGAACATCATCAACCTTTTGAAGAAATAAAACGGATTAATGCAAACGGAGCATAAAAATGGTTAGTAATACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGGGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                                                                                                                                                          <  4:438998/273‑1 (MQ=255)
                                                                          ggCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAACGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                                                                                                                                                                >  2:6927/1‑275 (MQ=255)
                                                                          ggCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                                                                                                                                                                >  4:291685/1‑275 (MQ=255)
                                                                            cagtGAACGAACATCTTCAACCTTTCGAAGAGGTAAAACTTATTAATGCAAACGTAGCAGAACAATGGTCAGCAAGCCAACCTAGGAAACTACTAGGTTATTCAGAGTATCGTTACTTAATAACTGTATTAACGCTCGCCAAAGACGCCTGTGAGAACAGTGGTCACACAATTGATGACGATTCCGCGGAGATCCGCGAAGTCGCCAAACTTGTCCCAAATCCCAACACGGTATTAGAGGACACCGTACTCCCCCCCTATTCCTTGTACCTTgt                                                                                                                                                                                               >  4:97757/4‑274 (MQ=255)
                                                                             aaTGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                   <  4:146761/237‑1 (MQ=255)
                                                                              aTGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCgaggag                                                                                                                                                                                                                                                                                  <  2:553493/189‑1 (MQ=255)
                                                                              aTGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTa                                                                                                                                                                                                     <  2:178323/266‑1 (MQ=255)
                                                                              aTGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGAGATCCTCGat                                                                                                                                                                                                                                                                         <  2:317994/198‑1 (MQ=255)
                                                                                gaacgaacATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAt                                                                                                                                                                                                                                       >  2:241433/1‑230 (MQ=255)
                                                                                                            aTAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTgg                                                                                                                                                                                                                                                           >  4:109560/1‑182 (MQ=255)
                                                                                                               anaCTGATTAATGCAAACGGAGCAGAACAATTGTCANCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAag                                                                                                                                                                                                                                            <  4:554451/194‑1 (MQ=255)
                                                                                                                      ttAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcag                                                                                                                                                      >  1:121846/1‑273 (MQ=255)
                                                                                                                                             tGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCa                                                                                                                                                                                                                                                >  1:25962/1‑160 (MQ=255)
                                                                                                                                                 cAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTc                                                                                                                                                                                                                                                                   >  6:789708/1‑137 (MQ=255)
                                                                                                                                                 cAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTc                                                                                                                                                                                                                                                                   <  5:789708/137‑1 (MQ=255)
                                                                                                                                                  aGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAt                                                                                                                                                                                                                                       >  1:67621/1‑164 (MQ=255)
                                                                                                                                                    cAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTccc                                                                                                                                                                                                                  >  2:570115/1‑183 (MQ=255)
                                                                                                                                                        aCAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCg                                                                                                                                                                                                                                                                           <  7:968133/122‑1 (MQ=255)
                                                                                                                                                           aCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTc                                                                                                                                                                                                                                                                            <  8:1044376/118‑1 (MQ=255)
                                                                                                                                                           aCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTc                                                                                                                                                                                                                                                                   <  7:637914/127‑1 (MQ=255)
                                                                                                                                                              tGGGAAACTACTGGGTTATTCCGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGAt                                                                                                                                                                                                                                                                                <  7:1079274/111‑1 (MQ=255)
                                                                                                                                                                gaaaactccggttttATTCATATTTTCGTCACTTTATACCTGTATTAACGCGGTGCAAAGAAGCCTGTGAGAACAGTNGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATTCCAAAAGAGCATTAACAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCt                                                                                                          <  3:121846/262‑1 (MQ=255)
                                                                                                                                                                       actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAagag                                                                                                                                                                                                                                          >  6:195526/1‑140 (MQ=255)
                                                                                                                                                                        ctGGGTTATTCAAAGGATCGTCCCTTTATACCTGTATTAACCCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGa                                                                                                  <  4:58570/275‑1 (MQ=255)
                                                                                                                                                                         tGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGa                                                                                                                                                                                                                                <  8:610825/148‑1 (MQ=255)
                                                                                                                                                                                ttCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCg                                                                                                                                                                         >  4:361601/1‑196 (MQ=255)
                                                                                                                                                                                ttCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAt                                                                                                                                           >  4:318627/1‑226 (MQ=255)
                                                                                                                                                                                 tCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATgg                                                                                                                                                                                                                                                                     >  5:1064631/1‑103 (MQ=255)
                                                                                                                                                                                    gagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                   >  8:1080127/1‑134 (MQ=255)
                                                                                                                                                                                    gagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTa                                                                                                                                            <  2:144260/221‑1 (MQ=255)
                                                                                                                                                                                     agTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCaa                                                                                     >  4:504464/1‑275 (MQ=255)
                                                                                                                                                                                     agTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCaa                                                                                     >  2:47428/1‑275 (MQ=255)
                                                                                                                                                                                            tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                   >  5:122607/1‑126 (MQ=255)
                                                                                                                                                                                            tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTAt                                                                                                                                                                                                             >  8:592821/1‑148 (MQ=255)
                                                                                                                                                                                             cACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCt                                                                                                                                                                                                                                                                             <  7:611463/83‑1 (MQ=255)
                                                                                                                                                                                                tttATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCg                                                                                                                                                                                                                                                                           <  8:121098/82‑1 (MQ=255)
                                                                                                                                                                                                tttATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTAGCGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcca                                                                                                                                                                                                         >  7:498736/1‑147 (MQ=255)
                                                                                                                                                                                                        tGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAcc                                                                                                                                                                                                   <  6:387815/146‑1 (MQ=255)
                                                                                                                                                                                                        tGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTg                                                                                                                                                                                                 >  7:599892/1‑148 (MQ=255)
                                                                                                                                                                                                          tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGTAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTTAAGACGAAAAGCGt                                                                >  2:390672/1‑275 (MQ=255)
                                                                                                                                                                                                           aTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGa                                                                                                                                                                                                                                                                          <  8:528210/72‑1 (MQ=255)
                                                                                                                                                                                                           aTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCAGTTCGAGGAGATCCTCGata                                                                                                                                                                                                                                                                        >  6:965775/1‑74 (MQ=255)
                                                                                                                                                                                                                cgcgcgCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAag                                                                                                                                                                                                                                            >  7:431675/1‑97 (MQ=255)
                                                                                                                                                                                                                cgcgcgCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAag                                                                                                                                                                                                                                            <  8:431675/97‑1 (MQ=255)
                                                                                                                                                                                                                cgcgcgCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAcc                                                                                                                                                                                                   <  5:479550/138‑1 (MQ=255)
                                                                                                                                                                                                                  cgcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGc                                                                                                                                                                                                                                         >  5:210930/1‑98 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAg                                                                                                                               >  2:442022/1‑202 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTga                                                      >  2:259283/1‑275 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTga                                                      >  4:390298/1‑275 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTga                                                      >  4:183503/1‑275 (MQ=255)
                                                                                                                                                                                                                      ccAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATg                                                                                                                                                                                                            <  6:909877/123‑1 (MQ=255)
                                                                                                                                                                                                                        aaaGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACtga                                                                                                                                            <  4:213444/185‑3 (MQ=255)
                                                                                                                                                                                                                         aagtagCCGTTGAGAAAAGTGGTCAAACAATGGAGGACCATTTCGGGGAGATCCTCGCTAAGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCGATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGGGAAGACGAAAAGCGTCTGTTTCTGAGAAAt                                                 <  4:412440/275‑1 (MQ=255)
                                                                                                                                                                                                                         aagaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAACTTGGCt                                                                                                                                                                                                                                                         <  8:716028/75‑1 (MQ=255)
                                                                                                                                                                                                                         aagaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                                                                                                                                                         >  7:716028/1‑75 (MQ=255)
                                                                                                                                                                                                                          agaagCCTGTGAGAGCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAc                                                                                                                                                                                                    >  8:844112/1‑127 (MQ=255)
                                                                                                                                                                                                                               ccTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTg                                           >  2:217767/1‑275 (MQ=255)
                                                                                                                                                                                                                                 tgtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATgg                                                                                                                                                                                                                                                                     >  6:476326/1‑55 (MQ=255)
                                                                                                                                                                                                                                  gtgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaaa                                                                                                                                                                                       <  7:1022946/132‑1 (MQ=255)
                                                                                                                                                                                                                                   tgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGATTTGAAAg                                       >  1:94364/1‑275 (MQ=255)
                                                                                                                                                                                                                                    gagaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                                                                                                                                                                                                            >  8:906548/1‑61 (MQ=255)
                                                                                                                                                                                                                                    gagaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                                                                                                                                                                                                            <  7:906548/61‑1 (MQ=255)
                                                                                                                                                                                                                                       aaCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTg                                                                                                                                                                                                 <  8:644985/117‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCa                                                                                                                                                                                                                                        >  7:84251/1‑74 (MQ=255)
                                                                                                                                                                                                                                           gTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATa                                 <  2:423106/273‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGGTCACACAATTGATGACCATTCCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaa                                                                                                                                                                                                                                    >  5:281166/1‑78 (MQ=255)
                                                                                                                                                                                                                                            tGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATaaca                              >  2:336115/1‑275 (MQ=255)
                                                                                                                                                                                                                                             ggTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                   >  8:756841/1‑77 (MQ=255)
                                                                                                                                                                                                                                             ggTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                   <  7:756841/77‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAg                                       >  4:339865/1‑265 (MQ=255)
                                                                                                                                                                                                                                                  cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAcc                                                                                                                                                                                                   <  7:202012/104‑1 (MQ=255)
                                                                                                                                                                                                                                                  cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggca                                                                                                                                                       <  8:247593/148‑1 (MQ=255)
                                                                                                                                                                                                                                                    caATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGAc                                                                                                                                                                                                                               >  7:1139806/1‑74 (MQ=255)
                                                                                                                                                                                                                                                    caATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCg                                                                                                                                                                                                                           >  5:971840/1‑78 (MQ=255)
                                                                                                                                                                                                                                                    caATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagg                                                                                                                                                     >  5:67023/1‑148 (MQ=255)
                                                                                                                                                                                                                                                      aTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                   >  6:119515/1‑68 (MQ=255)
                                                                                                                                                                                                                                                      aTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGaaaa                                                                    >  2:358075/1‑227 (MQ=255)
                                                                                                                                                                                                                                                      aTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGt                    >  2:366729/1‑275 (MQ=255)
                                                                                                                                                                                                                                                           tgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAgg               <  2:291686/275‑1 (MQ=255)
                                                                                                                                                                                                                                                            gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGt                                                                                                                                                                 <  6:265008/128‑1 (MQ=255)
                                                                                                                                                                                                                                                               cATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                                                                                                                                                          >  5:426083/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                atatCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaaa                                                                                                                                                                                       >  6:608981/4‑102 (MQ=38)
                                                                                                                                                                                                                                                                atatCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaaa                                                                                                                                                                                       <  5:608981/99‑1 (MQ=39)
                                                                                                                                                                                                                                                                 tttCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGcagcac         >  4:264080/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 tatCGAGGATCTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGc                                                                                                                                                                                                                >  5:924512/3‑76 (MQ=38)
                                                                                                                                                                                                                                                                  ttCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCTAAAAGAGCATTAAAAGACAt                                                                                                                                                                                                                             >  6:51145/1‑62 (MQ=255)
                                                                                                                                                                                                                                                                  ttCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACAt                                                                                                                                                                                                                             <  5:51145/62‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCGCGGAGATCCTCGATATGGTCAANATTGGCNCAAATACCAAAAGGGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACgg                                                                                                                      <  3:26971/164‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cGAGGATCTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt                                                                                                                                          <  7:932731/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cGAGGATCTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt                                                                                                                                          >  8:932731/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                      aggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                                                                                                                                                                >  6:1060916/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                      aggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTAt                                                                                                                                   <  6:67023/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ggagTTCCTCGATATGGTCAAAATTGGATAAAATGCCAAAAGAGAATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACGTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAg   <  4:472853/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gagaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCg                                                                                                                                                          <  8:360520/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gagaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCg                                                                                                                                                          >  7:360520/1‑123 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                
ACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAGC  >  NC_000913/3817682‑3818220

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