Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 569,161 A→C 5.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 5.1% 408.0 / 12.9 120D87A (GAC→GCC) ‡ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base A (56/57);  new base C (4/2);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 123 total aligned reads displayed.

TATCTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTT  >  NC_000913/569016‑569309
                                                                                                                                                 |                                                                                                                                                    
  tCTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGa                                                                                                                                                  >  2:296567/1‑148 (MQ=255)
  tCTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGAc                                                                                                                                                 >  3:179401/1‑149 (MQ=255)
      gATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTc                                                                                                                                             <  3:341912/149‑1 (MQ=255)
      gATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                                                                              >  2:146397/1‑148 (MQ=255)
      gATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                                                                              <  1:664684/148‑1 (MQ=255)
      gATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACttt                                                                                                                                              >  2:255374/1‑148 (MQ=255)
       aTACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCAGAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                                                                             >  1:625113/1‑148 (MQ=255)
        tACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCaa                                                                                                                                           >  3:898232/1‑149 (MQ=255)
        tACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                                                                            <  1:330532/148‑1 (MQ=255)
         aCCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                          >  3:145131/1‑149 (MQ=255)
               aCATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaaa                                                                                                                                    <  3:595894/149‑1 (MQ=255)
                cATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGaaaaa                                                                                                                                   <  3:641938/149‑1 (MQ=255)
                  tatGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGATGAAAGTCTGGACAGACTTTCAAGAGAAAAAGt                                                                                                                                 >  4:632887/1‑149 (MQ=255)
                    tGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTc                                                                                                                                <  1:294733/148‑1 (MQ=255)
                    tGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTc                                                                                                                                <  2:424698/148‑1 (MQ=255)
                    tGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCg                                                                                                                               <  4:319539/149‑1 (MQ=255)
                     gAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCg                                                                                                                               >  1:288579/1‑148 (MQ=255)
                     gAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCg                                                                                                                               <  1:746714/148‑1 (MQ=255)
                         aTTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGTACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGa                                                                                                                           >  2:788557/1‑148 (MQ=255)
                            tAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGTCCAGACTTTCa                                                                                                                                            >  2:279553/1‑128 (MQ=255)
                             aGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGa                                                                                                                      >  4:541854/1‑149 (MQ=255)
                               gATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAtt                                                                                                                    <  4:499509/149‑1 (MQ=255)
                               gATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                                                                     >  1:771479/1‑148 (MQ=255)
                                  tAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGa                                                                                                                  >  2:596460/1‑148 (MQ=255)
                                  tAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGa                                                                                                                  >  2:798411/1‑148 (MQ=255)
                                      tGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                          <  2:436327/120‑1 (MQ=255)
                                      tGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAag                                                                                                                                          >  1:436327/1‑120 (MQ=255)
                                         aTTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGa                                                                                                          >  4:464436/1‑149 (MQ=255)
                                          ttCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                                                                             <  4:760813/113‑1 (MQ=255)
                                          ttCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTc                                                                                                                                             >  3:760813/1‑113 (MQ=255)
                                           tCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAGGAGCAAAAGTCGGTGAAGCGATTGAACGTCTGaa                                                                                                         >  2:738074/1‑148 (MQ=255)
                                           tCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaa                                                                                                         >  2:145810/1‑148 (MQ=255)
                                           tCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGaa                                                                                                         <  2:567027/148‑1 (MQ=255)
                                             aaaGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAAtt                                                                                                      >  3:8161/1‑149 (MQ=255)
                                               aGGAAAGCATGCACAATCAGGAGTTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                                                                            >  2:81579/1‑109 (MQ=255)
                                               aGGAAAGCATGCACAATCAGGAGTTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCa                                                                                                                                            <  1:81579/109‑1 (MQ=255)
                                                ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGa                                                                                                    <  2:226616/148‑1 (MQ=255)
                                                ggAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGa                                                                                                    >  2:653290/1‑148 (MQ=255)
                                                   aaGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGAttt                                                                                                 <  1:379503/148‑1 (MQ=255)
                                                       aTGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAAt                                                                                            >  3:192075/1‑149 (MQ=255)
                                                          cacaATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCa                                                                                          <  2:851596/148‑1 (MQ=255)
                                                              aTCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTa                                                                                     <  4:179401/149‑1 (MQ=255)
                                                               tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGa                                                                                                    <  2:890544/133‑1 (MQ=255)
                                                               tCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGa                                                                                                    >  1:890544/1‑133 (MQ=255)
                                                               tCAGGAGCTTTTTCGAAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATAACGGTa                                                                                     <  2:517509/148‑1 (MQ=255)
                                                                   gAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGa                                                                                 <  1:841603/148‑1 (MQ=255)
                                                                    aGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGAt                                                                                <  1:541476/148‑1 (MQ=255)
                                                                    aGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATg                                                                               >  3:675584/1‑149 (MQ=255)
                                                                      cTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAAAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACAGTATTGATGtt                                                                             <  3:541854/149‑1 (MQ=255)
                                                                        ttttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTtat                                                                           <  4:22977/149‑1 (MQ=255)
                                                                          ttCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATaa                                                                         <  4:233991/149‑1 (MQ=255)
                                                                             ggAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGa                                                                                                                                                  <  3:913128/73‑1 (MQ=255)
                                                                             ggAATTTTTAGATGCGATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGa                                                                                                                                                  >  4:913128/1‑73 (MQ=255)
                                                                                aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                                    >  2:708772/1‑148 (MQ=255)
                                                                                aTTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCtt                                                                    <  2:625113/148‑1 (MQ=255)
                                                                                 tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                                  <  3:908343/149‑1 (MQ=255)
                                                                                 tttttAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTg                                                                  >  3:470390/1‑149 (MQ=255)
                                                                                          gCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                                                                     <  1:234567/89‑1 (MQ=255)
                                                                                          gCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGAt                                                                                                                     >  2:234567/1‑89 (MQ=255)
                                                                                           cTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATAc                                                         <  2:958123/148‑1 (MQ=255)
                                                                                            tataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACa                                                        <  1:243159/148‑1 (MQ=255)
                                                                                             ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGt                                                      <  4:102436/149‑1 (MQ=255)
                                                                                             ataGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCCGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAg                                                       >  1:608621/1‑148 (MQ=255)
                                                                                              taGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTTAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTc                                                     <  3:313070/149‑1 (MQ=255)
                                                                                                 agCATGGTTATATATTGCCAGGAAGTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATGACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCtaa                                                  <  3:297663/149‑2 (MQ=255)
                                                                                                  gCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCtat                                                  >  1:501755/1‑148 (MQ=255)
                                                                                                   cATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCtata                                                 <  1:422415/148‑1 (MQ=255)
                                                                                                      ggTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAatat                                             <  3:340959/149‑1 (MQ=255)
                                                                                                       gTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATTACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAatat                                             <  1:733656/148‑1 (MQ=255)
                                                                                                       gTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAatat                                             >  2:366790/1‑148 (MQ=255)
                                                                                                        ttATATATTGCCAGGAACTATATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTg                                           >  3:318122/1‑149 (MQ=255)
                                                                                                        ttATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATg                                                                               >  4:283622/1‑113 (MQ=255)
                                                                                                        ttATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATg                                                                               <  3:283622/113‑1 (MQ=255)
                                                                                                            atatTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTTAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGACTCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGActc                                       <  3:668444/149‑1 (MQ=255)
                                                                                                               ttGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCttt                                    <  3:366349/149‑1 (MQ=255)
                                                                                                               ttGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCtt                                     <  1:255374/148‑1 (MQ=255)
                                                                                                                  ccAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGaa                                 <  4:145131/149‑1 (MQ=255)
                                                                                                                   cAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGAct                                        >  3:869/1‑141 (MQ=255)
                                                                                                                   cAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGAct                                        <  4:869/141‑1 (MQ=255)
                                                                                                                   cAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAAt                                >  3:423064/1‑149 (MQ=255)
                                                                                                                   cAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAAt                                >  4:833236/1‑149 (MQ=255)
                                                                                                                     ggAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATg                               <  1:780944/148‑1 (MQ=255)
                                                                                                                     ggAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATg                               <  1:622118/148‑1 (MQ=255)
                                                                                                                       aaCTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGCAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAg                             >  2:168960/1‑148 (MQ=255)
                                                                                                                          tACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCat                         >  4:267796/1‑149 (MQ=255)
                                                                                                                           aCATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCTCGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCata                        >  4:479379/1‑149 (MQ=255)
                                                                                                                            cATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCata                        >  2:155152/1‑148 (MQ=255)
                                                                                                                                gTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCa                    <  1:348270/148‑1 (MQ=255)
                                                                                                                                gTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCa                    >  2:47363/1‑148 (MQ=255)
                                                                                                                                gTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCAt                   <  4:470390/149‑1 (MQ=255)
                                                                                                                                 ttAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACa                                                        <  4:555480/111‑1 (MQ=255)
                                                                                                                                 ttAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACa                                                        >  3:555480/1‑111 (MQ=255)
                                                                                                                                   aGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATtaa                <  4:898232/149‑1 (MQ=255)
                                                                                                                                    cttGAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATtaa                >  1:894040/2‑148 (MQ=255)
                                                                                                                                       gAAAGTCTGTCCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATtaataa             >  1:941370/1‑148 (MQ=255)
                                                                                                                                              tGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATa     <  4:169735/149‑1 (MQ=255)
                                                                                                                                               ggACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATa     <  1:218788/148‑1 (MQ=255)
                                                                                                                                                 aCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACtt  >  3:603369/1‑149 (MQ=255)
                                                                                                                                                 aCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACt   >  2:381729/1‑148 (MQ=255)
                                                                                                                                                 aCAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACt   >  1:15019/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TATCTTGATACCATAACATATGAAGATTTAGGATTAAGTGCATTCAAAGGAAAGCATGCACAATCAGGAGCTTTTTCGGAATTTTTAGATGCTATAGAGCATGGTTATATATTGCCAGGAACTACATTGTTAGTTGAAAGTCTGGACAGACTTTCAAGAGAAAAAGTCGGTGAAGCGATTGAACGTCTGAAATTGATTTTGAATCACGGTATTGATGTTATAACTCTTTGCGACAATACAGTCTATAATATTGACTCTTTGAATGAGCCATATTCATTAATAAAAGCCATACTT  >  NC_000913/569016‑569309

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