Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,623 C→A 9.1% 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 9.1% 1281.6 / 87.0 374V640L (GTA→TTA) yfjKradiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base C (171/169);  new base A (14/20);  total (185/189)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-04

Only 100 of 379 total aligned reads displayed.

AATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTCTGG  >  NC_000913/2761353‑2761895
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                
actctcgtcTAATACCCTGCTCTCGCCTTATAAACATCTTCTCCACTTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGCTACTGGCTCAACTCATCCACCTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGGCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCg                                                                                                                                                                                                                                                                              <  4:403212/271‑1 (MQ=255)
       tcTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCCAAGGCATTGACACATGTGGGAAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATGGATTTTGTNTAGTTNNCCCCCCTTTACaagcttttttttta                                                                                                                                                                                                                                                                    >  3:19243/1‑265 (MQ=255)
                 aGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAgtg                                                                                                                                                                                                                                                                              >  1:211879/1‑256 (MQ=255)
                 aGCGCGCGCCTGATAAACATCTGATCAACGTACCAAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCAAAGGCATTGAAACCTGGTGGTAGATGGAAACTTTCGAAAGCCAAGTTGACATGTGTGTAATCGACTTTGCTGCGTATACCATGCTTTACAGCGTGGGGTTTAACTAcct                                                                                                                                                                                                                                                              >  4:146410/1‑272 (MQ=255)
                                aaCATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                                                                                                                                                               <  2:100933/274‑1 (MQ=255)
                                        aTCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGGGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCTTATACCACGCTTTACAGCGTGGTGTTTAACTATATCCTCCATCAAGCTGAGCACCTtag                                                                                                                                                                                                                                      >  4:57010/1‑275 (MQ=255)
                                                                                                  aaCTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGa                                                                                                                                                                            <  2:647824/275‑1 (MQ=255)
                                                                                                    cTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                                                                                                                                                  >  4:88735/1‑171 (MQ=255)
                                                                                                        tCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCg                                                                                                                                                                                                                                                                              <  2:88735/171‑1 (MQ=255)
                                                                                                                     cgtcgGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCg                                                                                                                                                                                                                                                                              >  1:148323/1‑158 (MQ=255)
                                                                                                                     cgtcgGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCg                                                                                                                                                                                                                                                                              <  3:148323/158‑1 (MQ=255)
                                                                                                                                tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCGAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa                                                                                                                                              <  2:146412/275‑1 (MQ=255)
                                                                                                                                  cAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggt                                                                                                                                                                                                                                                                          <  8:1531346/149‑1 (MQ=255)
                                                                                                                                         tCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCtt                                                                                                                                                                                                    >  2:340111/1‑212 (MQ=255)
                                                                                                                                                gTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcc                                                                                                                                                                                                                                                            >  6:454453/1‑149 (MQ=255)
                                                                                                                                                     tgtAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCTTTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGACATCCACTTCCTTGCATAACCCATCTg                                                                                                                         >  4:298329/1‑275 (MQ=255)
                                                                                                                                                       tAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGtt                                                                                                                                                                                                                                                                       <  2:721921/131‑1 (MQ=255)
                                                                                                                                                       tAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                                                                                                                                                             >  8:314628/1‑141 (MQ=255)
                                                                                                                                                         aCGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAg                                                                                                                                                                                                                                                   >  8:1049429/1‑149 (MQ=255)
                                                                                                                                                          cGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGc                                                                                                                                                                                                                                                  <  5:1383476/149‑1 (MQ=255)
                                                                                                                                                             cTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                                                                                                                                                               <  6:163069/149‑1 (MQ=255)
                                                                                                                                                               gAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGaa                                                                                                                                                                           >  2:552351/1‑215 (MQ=255)
                                                                                                                                                                aaTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                            >  7:499135/1‑117 (MQ=255)
                                                                                                                                                                aaTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                                                                                                                                                            <  8:499135/117‑1 (MQ=255)
                                                                                                                                                                 aTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCAc                                                                                                                                                                                                                                           >  6:769125/1‑149 (MQ=255)
                                                                                                                                                                    ggAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCtt                                                                                                                                                                                                                                        <  7:982193/149‑1 (MQ=255)
                                                                                                                                                                         gCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                                                                                                                                                         <  4:191910/111‑1 (MQ=255)
                                                                                                                                                                         gCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggta                                                                                                                                                                                                                                   >  7:1585712/1‑149 (MQ=255)
                                                                                                                                                                            aaTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                                                                                                                                                                >  6:1124604/1‑149 (MQ=255)
                                                                                                                                                                             aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                                 >  7:279578/1‑147 (MQ=255)
                                                                                                                                                                             aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                                 <  8:279578/147‑1 (MQ=255)
                                                                                                                                                                             aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGt                                                                                                                                                                                                                               >  5:1652976/1‑149 (MQ=255)
                                                                                                                                                                              tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                                                                                                                                                       >  5:1175020/1‑140 (MQ=255)
                                                                                                                                                                              tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                                                                                                                                                       <  6:1175020/140‑1 (MQ=255)
                                                                                                                                                                                   tctAAGGCATTGACCCCTGGTGGTAGATGTGAACTTTCCAACGCCATTTTGACAGGTGTTTAATCGACTTTTTTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAg                                                                                            <  4:329527/274‑1 (MQ=255)
                                                                                                                                                                                    ctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGGATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCg                                                                                                                                                                                                                        >  6:1552849/1‑149 (MQ=255)
                                                                                                                                                                                     tAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCa                                                                                                                                                                                                                                            >  6:240153/1‑128 (MQ=255)
                                                                                                                                                                                     tAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGa                                                                                                                                                                                                                       >  5:1641394/1‑149 (MQ=255)
                                                                                                                                                                                      aaGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                                                                                                                                                      <  8:1269006/149‑1 (MQ=255)
                                                                                                                                                                                        ggCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAgg                                                                                                                                                                                                                    >  8:160235/1‑149 (MQ=255)
                                                                                                                                                                                           aTTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAt                                                                                                                                                                                                                 >  5:471040/1‑149 (MQ=255)
                                                                                                                                                                                            ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAtt                                                                                                                                                                                                                >  6:57611/1‑149 (MQ=255)
                                                                                                                                                                                            ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTTCATAATCGATCTGATTTTTCAGATCTTCCGAATAACTTTCAGCACTTAAAt                                                                                   >  2:409490/1‑274 (MQ=255)
                                                                                                                                                                                             tGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                                                                                                                                                               <  8:1143280/149‑1 (MQ=255)
                                                                                                                                                                                             tGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTACTg                                                                                                                                                                                                               <  5:42157/149‑1 (MQ=255)
                                                                                                                                                                                             tGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGACGAAGGTATTg                                                                                                                                                                                                               <  5:1274520/149‑1 (MQ=255)
                                                                                                                                                                                              gACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                                                                                                                                                      <  8:6664/141‑1 (MQ=255)
                                                                                                                                                                                              gACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                                                                                                                                                      >  7:6664/1‑141 (MQ=255)
                                                                                                                                                                                               acacCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGcc                                                                               <  2:197343/275‑1 (MQ=255)
                                                                                                                                                                                                   cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAg                                                                                                                                                                                                                                              >  4:463733/1‑112 (MQ=255)
                                                                                                                                                                                                   cTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATg                                                                                                                                                                                                         >  7:169156/1‑149 (MQ=255)
                                                                                                                                                                                                             aTGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAAt                                                                                                                                                                                               >  8:139430/1‑149 (MQ=255)
                                                                                                                                                                                                              tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCaaa                                                                                                                                                                                                <  7:370283/147‑1 (MQ=255)
                                                                                                                                                                                                               ggAATCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATATCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                                                                                                                                                             >  8:910941/1‑149 (MQ=255)
                                                                                                                                                                                                               ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                                                                                                                                                                    <  7:330884/78‑1 (MQ=255)
                                                                                                                                                                                                               ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                                                                                                                                                             >  6:714717/1‑149 (MQ=255)
                                                                                                                                                                                                                  aaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAATGTATTGGTTATGAGCTTCAAATCAttg                                                                                                                                                                                          >  7:1014302/1‑149 (MQ=255)
                                                                                                                                                                                                                  aaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAttg                                                                                                                                                                                          >  5:1587023/1‑149 (MQ=255)
                                                                                                                                                                                                                    cTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                                                                                                                                                             <  2:604756/144‑1 (MQ=255)
                                                                                                                                                                                                                     tttcaaactaaaaTTTGAAATTTGTTNNCTCGTNTTTGCTGCGTAAACCACGCTTTACAGCGTGGGGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCg                                                     <  3:7422/268‑1 (MQ=255)
                                                                                                                                                                                                                         gAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAggg                                                                                                                                                                                   >  6:894640/1‑149 (MQ=255)
                                                                                                                                                                                                                            cGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCaa                                                                                                                                                                                <  8:169156/149‑1 (MQ=255)
                                                                                                                                                                                                                              ccAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGtata                                                                                                                                                                                                                            <  5:1338829/103‑1 (MQ=255)
                                                                                                                                                                                                                              ccAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATg                                                                                                                                                                              >  7:142941/1‑149 (MQ=255)
                                                                                                                                                                                                                                             tgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGaa                                                                                                                                                                           >  4:627043/1‑137 (MQ=255)
                                                                                                                                                                                                                                             tgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa                                                                                                                                              >  1:63387/1‑166 (MQ=255)
                                                                                                                                                                                                                                               tAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCGAGTTTTGCGATCATGACATCCGt                                                    <  2:288011/254‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGaaa                                                                                                                                                                          >  8:776109/1‑134 (MQ=255)
                                                                                                                                                                                                                                                 aTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTccc                                                                                                                                                           >  7:528341/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 aTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaa                                                                                                                                                   <  4:713250/157‑1 (MQ=255)
                                                                                                                                                                                                                                                  tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtagg                                                                                                                                                                                                                                     <  6:347522/74‑1 (MQ=255)
                                                                                                                                                                                                                                                  tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtagg                                                                                                                                                                                                                                     >  5:347522/1‑74 (MQ=255)
                                                                                                                                                                                                                                                  tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTcccc                                                                                                                                                          <  7:898131/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTcccc                                                                                                                                                          <  5:1124604/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    gACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCtt                                                                                                                                                        <  8:350359/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCa                                                                                                                                                      <  7:417032/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       tttGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGtt                                                                                                                                                                                                            <  7:1355604/94‑1 (MQ=255)
                                                                                                                                                                                                                                                       tttGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCt                                                                                                                                                         <  5:195392/145‑1 (MQ=255)
                                                                                                                                                                                                                                                         tgctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTTTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTAAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCAACTGAAc                     >  4:232599/1‑275 (MQ=255)
                                                                                                                                                                                                                                                          gctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaaa                                                                                                                                                  <  8:8857/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           ctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCGCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa                                                                                                                                              <  2:428175/152‑1 (MQ=255)
                                                                                                                                                                                                                                                            tgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATc                                                                                                                                                >  7:1282824/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                tataCCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCCCTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAAc          >  1:107723/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                 acaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTg                                                                                                                                           >  5:241224/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                  caCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                           >  5:992161/2‑100 (MQ=255)
                                                                                                                                                                                                                                                                  caCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAagag                                                                                                                                                                       <  8:545969/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTc                                                                                                                                         <  5:800593/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt                                                                                                                                       >  7:785819/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt                                                                                                                                       <  8:1337378/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt                                                                                                                                       >  5:722368/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt                                                                                                                                       <  8:427925/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt                                                                                                                                       >  6:1280641/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTc                                                                                                                                                                                                   <  6:1655702/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg                                                                                                                                      <  7:566925/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aCGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                                                                                                                                                     <  5:865786/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tttCCAGCGTGGTGTTTAACTACAGCCTCCATCAAGCTGAGCACCTTAGTGAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                                                                                                                                                      <  3:36295/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                                                                                                                                                           <  6:1337518/92‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGc                                                                                                                                                                                                      <  6:1436625/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                             aCAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATc                                                                                                                               <  6:834260/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cAGCGTGGTGTTTAACTACATCCTGCATCAAGCTGAGCACGTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGACCTTCAAATCATGGTTGAGGGCAATGGAAAGAGTAAGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAAtag                                                                                                                              >  7:1158308/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                
AATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTCTGG  >  NC_000913/2761353‑2761895

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