Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,280 C→G 6.0% intergenic (+33/‑35) renD → / → emrE pseudogene, DLP12 prophage;Phage or Prophage Related/DLP12 prophage; multidrug resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913568,2800CG 6.0% 423.5 / 11.6 130intergenic (+33/‑35)renD/emrEpseudogene, DLP12 prophage;Phage or Prophage Related/DLP12 prophage; multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base C (61/61);  new base G (6/2);  total (67/63)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 133 total aligned reads displayed.

AGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTT  >  NC_000913/568131‑568429
                                                                                                                                                      |                                                                                                                                                      
agagAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACAGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTAc                                                                                                                                                        >  2:1261955/1‑150 (MQ=255)
   gaagaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAAt                                                                                                                                                     <  1:456560/150‑1 (MQ=255)
    aagaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCTCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAAtt                                                                                                                                                    <  1:476490/150‑1 (MQ=255)
     agaaTGGGACGGAACGGTGAGAACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtc                                                                                                                                                   <  1:1631683/150‑1 (MQ=255)
     agaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGACGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGCAGAGTTTAGTTCAAACTGTAGTACAAtta                                                                                                                                                   >  2:722719/1‑149 (MQ=255)
      gaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtct                                                                                                                                                  >  2:1586099/1‑150 (MQ=255)
      gaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtct                                                                                                                                                  >  1:1280326/1‑150 (MQ=255)
       aaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtctc                                                                                                                                                <  2:542308/150‑1 (MQ=255)
        atcgcaaggacACGTGAGCACGAATCTGATTTTTAAGGAGGGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtctct                                                                                                                                               <  1:1222893/145‑1 (MQ=255)
          gggAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTcc                                                                                                                                             <  2:730759/150‑1 (MQ=255)
          gggAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTcc                                                                                                                                             <  1:715643/150‑1 (MQ=255)
           ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCTCGATGAAACGGGTATTCATGGTATATAAAAGAACATCAATGGGAACACACTGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCa                                                                                                                                            <  1:87406/150‑1 (MQ=255)
           ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGTTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCa                                                                                                                                            >  1:271230/1‑150 (MQ=255)
            gaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAg                                                                                                                                           >  1:1382570/1‑150 (MQ=255)
              aggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGtt                                                                                                                                         <  1:506950/150‑1 (MQ=255)
                gaaAGGTGAGCACGAATCTGATTTTTAAGTAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTg                                                                                                                                       <  1:989209/150‑1 (MQ=255)
                gaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTg                                                                                                                                       >  2:1006889/1‑150 (MQ=255)
                 aaaGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGa                                                                                                                                      <  2:760999/150‑1 (MQ=255)
                   aGGTGATCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGTGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAAc                                                                                                                                    <  2:1497990/150‑1 (MQ=255)
                    ggTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACa                                                                                                                                   <  2:1537677/150‑1 (MQ=255)
                    ggTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTAGATAAAAGAACATCAATGGGACCAGAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGaagc                                                                                                                                   >  2:1824668/1‑148 (MQ=255)
                     gTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAg                                                                                                                                  >  1:827435/1‑150 (MQ=255)
                      tGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTAGTGAGGGTATATAAAAGAACATCAATGGGAACCCAATGATGAAACAGGTGAGTTGAGTTAAAACTGTAGTACAATTATCTCCAGTTTGAACAgg                                                                                                                                 >  1:1582265/1‑150 (MQ=255)
                        ttcaataaTCTGATTTTTAAGGAGTCTAGCCAGAGTGCCGCGATTAAACGTGTATTGTGGGTATATAAAAGAAGATCAATGGGAACACTATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCAAGTTTGAACAGGaa                                                                                                                               <  2:730805/144‑1 (MQ=255)
                         gCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGaaa                                                                                                                              >  1:1061934/1‑150 (MQ=255)
                            cGAATCTGATTTTTTAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCATACTGTAGTACAATTCTCTCCAGTGTGAACAGGAAAGaa                                                                                                                           >  2:1148696/1‑150 (MQ=255)
                             gAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAat                                                                                                                          >  2:984011/1‑150 (MQ=255)
                                 cTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGTGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGc                                                                                                                      <  1:459265/150‑1 (MQ=255)
                                 cTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGc                                                                                                                      >  1:1658315/1‑150 (MQ=255)
                                 cTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGc                                                                                                                      <  1:1263957/150‑1 (MQ=255)
                                  tGATTTTTAGGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCt                                                                                                                     >  2:985551/1‑150 (MQ=255)
                                  tGATTTTTAAGGAGTGGCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCt                                                                                                                     <  2:176812/150‑1 (MQ=255)
                                   gATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTa                                                                                                                    <  2:1478636/150‑1 (MQ=255)
                                     ttttttAGGAGTGTCGCCAGAGGGCCGAGATGAAACGTGTATTGAGGTAATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                  <  1:703567/150‑1 (MQ=255)
                                     tttttAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                  <  1:1284484/150‑1 (MQ=255)
                                     tttttAAGGAGTGTCACCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                  >  1:391228/1‑150 (MQ=255)
                                      ttttAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGAATATAAAGGAACTGAATGGGAACCCCAAGGAGCAACCCGGGAGTTGAATTAAAAAATGATGTACATTAATATCAAGTTTAAAAAGGAAAGAaaatcttatta                                                                                                                 >  2:1093663/1‑143 (MQ=255)
                                             gAGTGTCGCCAGAGTGCCGCGATGAAACGGGTTTTGAGTGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTta                                                                                                          <  1:1777031/150‑1 (MQ=255)
                                              aGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTtat                                                                                                         >  1:1556792/1‑150 (MQ=255)
                                                   cGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTCGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTa                                                                                                    <  1:862176/150‑1 (MQ=255)
                                                    gCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTAt                                                                                                   >  2:836494/1‑150 (MQ=255)
                                                    gCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTAt                                                                                                   >  2:635736/1‑150 (MQ=255)
                                                         agTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtgg                                                                                              <  2:1471022/150‑1 (MQ=255)
                                                            gCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtggtgg                                                                                           >  1:1252084/1‑150 (MQ=255)
                                                              cgcgATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGgtggtg                                                                                         <  2:1382570/150‑1 (MQ=255)
                                                               gcgATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGc                                                                                        <  1:1006889/150‑1 (MQ=255)
                                                               gcgATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGc                                                                                        >  1:1255150/1‑150 (MQ=255)
                                                                    gAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATAc                                                                                   <  2:873693/150‑1 (MQ=255)
                                                                        cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                               >  1:1412814/1‑150 (MQ=255)
                                                                        cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                               >  2:787864/1‑150 (MQ=255)
                                                                        cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGTGCAATACTTGCa                                                                              <  2:1768401/150‑1 (MQ=255)
                                                                        cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAAATCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                               <  1:879611/150‑1 (MQ=255)
                                                                                gAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCa                                                                       >  2:591920/1‑150 (MQ=255)
                                                                                gAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGAAGAGGtaa                                                                       >  2:221484/1‑148 (MQ=255)
                                                                                  gtgTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCAtt                                                                     <  1:48887/148‑1 (MQ=255)
                                                                                   ggTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTg                                                                    >  1:242954/1‑150 (MQ=255)
                                                                                    gTATATAAAAGAACATAAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTgg                                                                   <  1:926001/150‑1 (MQ=255)
                                                                                        ttaAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACa                                                               >  1:702295/2‑150 (MQ=255)
                                                                                        ataAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACa                                                               >  1:925939/1‑150 (MQ=255)
                                                                                         taAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACaa                                                              <  1:462561/150‑1 (MQ=255)
                                                                                           aaaGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAAcc                                                            >  2:286950/1‑150 (MQ=255)
                                                                                                 cATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAAtt                                                      >  1:642839/1‑149 (MQ=255)
                                                                                                  aTCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGa                                                     >  1:1614893/1‑150 (MQ=255)
                                                                                                    cAATGGTAACACATTGATGAAACAGGTGATTTTTGTTCAAACTGTAGTACAATTGTCTCCAGTTTGAACAGGAAAGAATAGGCTATGAACCCTTTTATTTATCTTGGTGGTGCAATACTTGCAGAGTTCATTGGTACAACCTTAATGAAg                                                   <  1:596935/150‑1 (MQ=255)
                                                                                                     aaTGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGt                                                  <  2:866081/150‑1 (MQ=255)
                                                                                                       tGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAAACTTAATGAAGttt                                                >  2:796684/1‑150 (MQ=255)
                                                                                                       tGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGATATGAAACCTTATATTTATCTTGGTGGTGCAATACTAGCAGAGATCATTGGTACAACCTTAATGAAtttt                                                >  1:242572/1‑150 (MQ=255)
                                                                                                         ggAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTc                                              <  1:1261955/150‑1 (MQ=255)
                                                                                                                aaTGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATTCTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAggg                                       >  2:919865/1‑149 (MQ=255)
                                                                                                                  tgatgaAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTGCAGTTTGAACAGGAAAGAATATGCAATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGttt                                     <  2:1393499/150‑1 (MQ=255)
                                                                                                                   gatgaAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGtttt                                    >  1:319544/1‑150 (MQ=255)
                                                                                                                       aaaCAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTacac                                >  2:545827/1‑150 (MQ=255)
                                                                                                                        aaCAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACg                               <  2:1657991/150‑1 (MQ=255)
                                                                                                                             gTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTAt                          >  2:242924/1‑150 (MQ=255)
                                                                                                                               tagttgagttCAATCTGTAGTACTATTCTCTCCAGTTTGAACAGGATAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTGATTGGTACATCCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATgg                        <  1:956812/149‑1 (MQ=255)
                                                                                                                               tagttgagttCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATgg                        <  2:694464/149‑1 (MQ=255)
                                                                                                                               gagttgagttCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATgg                        >  1:1506047/1‑150 (MQ=255)
                                                                                                                                 gttgagttCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGcc                      >  2:1531345/1‑150 (MQ=255)
                                                                                                                                  ttgagttCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCa                     >  1:352350/1‑150 (MQ=255)
                                                                                                                                  tcctgttCAAACTGGAGAGAATTGTACTACAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACATTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCa                     >  1:1621522/5‑150 (MQ=38)
                                                                                                                                   tgagttCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCAt                    >  2:837662/1‑150 (MQ=255)
                                                                                                                                   cctgttCAAACTGGAGAGAATTGTACTACAGATTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACAAGGTTATGGCCaa                    >  1:1223289/4‑149 (MQ=255)
                                                                                                                                    ctgttCAAACTGGAGAGAATTGTACTACAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATc                   >  1:1013954/3‑150 (MQ=38)
                                                                                                                                    ctgttCAAACTGGAGAGAATTGTACTACAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGACATc                   >  1:1149866/3‑150 (MQ=38)
                                                                                                                                      gttCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTg                 >  1:808669/1‑150 (MQ=255)
                                                                                                                                      gttCAAACTGGAGAGAATTGTACTACAGTTTGAACAGGAAAGAAAATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAAAGAAGTTTTCAGAAGGTTTTACACTCTTATGGCCATCTg                 >  2:1447126/1‑150 (MQ=255)
                                                                                                                                           aaaCTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTgg             >  2:841864/1‑150 (MQ=255)
                                                                                                                                           aaaCTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATAAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTgg             <  2:1043879/150‑1 (MQ=255)
                                                                                                                                             aCTGTAGTACAATTCTCTCCATTTTGAACAGGAAAGAATATGATATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTa           <  2:449378/150‑1 (MQ=255)
                                                                                                                                              cTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTAc          <  2:1026337/150‑1 (MQ=255)
                                                                                                                                              cTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTAc          <  1:836494/150‑1 (MQ=255)
                                                                                                                                                gtagtaCAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACaa        <  2:1061934/150‑1 (MQ=255)
                                                                                                                                                   gtaCAATTCTCTCCATTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta     <  1:107747/150‑1 (MQ=255)
                                                                                                                                                   gtaCAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta     >  1:1572842/1‑150 (MQ=255)
                                                                                                                                                   gtaCAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta     >  1:1680725/1‑150 (MQ=255)
                                                                                                                                                   gtaCAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACAGTTATGGCCATCTGTTGGTACAatta     >  2:988326/1‑150 (MQ=255)
                                                                                                                                                     aCAATTCTCTCCAGTTTGTACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGTTTATGGCCATCTGTTGGTACAattatt   <  2:1653858/150‑1 (MQ=255)
                                                                                                                                                     aCAATTCTCTCCAGTTTGAACAGGAATGAATATGCTATGATCCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAattatt   <  2:1690685/150‑1 (MQ=255)
                                                                                                                                                      cAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTAttt  <  1:635736/150‑1 (MQ=255)
                                                                                                                                                      cAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTAttt  >  1:99985/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTT  >  NC_000913/568131‑568429

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