Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,161 A→C 6.1% D87A (GAC→GCC) ‡ ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913569,1610AC 6.1% 588.7 / 25.0 165D87A (GAC→GCC) ‡ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base A (79/76);  new base C (5/5);  total (84/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-02

Only 100 of 165 total aligned reads displayed.

ACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGA  >  NC_000913/569025‑569429
                                                                                                                                        |                                                                                                                                                                                                                                                                             
aCCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                                                                                                                                                                                                       >  4:75882/1‑193 (MQ=255)
    tAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                                                                                                                                                                                                   <  2:75883/193‑1 (MQ=255)
        atatGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACa                                                                                                                                                                                                                                                                             <  4:567376/131‑1 (MQ=255)
          atGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                                                                                                                                                                                                       >  2:479297/1‑183 (MQ=255)
            gAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                                                                                                                                                                                                       >  4:332259/1‑149 (MQ=255)
            gAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTa                                                                                                                                                                                                              >  4:511140/1‑190 (MQ=255)
            gAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc                                                                                                                                                                                       >  2:210300/1‑213 (MQ=255)
            gAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCata                                                                                                                                                >  2:255262/1‑251 (MQ=255)
             aaGAGTTAGGATTAAGTGCATTCAAAGGAAAGTATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCgg                                                                                                                                                                                                                                                       >  5:320755/1‑148 (MQ=255)
              aGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                                                                                                                                                                                                   <  4:479296/183‑1 (MQ=255)
                aTTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                                                                                                                                                          <  2:511141/190‑1 (MQ=255)
                aTTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                                                                                                                                                   <  4:210299/213‑1 (MQ=255)
                  ttAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                                                                                                                                                                                                     <  5:927658/129‑1 (MQ=255)
                  ttAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                                                                                                                                                                                                     >  6:927658/1‑129 (MQ=255)
                  ttAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc                                                                                                                                                                                       >  1:117085/1‑207 (MQ=255)
                   tAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACttt                                                                                                                                                                                                                                                                       <  8:254121/126‑1 (MQ=255)
                   tAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTCTTGATGTTATAACTCTTTGCGCCAATACATTCTATCATATTTACTCTTTGAATGGCCCATATTCATTAATACAAGCCATACTtatcactcaa                                                                                                                 >  4:138050/1‑268 (MQ=255)
                      taTTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                                                                                                                                                   <  3:117085/206‑1 (MQ=255)
                      gATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                                                                                                                                                           <  4:219684/199‑1 (MQ=255)
                       aTTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACtt                                                                                                                                                                                                                                                                        <  6:247683/121‑1 (MQ=255)
                          aaGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                                                                                                                                                                                                          <  7:10844/148‑1 (MQ=255)
                          aaGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTataa                                                                                                                                                                         >  2:73598/1‑213 (MQ=255)
                                aTTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAGTTTTTAGATGCTATAGAGCATGGTTATATATTTCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                                                                                                                                                                                                       >  1:25503/1‑161 (MQ=255)
                                aTTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                                                                                                                                                   >  2:245720/1‑197 (MQ=255)
                                  tCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACt                                                                                                                                                                                                                                                                         >  5:412931/1‑109 (MQ=255)
                                   cAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                                                                                                                                                                                                     <  5:296669/112‑1 (MQ=255)
                                     ggacagaGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCGATAATTTTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAg                                                                                                <  2:487793/268‑1 (MQ=255)
                                     aaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                                                                                                                                                   >  7:290225/1‑112 (MQ=255)
                                             cATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAg                                                                                                                                                                                                                                                                            >  7:1113035/1‑95 (MQ=255)
                                             cATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAg                                                                                                                                                                                                                                                                            <  8:1113035/95‑1 (MQ=255)
                                               tGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                                                                                                                                                                                               <  5:832341/122‑1 (MQ=255)
                                                gCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc                                                                                                                                                                                       >  2:552389/1‑177 (MQ=255)
                                                    aaTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                                                                                                                                                   <  4:552388/177‑1 (MQ=255)
                                                      tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGAAAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGATATAACTCTTTGCGACAATCAAGTCTATAATATTGACTCTTTGAAGGAGCCTTTTTCATTAATAAAAGCCATTATTATAGCCAAAAGGGCAAATGAAAAAAGCGAATaaaatttcgttcc                                                                             >  4:487792/1‑265 (MQ=255)
                                                                tttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGa                                                                                                                                                                                                                       >  2:39975/1‑129 (MQ=255)
                                                                    ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                                                                                                                                                                                                   <  4:39974/129‑1 (MQ=255)
                                                                     gAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAgaga                                                                                                                                                                                                                                                                >  6:401741/1‑83 (MQ=255)
                                                                                    ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAtt                                                                                                                                                                                                            <  8:156089/120‑1 (MQ=255)
                                                                                     taGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                                                                                                                                                          <  5:650915/121‑1 (MQ=255)
                                                                                     taGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACaa                                                                                                                                                                                     >  3:66248/1‑142 (MQ=255)
                                                                                        agCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGtt                                                                                                                                                                                                      >  2:568917/1‑122 (MQ=255)
                                                                                         gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                                                                                                                                                          >  8:836923/1‑117 (MQ=255)
                                                                                         gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACa                                                                                                                                                                                 <  1:66248/142‑1 (MQ=255)
                                                                                          cATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAtt                                                                                                                                                                                                                               >  3:29043/1‑95 (MQ=255)
                                                                                            tGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATaa                                                                                                                                                                                                  <  4:568916/122‑1 (MQ=255)
                                                                                             ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                                                                                                                                                                                            >  2:552061/1‑127 (MQ=255)
                                                                                             ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGGCAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAACGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGAAAGAAAAAACGGCAGGATGCACt                                       >  1:91687/1‑275 (MQ=255)
                                                                                             ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGAc                                                                                                                                                                                       >  4:3805/1‑132 (MQ=255)
                                                                                             ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACaa                                                                                                                                                                                     >  4:547442/1‑134 (MQ=255)
                                                                                             ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACaa                                                                                                                                                                                     >  4:247461/1‑134 (MQ=255)
                                                                                             ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACANTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCAATAATTATAGCACAAAGGCCAAATGAAGAAAGCNAGATAAAGTCAAGTCGGGTTAAATCNAAAGGGAAAAAAAAAAAACAGGAAGCaat                                       >  3:97256/1‑273 (MQ=255)
                                                                                                tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                                                                                                                                                                                            >  7:1143514/1‑92 (MQ=255)
                                                                                                tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                                                                                                                                                                                                   <  7:587227/117‑1 (MQ=255)
                                                                                                tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtata                                                                                                                                                                                                   >  8:587227/1‑117 (MQ=255)
                                                                                                tatataTTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                                                                                                                                                   <  2:153711/133‑1 (MQ=255)
                                                                                                 atatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGa                                                                                                                                                                                        <  4:552060/127‑1 (MQ=255)
                                                                                                 atatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACa                                                                                                                                                                                 <  2:547443/134‑1 (MQ=255)
                                                                                                 atatatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATTATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAaagtcaag                                                                                 <  2:96065/230‑1 (MQ=255)
                                                                                                    tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTc                                                                                                                                                                                                                                                         >  6:586524/1‑59 (MQ=255)
                                                                                                    tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTc                                                                                                                                                                                                                                                         <  5:586524/59‑1 (MQ=255)
                                                                                                    tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa                                                                                                                                                                                      >  7:23335/1‑126 (MQ=255)
                                                                                                    tatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa                                                                                                                                                                                      <  8:23335/126‑1 (MQ=255)
                                                                                                     atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGa                                                                                                                                                                                                                                                    >  5:322789/1‑63 (MQ=255)
                                                                                                     atTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGa                                                                                                                                                                                                                                                    <  6:322789/63‑1 (MQ=255)
                                                                                                       tGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTg                                                                                                                                                                                                                              >  6:578455/1‑83 (MQ=255)
                                                                                                       tGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTAt                                                                                                                                                                                                             >  5:582146/1‑100 (MQ=255)
                                                                                                       tGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCttt                                                                                                                                                            >  6:691491/1‑148 (MQ=255)
                                                                                                         ccAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTa                            >  2:27168/1‑274 (MQ=255)
                                                                                                           aGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                                                                                                                                                   <  2:396571/122‑1 (MQ=255)
                                                                                                           aGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                                                                                                                                                   <  1:33396/122‑1 (MQ=255)
                                                                                                           aGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAAt                                                                                                                                                        <  4:69356/148‑1 (MQ=255)
                                                                                                              aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                                                                                                                                                             >  6:101984/1‑109 (MQ=255)
                                                                                                              aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                                                                                                                                                             <  5:21450/109‑1 (MQ=255)
                                                                                                              aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                                                                                                                                                             <  5:101984/109‑1 (MQ=255)
                                                                                                                ttaaaTTGTTAATTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATg                    <  4:27167/271‑1 (MQ=255)
                                                                                                                  aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                                                                                                                                                                                               >  8:384953/1‑55 (MQ=255)
                                                                                                                  aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                                                                                                                                                                                              >  6:2091/1‑56 (MQ=255)
                                                                                                                  aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                                                                                                                                                           >  8:1078246/1‑107 (MQ=255)
                                                                                                                   cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                                                                                                                                                             <  7:847743/104‑1 (MQ=255)
                                                                                                                     ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACg                                                                                                                                                                                                                 <  7:688215/82‑1 (MQ=255)
                                                                                                                     ttGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCtata                                                                                                                                                                          >  7:225266/1‑121 (MQ=255)
                                                                                                                     ttGTTAGTTGAAAGTACGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGaa                                                                                                                                                                                                                                          >  6:29998/1‑57 (MQ=255)
                                                                                                                        ttAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGtt                                                                          >  4:83127/1‑213 (MQ=255)
                                                                                                                         tAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg                                                                                                                                                                                <  6:945731/111‑1 (MQ=255)
                                                                                                                           cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                                                                                                                                                           >  7:782756/2‑98 (MQ=255)
                                                                                                                           cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                                                                                                                                                           <  8:782756/97‑1 (MQ=255)
                                                                                                                            ttGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACa                                                                                                                                                                                      <  6:640053/102‑1 (MQ=255)
                                                                                                                            ttGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATa                                                                                                                                                                                   <  2:356652/105‑1 (MQ=255)
                                                                                                                             tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                                                                                                                                                                                            >  8:84970/1‑63 (MQ=255)
                                                                                                                             tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAt                                                                                                                                                                                                                            <  7:84970/63‑1 (MQ=255)
                                                                                                                             tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                                                                                                                                                                                            <  8:695534/95‑1 (MQ=255)
                                                                                                                             tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCttt                                                                                                                                                                                            >  7:695534/1‑95 (MQ=255)
                                                                                                                             tGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACa                                                                                                                                                                                 <  4:40720/106‑1 (MQ=255)
                                                                                                                               aaaGTCTGGACAGTCTTTCAAGAGAAAAAGTGGGTGAAGGGATTGAATGTCTGAAATTGATTTTGCATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATaaa                                                                                                                                    <  8:600624/148‑1 (MQ=255)
                                                                                                                               aaaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                                                                                                                                                             <  6:981770/92‑1 (MQ=255)
                                                                                                                                aaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc                                                                                                                                                                                  >  5:23305/1‑102 (MQ=255)
                                                                                                                                aaGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc                                                                                                                                                                                  <  6:23305/102‑1 (MQ=255)
                                                                                                                                   tCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                                                                                                                                                          >  7:50817/1‑75 (MQ=255)
                                                                                                                                    cTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                                                                                                                                                                                                   >  5:956201/1‑65 (MQ=255)
                                                                                                                                       gACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAAt                                                                                                                                                                                    <  8:1142909/93‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                                                                                                                                             
ACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTTATAGCACAAAGGGCAAATGAAGAAAGCGAGATAAAGTCAAGTCGGGTTAAATTATCATGGAAGAAAAAACGGCAGGATGCACTGGAATCAGGTACGATTATGACGGCGTCTTGTCCGAGA  >  NC_000913/569025‑569429

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: