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.

AGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTAT  >  NC_000913/568136‑568434
                                                                                                                                                      |                                                                                                                                                      
agaaTGGGACGGAACGGTGAGAACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtc                                                                                                                                                        <  1:1631683/150‑1 (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)
    tGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtctctc                                                                                                                                                   >  2:563869/1‑150 (MQ=255)
     gggAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTcc                                                                                                                                                  <  1:715643/150‑1 (MQ=255)
     gggAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTcc                                                                                                                                                  <  2:730759/150‑1 (MQ=255)
      ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCTCGATGAAACGGGTATTCATGGTATATAAAAGAACATCAATGGGAACACACTGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCa                                                                                                                                                 <  1:87406/150‑1 (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)
              aGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAAAATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAAc                                                                                                                                         <  2:1505068/150‑1 (MQ=255)
               ggTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTAGATAAAAGAACATCAATGGGACCAGAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGaagc                                                                                                                                        >  2:1824668/1‑148 (MQ=255)
                gTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAg                                                                                                                                       >  1:827435/1‑150 (MQ=255)
                 tGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGGAGAGAATTGTACTACAGTTTGAActca                                                                                                                                      <  1:181693/150‑4 (MQ=38)
                 tGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTAGTGAGGGTATATAAAAGAACATCAATGGGAACCCAATGATGAAACAGGTGAGTTGAGTTAAAACTGTAGTACAATTATCTCCAGTTTGAACAgg                                                                                                                                      >  1:1582265/1‑150 (MQ=255)
                   ttcaataaTCTGATTTTTAAGGAGTCTAGCCAGAGTGCCGCGATTAAACGTGTATTGTGGGTATATAAAAGAAGATCAATGGGAACACTATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCAAGTTTGAACAGGaa                                                                                                                                    <  2:730805/144‑1 (MQ=255)
                    gCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGaaa                                                                                                                                   >  1:1061934/1‑150 (MQ=255)
                    gCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGaaa                                                                                                                                   >  1:595339/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)
                                ttttttAGGAGTGTCGCCAGAGGGCCGAGATGAAACGTGTATTGAGGTAATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                       <  1:703567/150‑1 (MQ=255)
                                tttttAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATg                                                                                                                       <  1:1284484/150‑1 (MQ=255)
                                 ttttAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGAATATAAAGGAACTGAATGGGAACCCCAAGGAGCAACCCGGGAGTTGAATTAAAAAATGATGTACATTAATATCAAGTTTAAAAAGGAAAGAaaatcttatta                                                                                                                      >  2:1093663/1‑143 (MQ=255)
                                    tAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAAcc                                                                                                                   <  1:984011/150‑1 (MQ=255)
                                        gAGTGTCTCCAGAGTGCCGCGATGATACGTGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTTAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTta                                                                                                               <  1:91567/150‑1 (MQ=255)
                                        gAGTGTCGCCAGAGTGCCGCGATGAAACGGGTTTTGAGTGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTta                                                                                                               <  1:1777031/150‑1 (MQ=255)
                                         aGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTtat                                                                                                              >  1:1556792/1‑150 (MQ=255)
                                               gCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTAt                                                                                                        >  2:836494/1‑150 (MQ=255)
                                               gCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTAt                                                                                                        >  2:635736/1‑150 (MQ=255)
                                                    agTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtgg                                                                                                   <  2:1471022/150‑1 (MQ=255)
                                                      tGCCGCGATTAAACGGGTATTGAGGGTATATAAAATAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtggtg                                                                                                 <  1:1647304/150‑1 (MQ=255)
                                                       gCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTtggtgg                                                                                                >  1:1252084/1‑150 (MQ=255)
                                                         cgcgATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGgtggtg                                                                                              <  2:1382570/150‑1 (MQ=255)
                                                          gcgATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGc                                                                                             <  2:1015290/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)
                                                                 aaCGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAAATCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACtt                                                                                      >  2:707651/1‑150 (MQ=255)
                                                                   cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                                    >  2:787864/1‑150 (MQ=255)
                                                                   cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                                    >  1:1412814/1‑150 (MQ=255)
                                                                   cGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAAATCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGc                                                                                    <  1:879611/150‑1 (MQ=255)
                                                                        aTTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAgg                                                                               <  2:1694323/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)
                                                                                   ataAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACa                                                                    >  1:1601295/1‑150 (MQ=255)
                                                                                    taAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACaa                                                                   <  1:462561/150‑1 (MQ=255)
                                                                                    taAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACaa                                                                   >  1:644448/1‑150 (MQ=255)
                                                                                            cATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAAtt                                                           >  1:642839/1‑149 (MQ=255)
                                                                                            cATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATg                                                           >  2:1464131/1‑150 (MQ=255)
                                                                                             aTCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGa                                                          >  1:1614893/1‑150 (MQ=255)
                                                                                                aaTGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGt                                                       <  2:866081/150‑1 (MQ=255)
                                                                                                  tGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAAACTTAATGAAGttt                                                     >  2:796684/1‑150 (MQ=255)
                                                                                                  tGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGATATGAAACCTTATATTTATCTTGGTGGTGCAATACTAGCAGAGATCATTGGTACAACCTTAATGAAtttt                                                     >  1:242572/1‑150 (MQ=255)
                                                                                                    ttgcgccaATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTc                                                   <  2:1589335/144‑1 (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)
                                                                                                                    aCAGGTGAGGTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACgg                                   <  2:1087638/150‑1 (MQ=255)
                                                                                                                     cAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGt                                  <  1:1528156/150‑1 (MQ=255)
                                                                                                                       ggTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTTCACGGttt                                >  2:204049/1‑149 (MQ=255)
                                                                                                                        gTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTAt                               >  2:242924/1‑150 (MQ=255)
                                                                                                                         tgagttgagtTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATg                              <  1:297750/150‑1 (MQ=255)
                                                                                                                          tagttgagttCAATCTGTAGTACTATTCTCTCCAGTTTGAACAGGATAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTGATTGGTACATCCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATgg                             <  1:956812/149‑1 (MQ=255)
                                                                                                                            gttgagttCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGcc                           >  2:1531345/1‑150 (MQ=255)
                                                                                                                             ttgagttCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCa                          >  1:352350/1‑150 (MQ=255)
                                                                                                                              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)
                                                                                                                                 gttCAAACTGGAGAGAATTGTACTACAGTTTGAACAGGAAAGAAAATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAAAGAAGTTTTCAGAAGGTTTTACACTCTTATGGCCATCTg                      >  2:1447126/1‑150 (MQ=255)
                                                                                                                                     cAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTg                   >  1:425169/1‑150 (MQ=255)
                                                                                                                                      aaaCTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTgg                  >  1:84175/1‑150 (MQ=255)
                                                                                                                                      aaaCTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTgg                  >  2:841864/1‑150 (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)
                                                                                                                                           gtagtaCAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGATGGTACaa             >  2:126440/1‑150 (MQ=255)
                                                                                                                                             agtaGAATTCTCTCCAGTTTGAACAGGAAAGAATTTGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatt           <  1:722719/150‑1 (MQ=255)
                                                                                                                                             agtaCAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatt           >  2:385352/1‑150 (MQ=255)
                                                                                                                                              gtaCAATTCTCTCCATTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta          <  1:107747/150‑1 (MQ=255)
                                                                                                                                              gtaCAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAatta          >  1:1572842/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)
                                                                                                                                                     tctctcCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTa   <  1:1464131/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACAGGAAAGAATATGCTATGAACCCTTATATTTATCTTGGTGGTGCAATACTTGCAGAGGTCATTGGTACAACCTTAATGAAGTTTTCAGAAGGTTTTACACGGTTATGGCCATCTGTTGGTACAATTATTTGTTAT  >  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: