Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,876,559 G→T 20.5% intergenic (‑230/‑22) cysD ← / → iap sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,876,5590GT20.5% 163.1 / 4.2 70intergenic (‑230/‑22)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (13/41);  new base T (14/2);  total (27/43)
Fisher's exact test for biased strand distribution p-value = 7.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCAGAAATGTTATCTGCCGATTATATTCGCCAACAGTTTCAGCAAATGGGTTATCGCAGT  >  NC_000913/2876322‑2876796
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
tGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGGCATCAGAAGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGc                                                                                                                                                                                                                                                    <  1:320082/239‑1 (MQ=255)
tGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGc                                                                                                                                                                                                                                                    <  1:879711/239‑1 (MQ=255)
 tatcccataCCCTTCCTTTGGAATTCCCCTTATTTCTTGCCATCAGATGTTTCGACTATAGGGTGCGTAAGAGAAAGAATGAAATTACCAATTAGAATGAGTAGTCCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAATTGCa                                                                                                                                                                                                                                                  <  2:1027369/231‑1 (MQ=255)
 taTCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                                                                                                                                                                                  <  2:1079632/238‑1 (MQ=255)
 gATCCATAACCGTTCCTTTGCAATTCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                                                                                                                                                                                  <  1:239489/239‑1 (MQ=255)
 gATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                                                                                                                                                                                  <  1:302078/239‑1 (MQ=255)
  aTCCATAACCGTTCCTTTGCAATACCGCTATTGTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACAAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCAt                                                                                                                                                                                                                                                 <  1:401126/239‑1 (MQ=255)
   tCCATAACCGTTCCTTTGCAATGCCGCTTTTTTCTTGGCAACAGATGTTTCGACTATAGGGGGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATg                                                                                                                                                                                                                                                <  2:152663/239‑1 (MQ=255)
   tCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGTCGTTTAGGAAAATTTGAAATTGTTTTAGGGCAGCGGCAGGTTCAAACTATGGGGGTAAAAAAAAATGCagg                                                                                                                                                                                                                                                >  2:401126/1‑237 (MQ=255)
   tCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATg                                                                                                                                                                                                                                                >  1:131087/1‑239 (MQ=255)
    ccATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATgg                                                                                                                                                                                                                                               <  2:131087/239‑1 (MQ=255)
       tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGGCTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                                                                                                                                                                                  <  2:299029/233‑1 (MQ=255)
       tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGTGAGCGTAAGAGAACAAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCCGTTTCATACTATGGCGGGAAAAAAAATTTTATggaaa                                                                                                                                                                                                                                            >  2:434135/1‑236 (MQ=255)
       tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTTCa                                                                                                                                                                                                                                                  >  1:299029/1‑233 (MQ=255)
             ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAGATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGGAAAAAAAATTTTCAGGGATTTTAAg                                                                                                                                                                                                                                       >  2:407566/1‑239 (MQ=255)
             ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGGAAAAAAAATTTTTAGGGTATTTaga                                                                                                                                                                                                                                       >  2:580272/1‑237 (MQ=255)
              tCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGGGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGAAGTGGCGGTTTAATATTATCGGGGTTAAAAAAATTGTATggggttttagaa                                                                                                                                                                                                                                     >  2:635290/1‑229 (MQ=255)
                   tGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGGAAAAACAATTTCAGGGAAttttgggggtc                                                                                                                                                                                                                                  >  1:323547/1‑229 (MQ=255)
                            gCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAACAAAATTTTCATGGGTTTTAAGGACCCACCTTGtttt                                                                                                                                                                                                                      >  1:351503/1‑238 (MQ=255)
                             cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAATTGCATGGGATTTAATGTGTCACTATGTGTTct                                                                                                                                                                                                                    >  1:290206/1‑238 (MQ=255)
                                   tcctGCCAATAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTg                                                                                                                                                                                                             <  1:407566/236‑1 (MQ=255)
                                   tCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTg                                                                                                                                                                                                             <  1:770629/239‑1 (MQ=255)
                                       gCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                                                                                                                                                                                  >  1:280623/1‑201 (MQ=255)
                                       gCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                                                                                                                                                                                  <  2:280623/201‑1 (MQ=255)
                                          atgaagtGTTTCGACTTGAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCAcc                                                                                                                                                                                                      <  1:717300/233‑1 (MQ=255)
                                          aTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAACTTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGGAAAAAAAATTGCATGGGATTTAAGGACTCCCCTTGTTTTTCGCATTGGGCCAc                                                                                                                                                                                                       >  1:151044/1‑238 (MQ=255)
                                                   ttCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTg                                                                                                                                                                                             <  1:531395/239‑1 (MQ=255)
                                                               gAGGGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCt                                                                                                                                                                                   <  1:434135/237‑1 (MQ=255)
                                                                  cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTTAAAAAAATTTCATGGGATTTAAGGACTCACTATGTTTTTCGCATTTGGCCACCCCACCCCTCCCCCGGGGCCCgcgg                                                                                                                                                                              >  2:214100/1‑236 (MQ=255)
                                                                  cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAATTGCATGGGATTTTAGGACTTCCTATGTTTTCCGCATTGCGCCCCCGTACCCCTGCCCTGgccgccggcg                                                                                                                                                                              >  1:204685/1‑231 (MQ=255)
                                                                  cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGGAAAAAAATTTTCAATGGTTTTAAAGAGACACTATGGTTTCCGCATTTTGCCACCGTGCCCCTGCTCTGGCGCTCggag                                                                                                                                                                              >  1:1036374/1‑237 (MQ=255)
                                                                   gTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTTCATGGGTTTTTAGGAACAACTATTTTTTCCCCACTTGCCCCCCCTTCCGGTTCCCTGCCGCTctctgt                                                                                                                                                                             >  1:366388/1‑234 (MQ=255)
                                                                    tAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGt                                                                                                                                                                             <  1:739030/238‑1 (MQ=255)
                                                                    tAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGt                                                                                                                                                                             <  1:279952/238‑1 (MQ=255)
                                                                    tAAGAGAACAAATGAAATTCCCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGt                                                                                                                                                                             <  1:467878/238‑1 (MQ=255)
                                                                            cGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTTTATGGGATTTAAGGACTCACTATGTTTTTCGCATTTGCCCACCCTACCCCTGCCCTGGCGCCCGGGGTTGGCttttt                                                                                                                                                                   >  1:297289/1‑237 (MQ=255)
                                                                            cGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTTCATGGTATTTAAGGACTCACTATGTTTTTCGCCTTGCGACACCGTACCGCTGCCCTGGGGCCCGGCGTATGGttttt                                                                                                                                                                   >  1:1061269/1‑237 (MQ=255)
                                                                            cGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTTCCCATTTCGCCACCGTACCCCTGCCCTGGGGCCCGCCGGTTGCTTTAt                                                                                                                                                                   >  1:1016333/1‑239 (MQ=255)
                                                                                       aCCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAATTGCATGGGATTTAAGAGCCCACTTTGTTTTCCGCATTGCGCCACCCTACCGCTGCACTTGCGGCCTGCGTCTGGTTTGTTTTCCCCCTa                                                                                                                                                         >  1:808818/1‑238 (MQ=255)
                                                                                        ccAATTAAAATTGTTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTGAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTaca                                                                                                                                                       <  2:808818/239‑1 (MQ=255)
                                                                                         cAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTacac                                                                                                                                                      <  2:263377/239‑1 (MQ=255)
                                                                                            ttAGAATGAGCAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGAAAGATTTGAAATTGTTTTACAGCAGCGGCAGTTTCATACTATGGCGGGTCAAAAACTTTGATGGTCTTTTAGGAACCACCTATTTTTGCGACAGTGGCCCACCAACCGCCGCCCCTGCGCCCCGCGTGGTGCTTTATTccccccaccccac                                                                                                                                                     >  2:597701/1‑231 (MQ=255)
                                                                                                 aTGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCg                                                                                                                                                <  1:706719/237‑1 (MQ=255)
                                                                                                      tagtTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAATTTCATGGTATTTCAGGAGTCCATATGTTTTCTGCAATGCGGCCCCGTAACGCTTCCCTGGGCCCCGGCCTTGTCTTTATTTTCCCCCTCCAccccccccccccac                                                                                                                                         >  2:30152/1‑229 (MQ=255)
                                                                                                         ttCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGGATTTAAGGACTCCCTATGTTTTTCGCATTTGGCCACCGTACCGCTGCCCTGGGGCTCGGGGTATGCTTTTTTTTCCCCCTAAACGCCTCGGCAccccaa                                                                                                                                       >  1:232661/1‑234 (MQ=255)
                                                                                                          tCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATg                                                                                                                                                                                                                                                <  1:1002367/136‑1 (MQ=255)
                                                                                                                             ttGGCAAAGTTAAATTCAAAAAGTTCTTAAGGCCCCGGATTTCGGGCTTTTAGGAAAATTTGAAATTGTTTTAGCGGGACGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTCCCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACt                                                                                                                  <  2:522653/239‑1 (MQ=255)
                                                                                                                             ttGGCAAAGCTAATATCAAAAAGTGCTTAAGGCCCCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACt                                                                                                                  <  1:233499/239‑1 (MQ=255)
                                                                                                                              tGGAAAAGGTAATATCAAAAAGTGCTTAAGGCACCGGATTCCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTc                                                                                                                 <  1:214100/239‑1 (MQ=255)
                                                                                                                                   aaGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCa                                                                                                           <  1:580272/239‑1 (MQ=255)
                                                                                                                                      cTAAAAGCAAAAAGGGCGTAAGGCCACGGATGTCGGGCGTTAAGGAAGATTTGAAATTGTTTTAGCGGAGCGGCAGTTTCATACTAGGGCGGTAAAAAAATTTGCATGGTATTGAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTACTAATACTc                                                                                                                 <  1:370393/231‑1 (MQ=255)
                                                                                                                                        aaTATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACa                                                                                                      <  1:380493/239‑1 (MQ=255)
                                                                                                                                          tatCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTTCATGGTATTTAAGGACTCACTATGTTTTTCGCATTTCGCCACCCTACCGCTGCCCTGGCGCTCGGCGGATGCTTTTTTTTTCCCGTACACGCCTCGTCACCCAAACCTGGCGATTTTGCTAATACTCAGGCCCGCCata                                                                                                    >  1:849853/1‑239 (MQ=255)
                                                                                                                                                    tGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTa                                                                                                                                         <  1:200562/193‑1 (MQ=255)
                                                                                                                                                              aCCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTTCATGGGATTTAAGGGATCACTATTTTTTTCGCCTTTCGCCACCGTACCCCTTGCCCGGCCCTCGGCCTATGCTTTTTTTTCCCCCCACCCCCCCCCTCCACTCAACCCGGCGGTTTTTGTAAATCACAGGCAACGACAATTTTCCCCTTTTTTTCCggg                                                                                 >  2:950173/1‑239 (MQ=255)
                                                                                                                                                                cGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCCAAACCCGGCGATTTTGCTAATTCTCAGGCACGACATATTTCCACTTTCTTTCCTGGAc                                                                               >  1:82389/1‑238 (MQ=255)
                                                                                                                                                                      tCGGGCTTTTAGGAAGATTTGAAATTGGTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTGAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGAc                                                                        <  1:1064868/239‑1 (MQ=255)
                                                                                                                                                                              ttAGGAAGGTTTGAATTTGTTTTTCTCCAGCGGAAGTTTCATTCTAGGCCAGTAAAAAAATTTTCCTCGTATCTAAGGACTCACTACGTTTTCCCCATGGCCCAACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCTTACACGCCTCGTCTCCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTc                                                                <  2:150675/239‑1 (MQ=255)
                                                                                                                                                                              ttAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTc                                                                <  1:558572/239‑1 (MQ=255)
                                                                                                                                                                                 ggAAGTTTTGAAATTGTTTTAGCGCAGCGGGAGTTTCATACTAGGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTg                                                             <  1:255390/239‑1 (MQ=255)
                                                                                                                                                                                    agttttGAATTTGTTTTAGCGCATGGGCCGTTTCATGTTATGGCGTAAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCa                                                           <  2:1005384/235‑1 (MQ=255)
                                                                                                                                                                                    aGATTTTAAATTGTTTGAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATGTGCATGTTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCa                                                           <  1:491863/238‑1 (MQ=255)
                                                                                                                                                                                                 ttttAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCa                                                                                                           <  1:727776/177‑1 (MQ=255)
                                                                                                                                                                                                                   ttctttttttGGCGGTAAAAAAATTTGCATGGTATTTAAGGACCCACTATGTTTTCCGCATTGCCCCCCCGTCCCGCTGCCCTGGCGCCCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCAGAAATGTTATCTGCCGATTATATTCGCCa                              <  2:849853/227‑1 (MQ=255)
                                                                                                                                                                                                                              gCGGTAAAAAAATTTGCATGGTATTTAAGGACTGACTATGGTTTCCGAATTGCGCCCCCGATCCGCTGCCCTGGCGCTAGGCGTAAGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTTCTTTCTTTCCGGGACGCATGACCGGAACTCCTGCAGAAATGTTATCTGCCGATTATATTCGCCAACAGTTTCAGCaaa                <  1:94520/239‑1 (MQ=255)
                                                                                                                                                                                                                               cGGTAAAAAAATTTGCATGGTATTTAAGGACCCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCAGAAATGTTATCTGCCGATTATATTCGCCAACAGTTTCAGCaaa                <  1:252318/238‑1 (MQ=255)
                                                                                                                                                                                                                                 gggtaaaaaaTTTTCATGGTTTTTAAGGGCTCCACTTTTTTTTCCGCTTGGCGCCACCGTCCCGTTCCCCTGCCGTCCGGGTTATGCTTTTTTCCCCCGGTACACGCCTCGTCACCGAACCCTGGCGATTTTCATAATACACAGGCACGAAATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACCCCTGCAGAAATGTTATCTCCCGATTATATCCGCCAACAGTTGCAGaaaa                <  2:351503/234‑1 (MQ=255)
                                                                                                                                                                                                                                        aattttGCCTTGTATTTAAGGACTCCCTATGTTTTCCGCATTGCCCCCCCGTACCTGTGCCCGGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCAGAAATGTTATCTGCCGATTATATTCGCCAACAGTTTCAGCAAATGGGTTATc       <  2:290206/236‑1 (MQ=255)
                                                                                                                                                                                                                                         attttGCATGGTTTTTAGGGTCTCCCTTTGTTTTCCGCATTGCCCCCCCGTACCGCTGCCCTGGCGCTCGGCGTATGTTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCAGAAATGTTATCTGCCGATTATATTCGCCAACAGTTTCAGCAAATGGGTTATCg      <  2:232661/237‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCAGAAATGTTATCTGCCGATTATATTCGCCAACAGTTTCCGCAAATGGGGTaaccaagg  >  2:897739/1‑232 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
TGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTCGTCACCTAAACCTGGCGATTTTGCTAATACTCAGGCACGACATATTGCTACTTTCTTTCCGGGACGCATGACCGGAACTCCTGCAGAAATGTTATCTGCCGATTATATTCGCCAACAGTTTCAGCAAATGGGTTATCGCAGT  >  NC_000913/2876322‑2876796

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