Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,285 C→G 6.2% intergenic (+38/‑30) 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,2850CG 6.2% 400.1 / 10.7 123intergenic (+38/‑30)renD/emrEpseudogene, DLP12 prophage;Phage or Prophage Related/DLP12 prophage; multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base C (56/57);  new base G (6/2);  total (64/59)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTAT  >  NC_000913/568136‑568434
                                                                                                                                                      |                                                                                                                                                      
agaaTGGGACGGAACGGTGAGAACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATtc                                                                                                                                                        <  1:1631683/150‑1 (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)
    tGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATtct‑ctc                                                                                                                                                   >  2:563869/1‑150 (MQ=255)
     gggAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTcc                                                                                                                                                  <  1:715643/150‑1 (MQ=255)
     gggAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTcc                                                                                                                                                  <  2:730759/150‑1 (MQ=255)
      ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCTCGATGAAACGGGTATTCATGGTATATAAAAGAACATCAATGGGAACACACTGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCa                                                                                                                                                 <  1:87406/150‑1 (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)
              aGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAAAATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAAc                                                                                                                                         <  2:1505068/150‑1 (MQ=255)
               ggTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTAGATAAAAGAACATCAATGGGACCAGAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGaagc                                                                                                                                        >  2:1824668/1‑148 (MQ=255)
                gTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAg                                                                                                                                       >  1:827435/1‑150 (MQ=255)
                 tGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGGAG‑AGAATTGTACTACAGTTTGAActca                                                                                                                                      <  1:181693/150‑4 (MQ=38)
                 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)
                    gCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGaaa                                                                                                                                   >  1:595339/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)
                                ttttttAGGAGTGTCGCCAGAGGGCCGAGATGAAACGTGTATTGAGGTAATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                       <  1:703567/150‑1 (MQ=255)
                                tttttAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                       <  1:1284484/150‑1 (MQ=255)
                                 ttttAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGAATATAAAGGAAC‑TGAATGGGAACCCCAAGGAGCAACCCGGGAGTTGAATTAAAAAATGATGTACATTAAT‑ATCAAGTTTAAAAAGGAAAGAaaatcttatta                                                                                                                      >  2:1093663/1‑143 (MQ=255)
                                    tAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAAcc                                                                                                                   <  1:984011/150‑1 (MQ=255)
                                        gAGTGTCTCCAGAGTGCCGCGATGATACGTGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTTAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTta                                                                                                               <  1:91567/150‑1 (MQ=255)
                                        gAGTGTCGCCAGAGTGCCGCGATGAAACGGGTTTTGAGTGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTta                                                                                                               <  1:1777031/150‑1 (MQ=255)
                                         aGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTtat                                                                                                              >  1:1556792/1‑150 (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)
                                                      tGCCGCGATTAAACGGGTATTGAGGGTATATAAAATAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtggtg                                                                                                 <  1:1647304/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                                                                                             <  2:1015290/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)
                                                                 aaCGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAAATCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACtt                                                                                      >  2:707651/1‑150 (MQ=255)
                                                                   cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                                    >  2:787864/1‑150 (MQ=255)
                                                                   cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                                    >  1:1412814/1‑150 (MQ=255)
                                                                   cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAAATCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                                    <  1:879611/150‑1 (MQ=255)
                                                                        aTTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAgg                                                                               <  2:1694323/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)
                                                                                   ataAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACa                                                                    >  1:1601295/1‑150 (MQ=255)
                                                                                    taAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACaa                                                                   <  1:462561/150‑1 (MQ=255)
                                                                                    taAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACaa                                                                   >  1:644448/1‑150 (MQ=255)
                                                                                            cATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAAtt                                                           >  1:642839/1‑149 (MQ=255)
                                                                                            cATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATg                                                           >  2:1464131/1‑150 (MQ=255)
                                                                                             aTCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGa                                                          >  1:1614893/1‑150 (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)
                                                                                                    ttgcgccaATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTc                                                   <  2:1589335/144‑1 (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)
                                                                                                                    aCAGGTGAGGTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACgg                                   <  2:1087638/150‑1 (MQ=255)
                                                                                                                     cAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGt                                  <  1:1528156/150‑1 (MQ=255)
                                                                                                                       ggTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTTCACGGttt                                >  2:204049/1‑149 (MQ=255)
                                                                                                                        gTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTAt                               >  2:242924/1‑150 (MQ=255)
                                                                                                                         tgagttgagtT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATg                              <  1:297750/150‑1 (MQ=255)
                                                                                                                          tagttgagtt‑CAATCTGTAGTACTATTCT‑CTCCAGTTTGAACAGGATAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTGATTGGTACATCCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATgg                             <  1:956812/149‑1 (MQ=255)
                                                                                                                            gttgagtt‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGcc                           >  2:1531345/1‑150 (MQ=255)
                                                                                                                             ttgagtt‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCa                          >  1:352350/1‑150 (MQ=255)
                                                                                                                              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‑CAAACTGGAG‑AGAATTGTACTACAGTTTGAACAGGAAAGAAAATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAAAGAAGTTTTCAGAAGGTTTTACACTCTTATGGCCATCTg                      >  2:1447126/1‑150 (MQ=255)
                                                                                                                                     cAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTg                   >  1:425169/1‑150 (MQ=255)
                                                                                                                                      aaaCTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTgg                  >  1:84175/1‑150 (MQ=255)
                                                                                                                                      aaaCTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTgg                  >  2:841864/1‑150 (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)
                                                                                                                                           gtagtaCAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGATGGTACaa             >  2:126440/1‑150 (MQ=255)
                                                                                                                                             agtaGAATTCT‑CTCCAGTTTGAACAGGAAAGAATTTGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatt           <  1:722719/150‑1 (MQ=255)
                                                                                                                                             agtaCAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatt           >  2:385352/1‑150 (MQ=255)
                                                                                                                                              gtaCAATTCT‑CTCCATTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta          <  1:107747/150‑1 (MQ=255)
                                                                                                                                              gtaCAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta          >  1:1572842/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)
                                                                                                                                                     tct‑ctcCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTa   <  1:1464131/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTT‑CAAACTGTAGTACAATTCT‑CTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTAT  >  NC_000913/568136‑568434

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