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.

AGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTT  >  NC_000913/568131‑568429
                                                                                                                                                      |                                                                                                                                                      
agagAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACAGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTAc                                                                                                                                                        >  2:1261955/1‑150 (MQ=255)
   gaagaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAAt                                                                                                                                                     <  1:456560/150‑1 (MQ=255)
    aagaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCTCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAAtt                                                                                                                                                    <  1:476490/150‑1 (MQ=255)
     agaaTGGGACGGAACGGTGAGAACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATtc                                                                                                                                                   <  1:1631683/150‑1 (MQ=255)
     agaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGACGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGCAGAGTTTAGTT‑CAAACTGTAGTACAAtta                                                                                                                                                   >  2:722719/1‑149 (MQ=255)
      gaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATtct                                                                                                                                                  >  2:1586099/1‑150 (MQ=255)
      gaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATtct                                                                                                                                                  >  1:1280326/1‑150 (MQ=255)
       aaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATtct‑c                                                                                                                                                <  2:542308/150‑1 (MQ=255)
        atcgcaaggacACGTGAGCACGAATCTGATTTTTAAGGAGGGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATtct‑ct                                                                                                                                               <  1:1222893/145‑1 (MQ=255)
          gggAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTcc                                                                                                                                             <  2:730759/150‑1 (MQ=255)
          gggAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTcc                                                                                                                                             <  1:715643/150‑1 (MQ=255)
           ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCTCGATGAAACGGGTATTCATGGTATATAAAAGAACATCAATGGGAACACACTGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCa                                                                                                                                            <  1:87406/150‑1 (MQ=255)
           ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGTTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCa                                                                                                                                            >  1:271230/1‑150 (MQ=255)
            gaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAg                                                                                                                                           >  1:1382570/1‑150 (MQ=255)
              aggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGtt                                                                                                                                         <  1:506950/150‑1 (MQ=255)
                gaaAGGTGAGCACGAATCTGATTTTTAAGTAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTg                                                                                                                                       <  1:989209/150‑1 (MQ=255)
                gaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTg                                                                                                                                       >  2:1006889/1‑150 (MQ=255)
                 aaaGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGa                                                                                                                                      <  2:760999/150‑1 (MQ=255)
                   aGGTGATCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGTGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAAc                                                                                                                                    <  2:1497990/150‑1 (MQ=255)
                    ggTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACa                                                                                                                                   <  2:1537677/150‑1 (MQ=255)
                    ggTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTAGATAAAAGAACATCAATGGGACCAGAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGaagc                                                                                                                                   >  2:1824668/1‑148 (MQ=255)
                     gTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAg                                                                                                                                  >  1:827435/1‑150 (MQ=255)
                      tGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTAGTGAGGGTATATAAAAGAACATCAATGGGAACCCAATGATGAAACAGGTGAGTTGAGTT‑AAAACTGTAGTACAATTAT‑CTCCAGTTTGAACAgg                                                                                                                                 >  1:1582265/1‑150 (MQ=255)
                        ttcaataaTCTGATTTTTAAGGAGTCTAGCCAGAGTGCCGCGATTAAACGTGTATTGTGGGTATATAAAAGAAGATCAATGGGAACACTATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCAAGTTTGAACAGGaa                                                                                                                               <  2:730805/144‑1 (MQ=255)
                         gCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGaaa                                                                                                                              >  1:1061934/1‑150 (MQ=255)
                            cGAATCTGATTTTTTAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CATACTGTAGTACAATTCT‑CTCCAGTGTGAACAGGAAAGaa                                                                                                                           >  2:1148696/1‑150 (MQ=255)
                             gAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAat                                                                                                                          >  2:984011/1‑150 (MQ=255)
                                 cTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGTGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGc                                                                                                                      <  1:459265/150‑1 (MQ=255)
                                 cTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGc                                                                                                                      >  1:1658315/1‑150 (MQ=255)
                                 cTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGc                                                                                                                      <  1:1263957/150‑1 (MQ=255)
                                  tGATTTTTAGGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCt                                                                                                                     >  2:985551/1‑150 (MQ=255)
                                  tGATTTTTAAGGAGTGGCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCt                                                                                                                     <  2:176812/150‑1 (MQ=255)
                                   gATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTa                                                                                                                    <  2:1478636/150‑1 (MQ=255)
                                     ttttttAGGAGTGTCGCCAGAGGGCCGAGATGAAACGTGTATTGAGGTAATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                  <  1:703567/150‑1 (MQ=255)
                                     tttttAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                  <  1:1284484/150‑1 (MQ=255)
                                     tttttAAGGAGTGTCACCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                  >  1:391228/1‑150 (MQ=255)
                                      ttttAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGAATATAAAGGAAC‑TGAATGGGAACCCCAAGGAGCAACCCGGGAGTTGAATTAAAAAATGATGTACATTAAT‑ATCAAGTTTAAAAAGGAAAGAaaatcttatta                                                                                                                 >  2:1093663/1‑143 (MQ=255)
                                             gAGTGTCGCCAGAGTGCCGCGATGAAACGGGTTTTGAGTGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTta                                                                                                          <  1:1777031/150‑1 (MQ=255)
                                              aGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTtat                                                                                                         >  1:1556792/1‑150 (MQ=255)
                                                   cGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTCGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTa                                                                                                    <  1:862176/150‑1 (MQ=255)
                                                    gCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTAt                                                                                                   >  2:836494/1‑150 (MQ=255)
                                                    gCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTAt                                                                                                   >  2:635736/1‑150 (MQ=255)
                                                         agTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtgg                                                                                              <  2:1471022/150‑1 (MQ=255)
                                                            gCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtggtgg                                                                                           >  1:1252084/1‑150 (MQ=255)
                                                              cgcgATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGgtggtg                                                                                         <  2:1382570/150‑1 (MQ=255)
                                                               gcgATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGc                                                                                        <  1:1006889/150‑1 (MQ=255)
                                                               gcgATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGc                                                                                        >  1:1255150/1‑150 (MQ=255)
                                                                    gAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATAc                                                                                   <  2:873693/150‑1 (MQ=255)
                                                                        cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                               >  1:1412814/1‑150 (MQ=255)
                                                                        cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                               >  2:787864/1‑150 (MQ=255)
                                                                        cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGT‑GTGCAATACTTGCa                                                                              <  2:1768401/150‑1 (MQ=255)
                                                                        cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAAATCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                               <  1:879611/150‑1 (MQ=255)
                                                                                gAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCa                                                                       >  2:591920/1‑150 (MQ=255)
                                                                                gAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGAAGAGGtaa                                                                       >  2:221484/1‑148 (MQ=255)
                                                                                  gtgTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCAtt                                                                     <  1:48887/148‑1 (MQ=255)
                                                                                   ggTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTg                                                                    >  1:242954/1‑150 (MQ=255)
                                                                                    gTATATAAAAGAACATAAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTgg                                                                   <  1:926001/150‑1 (MQ=255)
                                                                                        ttaAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACa                                                               >  1:702295/2‑150 (MQ=255)
                                                                                        ataAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACa                                                               >  1:925939/1‑150 (MQ=255)
                                                                                         taAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACaa                                                              <  1:462561/150‑1 (MQ=255)
                                                                                           aaaGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAAcc                                                            >  2:286950/1‑150 (MQ=255)
                                                                                                 cATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAAtt                                                      >  1:642839/1‑149 (MQ=255)
                                                                                                  aTCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGa                                                     >  1:1614893/1‑150 (MQ=255)
                                                                                                    cAATGGTAACACATTGATGAAACAGGTGATTTTTGTT‑CAAACTGTAGTACAATTGT‑CTCCAGTTTGAACAGGAAAGAATAGGCTATGAACCCTTTTATTTATCTTGGTGGTGCAATACTTGCAGAGTTCATTGGTACAACCTTAATGAAg                                                   <  1:596935/150‑1 (MQ=255)
                                                                                                     aaTGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGt                                                  <  2:866081/150‑1 (MQ=255)
                                                                                                       tGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAAACTTAATGAAGttt                                                >  2:796684/1‑150 (MQ=255)
                                                                                                       tGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGATATGAAACCTTATATTTATCTTGGTGGTGCAATACTAGCAGAGATCATTGGTACAACCTTAATGAAtttt                                                >  1:242572/1‑150 (MQ=255)
                                                                                                         ggAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTc                                              <  1:1261955/150‑1 (MQ=255)
                                                                                                                aaTGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATTCTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAggg                                       >  2:919865/1‑149 (MQ=255)
                                                                                                                  tgatgaAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTGCAGTTTGAACAGGAAAGAATATGCAATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGttt                                     <  2:1393499/150‑1 (MQ=255)
                                                                                                                   gatgaAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGtttt                                    >  1:319544/1‑150 (MQ=255)
                                                                                                                       aaaCAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTacac                                >  2:545827/1‑150 (MQ=255)
                                                                                                                        aaCAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACg                               <  2:1657991/150‑1 (MQ=255)
                                                                                                                             gTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTAt                          >  2:242924/1‑150 (MQ=255)
                                                                                                                               tagttgagtt‑CAATCTGTAGTACTATTCT‑CTCCAGTTTGAACAGGATAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTGATTGGTACATCCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATgg                        <  1:956812/149‑1 (MQ=255)
                                                                                                                               tagttgagtt‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATgg                        <  2:694464/149‑1 (MQ=255)
                                                                                                                               gagttgagtt‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATgg                        >  1:1506047/1‑150 (MQ=255)
                                                                                                                                 gttgagtt‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGcc                      >  2:1531345/1‑150 (MQ=255)
                                                                                                                                  ttgagtt‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCa                     >  1:352350/1‑150 (MQ=255)
                                                                                                                                  tcctgtt‑CAAACTGGAG‑AGAATTGTACTACAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACATTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCa                     >  1:1621522/5‑150 (MQ=38)
                                                                                                                                   tgagtt‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCAt                    >  2:837662/1‑150 (MQ=255)
                                                                                                                                   cctgtt‑CAAACTGGAG‑AGAATTGTACTACAGATTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACAAGGTTATGGCCaa                    >  1:1223289/4‑149 (MQ=255)
                                                                                                                                    ctgtt‑CAAACTGGAG‑AGAATTGTACTACAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATc                   >  1:1013954/3‑150 (MQ=38)
                                                                                                                                    ctgtt‑CAAACTGGAG‑AGAATTGTACTACAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGACATc                   >  1:1149866/3‑150 (MQ=38)
                                                                                                                                      gtt‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTg                 >  1:808669/1‑150 (MQ=255)
                                                                                                                                      gtt‑CAAACTGGAG‑AGAATTGTACTACAGTTTGAACAGGAAAGAAAATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAAAGAAGTTTTCAGAAGGTTTTACACTCTTATGGCCATCTg                 >  2:1447126/1‑150 (MQ=255)
                                                                                                                                           aaaCTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTgg             >  2:841864/1‑150 (MQ=255)
                                                                                                                                           aaaCTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATAAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTgg             <  2:1043879/150‑1 (MQ=255)
                                                                                                                                             aCTGTAGTACAATTCT‑CTCCATTTTGAACAGGAAAGAATATGATATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTa           <  2:449378/150‑1 (MQ=255)
                                                                                                                                              cTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTAc          <  2:1026337/150‑1 (MQ=255)
                                                                                                                                              cTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTAc          <  1:836494/150‑1 (MQ=255)
                                                                                                                                                gtagtaCAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACaa        <  2:1061934/150‑1 (MQ=255)
                                                                                                                                                   gtaCAATTCT‑CTCCATTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta     <  1:107747/150‑1 (MQ=255)
                                                                                                                                                   gtaCAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta     >  1:1572842/1‑150 (MQ=255)
                                                                                                                                                   gtaCAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta     >  1:1680725/1‑150 (MQ=255)
                                                                                                                                                   gtaCAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACAGTTATGGCCATCTGTTGGTACAatta     >  2:988326/1‑150 (MQ=255)
                                                                                                                                                     aCAATTCT‑CTCCAGTTTGTACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGTTTATGGCCATCTGTTGGTACAattatt   <  2:1653858/150‑1 (MQ=255)
                                                                                                                                                     aCAATTCT‑CTCCAGTTTGAACAGGAATGAATATGCTATGATCCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAattatt   <  2:1690685/150‑1 (MQ=255)
                                                                                                                                                      cAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTAttt  <  1:635736/150‑1 (MQ=255)
                                                                                                                                                      cAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTAttt  >  1:99985/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTT  >  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: ‑