Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,623 C→A 7.8% V640L (GTA→TTA)  yfjK ← radiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,6230CA 7.8% 985.4 / 52.3 284V640L (GTA→TTA) yfjKradiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base C (122/139);  new base A (13/9);  total (135/149)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 288 total aligned reads displayed.

ATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTCTGG  >  NC_000913/2761350‑2761895
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                
aTCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACa                                                                                                                                                                                                                                                                                 >  4:322011/1‑275 (MQ=255)
                                                         tCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTTCATCCTCCCTCAAGATGAGCACCTTAGGTAGGGTATAGCCg                                                                                                                                                                                                                        >  1:15713/1‑275 (MQ=255)
                                                                      aaTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAATCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtag                                                                                                                                                                                                                                  >  2:82763/1‑252 (MQ=255)
                                                                       aTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACAACTGGTGGTAGATGGAAACTTGCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGGGTTTAACTAAATCCTCCATCAACCTGCACACCTTAGGGCGGGTATGACCGAAGGTATTGGGTAt                                                                                                                                                                                                          >  4:245669/1‑275 (MQ=255)
                                                                           ccTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                                                                                                                                                  >  2:40885/1‑199 (MQ=255)
                                                                                                      aCTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACAACGCTTTACa                                                                                                                                                                                                                                                                                 >  4:61814/1‑173 (MQ=255)
                                                                                                                  cAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                                                                                                                                                  >  2:309427/1‑160 (MQ=255)
                                                                                                                             gAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                                                                                                                                                  >  5:961270/1‑149 (MQ=255)
                                                                                                                                aaTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGc                                                                                                                                                                                                                                                                               >  8:1064945/1‑149 (MQ=255)
                                                                                                                                  ttCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                                                                                                                                                             >  8:753820/1‑149 (MQ=255)
                                                                                                                                  ttCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                                                                                                                                                             <  5:743193/149‑1 (MQ=255)
                                                                                                                                   tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                            <  7:1096234/149‑1 (MQ=255)
                                                                                                                                   tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                            <  5:1123205/149‑1 (MQ=255)
                                                                                                                                   tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                            <  7:881964/149‑1 (MQ=255)
                                                                                                                                   tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtg                                                                                                                                                                                                                                                                            <  7:605610/149‑1 (MQ=255)
                                                                                                                                     cAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggt                                                                                                                                                                                                                                                                          <  5:811178/149‑1 (MQ=255)
                                                                                                                                     cAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                            <  7:52611/147‑1 (MQ=255)
                                                                                                                                      aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                                                                                                                                                         <  6:797254/149‑1 (MQ=255)
                                                                                                                                      aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                                                                                                                                                         >  7:769831/1‑149 (MQ=255)
                                                                                                                                      aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                                                                                                                                                         >  6:402405/1‑149 (MQ=255)
                                                                                                                                      aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggta                                                                                                                                                                                                                                                                         <  5:1190484/149‑2 (MQ=255)
                                                                                                                                        cAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGtt                                                                                                                                                                                                                                                                       <  7:1419874/149‑1 (MQ=255)
                                                                                                                                         aaGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGttt                                                                                                                                                                                                                                                                      <  7:282386/149‑1 (MQ=255)
                                                                                                                                           gTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                                                                                                                                                                    <  8:869325/149‑1 (MQ=255)
                                                                                                                                            tCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAAc                                                                                                                                                                                                                                                                   >  7:311525/1‑149 (MQ=255)
                                                                                                                                                  aGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                                                                                                                                                             <  6:439218/149‑1 (MQ=255)
                                                                                                                                                   gTCGCTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcc                                                                                                                                                                                                                                                            >  8:954290/1‑149 (MQ=255)
                                                                                                                                                    tCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcct                                                                                                                                                                                                                                                           >  6:1012404/1‑149 (MQ=255)
                                                                                                                                                      ggTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGGGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCAATGTCGTGCATAATCGAtt                                                                                                                           >  4:233769/1‑274 (MQ=255)
                                                                                                                                                        tgtAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTGCAGCGTGGTGTTTAACTACATCCTCCAt                                                                                                                                                                                                                                                       >  8:935956/1‑149 (MQ=255)
                                                                                                                                                         gtAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATc                                                                                                                                                                                                                                                      <  8:1189796/149‑1 (MQ=255)
                                                                                                                                                          tAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtgg                                                                                                                                                                                                                                                                           <  8:183986/127‑1 (MQ=255)
                                                                                                                                                           aaCGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCGGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCaa                                                                                                                                                                                                                                                    >  6:1024573/1‑149 (MQ=255)
                                                                                                                                                              gTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                            <  5:1084020/122‑1 (MQ=255)
                                                                                                                                                               tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTg                                                                                                                                                                                                                                                >  5:776350/1‑149 (MQ=255)
                                                                                                                                                               tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTg                                                                                                                                                                                                                                                >  6:901392/1‑149 (MQ=255)
                                                                                                                                                                    aTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCAc                                                                                                                                                                                                                                           <  6:508594/149‑1 (MQ=255)
                                                                                                                                                                     tGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                                                                                                                                                               >  5:1359503/1‑144 (MQ=255)
                                                                                                                                                                      gggAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCt                                                                                                                                                                                                                                         >  6:1128028/1‑149 (MQ=255)
                                                                                                                                                                      gggAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCt                                                                                                                                                                                                                                         >  5:249941/1‑149 (MQ=255)
                                                                                                                                                                         aaTGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtag                                                                                                                                                                                                                                      <  6:391653/149‑1 (MQ=255)
                                                                                                                                                                              cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                                 <  8:405684/149‑1 (MQ=255)
                                                                                                                                                                               aaTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                                                                                                                                                                <  6:847198/149‑1 (MQ=255)
                                                                                                                                                                               aaTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                                                                                                                                                                >  7:199783/1‑149 (MQ=255)
                                                                                                                                                                                 tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTa                                                                                                                                                                                                                                                                 >  1:824/1‑114 (MQ=255)
                                                                                                                                                                                     ttctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                                                                                                                                                             <  3:824/114‑1 (MQ=255)
                                                                                                                                                                                          aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                                                                                                                                                     <  5:314169/149‑1 (MQ=255)
                                                                                                                                                                                          aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                                                                                                                                                     <  5:402405/149‑1 (MQ=255)
                                                                                                                                                                                              aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAt                                                                                                                                                                                                                 <  6:1400139/149‑1 (MQ=255)
                                                                                                                                                                                               ttaacccttggtggAGATGGAAACTTTCGAACGCCAATTTGACATGTTGGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAg                                                                                 <  4:9017/267‑1 (MQ=255)
                                                                                                                                                                                               ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAAc                                                                                                  >  2:258847/1‑259 (MQ=255)
                                                                                                                                                                                                tGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                                                                                                                                                      >  7:97233/1‑142 (MQ=255)
                                                                                                                                                                                                 gACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                            <  6:301878/87‑1 (MQ=255)
                                                                                                                                                                                                  acacCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGt                                                                                                                                                                                                             <  8:973619/149‑1 (MQ=255)
                                                                                                                                                                                                   cacCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAg                                                                                                                <  4:308357/241‑1 (MQ=255)
                                                                                                                                                                                                    acCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATGCCACGCTTTACAGAGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTa                                                                                                                                                                                                           >  8:695165/1‑149 (MQ=255)
                                                                                                                                                                                                    acCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                                                                                                                                                     >  6:1154877/1‑139 (MQ=255)
                                                                                                                                                                                                    acCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTa                                                                                                                                                                                                           <  7:505355/149‑1 (MQ=255)
                                                                                                                                                                                                     ccTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTAt                                                                                                                                                                                                          <  7:1168095/149‑1 (MQ=255)
                                                                                                                                                                                                     ccTGGCGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTAt                                                                                                                                                                                                          <  5:1012404/149‑1 (MQ=255)
                                                                                                                                                                                                      cTGGTTGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGAGTATACCCCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATg                                                                                                                                                                                                         >  6:1246104/1‑149 (MQ=255)
                                                                                                                                                                                                       tggtggTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                                                                                                                                                                        <  7:168789/149‑1 (MQ=255)
                                                                                                                                                                                                            gtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                                 <  8:383676/119‑1 (MQ=255)
                                                                                                                                                                                                            gtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                                 >  7:383676/1‑119 (MQ=255)
                                                                                                                                                                                                             tAGATGGAAACTTGCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGAGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGAACCTTAGGTAGGGTATAGCCGAAGGTAGTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAAGAACTTTCAGCACTTAAATAGCCGGCAAGTTTTg                                                                    <  2:248350/275‑1 (MQ=255)
                                                                                                                                                                                                               gATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGtt                                                                                                                                                                                                            <  8:523901/137‑1 (MQ=255)
                                                                                                                                                                                                               gATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCaaa                                                                                                                                                                                                <  6:965584/149‑1 (MQ=255)
                                                                                                                                                                                                                 tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCCGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATc                                                                                                                                                                                              >  5:703396/1‑149 (MQ=255)
                                                                                                                                                                                                                   gAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAt                                                                                                                                                                                            >  8:1235881/1‑149 (MQ=255)
                                                                                                                                                                                                                    aaaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                           <  6:249941/149‑1 (MQ=255)
                                                                                                                                                                                                                    aaaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                           >  6:1438804/1‑149 (MQ=255)
                                                                                                                                                                                                                    aaaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                           <  7:936516/149‑1 (MQ=255)
                                                                                                                                                                                                                      aCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTCAAATAGCCGGCAAGTTTTGCGATCATGa                                                           <  2:198836/275‑1 (MQ=255)
                                                                                                                                                                                                                        tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                           >  6:357060/1‑145 (MQ=255)
                                                                                                                                                                                                                          tCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAg                                                                                                                                                                                     <  7:724356/149‑1 (MQ=255)
                                                                                                                                                                                                                               cGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCaa                                                                                                                                                                                >  8:483278/1‑149 (MQ=255)
                                                                                                                                                                                                                                gCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAAt                                                                                                                                                                               >  8:1444271/1‑149 (MQ=255)
                                                                                                                                                                                                                                gCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAAt                                                                                                                                                                               <  8:1310080/149‑1 (MQ=255)
                                                                                                                                                                                                                                          cATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                                                                                                                                                     <  6:701238/149‑1 (MQ=255)
                                                                                                                                                                                                                                           aTGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                                                                                                                                                               >  5:683110/1‑106 (MQ=255)
                                                                                                                                                                                                                                             gtgtgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                                                                                                                                                     >  7:326324/1‑146 (MQ=255)
                                                                                                                                                                                                                                                    aTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                                                                                                                                                      <  7:741609/90‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTCGACTTTGCTGCATATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTccc                                                                                                                                                           >  7:104357/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         cTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCa                                                                                                                                                      >  8:27541/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            tgctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGaaa                                                                                                                                                                          >  8:622209/1‑126 (MQ=255)
                                                                                                                                                                                                                                                             gctgcGTATACTACTCTTTACAGCGTGGTGTTTAAATACAACCTCCAACAACCGGCGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATTATTGTTGAGGGCAATGGACAGAGTACGTTTCTCCCCTTCAGaca                                                                                                                                                  <  8:31589/149‑3 (MQ=255)
                                                                                                                                                                                                                                                             gctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                                 <  7:1209674/70‑1 (MQ=255)
                                                                                                                                                                                                                                                                gcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATcc                                                                                                                                               <  8:211291/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    acaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAgg                                                                                                                                                                                    >  5:1139003/3‑108 (MQ=255)
                                                                                                                                                                                                                                                                     caCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTGGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGt                                                                                                                                          >  5:605737/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                     caCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                                                                                                                                                     >  7:393816/2‑122 (MQ=255)
                                                                                                                                                                                                                                                                      aCCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTc                                                                                                                                         <  8:1409152/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTc                                                                                                                                         >  7:1350786/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       ccACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCg                                                                                                                                        >  6:620193/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt                                                                                                                                       >  7:490684/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg                                                                                                                                      >  7:13394/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           gCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGt                                                                                                                                                                                                             <  5:1285674/76‑1 (MQ=255)
                                                                                                                                                                                                                                                                             tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGaaa                                                                                                                                                                          >  5:364711/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                             tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATa                                                                                                                                  >  8:1098318/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                aCAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATc                                                                                                                               <  5:1406063/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                
ATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTCTGG  >  NC_000913/2761350‑2761895

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