Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,727 G→T 14.7% intergenic (‑282/‑107) ylcI ← / → nohD DUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,7270GT14.7% 467.5 / 19.1 180intergenic (‑282/‑107)ylcI/nohDDUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (81/65);  new base T (3/31);  total (84/96)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

GTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580579‑580876
                                                                                                                                                    |                                                                                                                                                     
gTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGt                                                                                                                                                      <  2:1739658/150‑1 (MQ=255)
  tgtgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtat                                                                                                                                                    <  1:2470514/150‑1 (MQ=255)
    tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattt                                                                                                                                                  >  1:1107439/1‑150 (MQ=255)
          gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCGGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            >  1:3885258/1‑150 (MQ=255)
          gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            >  1:2717820/1‑150 (MQ=255)
                 atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                     <  1:941660/150‑1 (MQ=255)
                   aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                   >  1:1952020/1‑150 (MQ=255)
                   aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccccc                                                                                                                                   >  1:1960532/1‑150 (MQ=255)
                       attatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaa                                                                                                                               <  2:3422429/150‑1 (MQ=255)
                           tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                           >  1:3661483/1‑150 (MQ=255)
                           tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                           >  1:441600/1‑150 (MQ=255)
                           tCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaaga                                                                                                                           <  1:2449342/150‑1 (MQ=255)
                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:2958733/150‑1 (MQ=255)
                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:3165072/150‑1 (MQ=255)
                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:3929694/150‑1 (MQ=255)
                              ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                        <  1:1422433/150‑1 (MQ=255)
                                  gtgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                    <  1:1034838/150‑1 (MQ=255)
                                  gtgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                    >  1:1814571/1‑150 (MQ=40)
                                  gtgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTGTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                    >  2:2288173/1‑150 (MQ=255)
                                    gggTCCTTTCCGGCGATCGGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTTAGGTGTTTGCGTTCTTCTTCGTCGTAACTTATTGTATTTATTTTATATACGCCCTGTAAAGAAAGGATACg                                                                                                                  >  2:416074/1‑150 (MQ=14)
                                     ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 >  2:3065671/1‑150 (MQ=39)
                                     ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 <  2:3276583/150‑1 (MQ=39)
                                     ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 <  2:973328/150‑1 (MQ=39)
                                      gTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGAc                                                                                                                >  2:2491317/1‑150 (MQ=39)
                                      gTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGAc                                                                                                                >  2:2318184/1‑150 (MQ=39)
                                        ccTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAg                                                                                                              >  2:822311/1‑150 (MQ=39)
                                        ccTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAg                                                                                                              >  2:2198516/1‑150 (MQ=39)
                                           ttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTg                                                                                                           >  2:1362462/1‑150 (MQ=39)
                                           ttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTg                                                                                                           >  2:3869326/1‑150 (MQ=39)
                                             ccGGCTATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCt                                                                                                         <  2:277860/150‑1 (MQ=255)
                                             ccGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCt                                                                                                         >  2:2623909/1‑150 (MQ=39)
                                               ggCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                       >  1:915439/1‑150 (MQ=39)
                                                 cGATCCGACTGGTTAGGGGGCGGCGACCTCGCGGTTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTTTTTCCTTTCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaaa                                                                                                     <  2:2378724/150‑1 (MQ=11)
                                                 cGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaaa                                                                                                     <  2:2915393/150‑1 (MQ=255)
                                                  gATCCGACAGGTTCGGGGGCGGCGACCTCGGGGGTTTTCGCTATTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTTTTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAg                                                                                                    <  2:1542575/150‑1 (MQ=11)
                                                     ccGACTGGTTAGGGGGCGGCGACCTCGGTGTTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCTTAACTTAATTTTTTTTTTTAAAAAACCCCCTGAAAAGAAAGCAAACGACAGGTGCTGAAAGCGa                                                                                                 <  2:1446435/150‑1 (MQ=11)
                                                       gACAGGTTCCGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTATGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               <  2:650299/150‑1 (MQ=255)
                                                       gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:329304/1‑150 (MQ=39)
                                                       gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               <  2:1267613/150‑1 (MQ=39)
                                                       gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:2472132/1‑150 (MQ=39)
                                                       gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:2472133/1‑150 (MQ=39)
                                                          aGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttt                                                                                            >  2:599219/1‑150 (MQ=39)
                                                            gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          <  2:2156843/150‑1 (MQ=40)
                                                            gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          >  2:3653823/1‑150 (MQ=40)
                                                               tcGGGGCGGCGACCTCGGGTGTTTTCTCTATTTATTAAAATTTTCGTGTTTAAGGTTTTACCGTTCTTCTTCGTCGTAAATTATTGTATTTATTTAAAATACCCCCTGAAAAGAATGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGc                                                                                       <  2:1807938/149‑1 (MQ=14)
                                                                cGGGGCGGCGACCTCGGGGGTGTTCGCTATTTATGAAATTTTTCCGGTTTATGGTGTTTCCGTTCTTCTTCGTCGTAACTTATTTTTTTTATTTAAAATACCCCCTGAAACGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGcc                                                                                      <  2:1144358/150‑1 (MQ=11)
                                                                 ggggCGCCGACTTCGGGTTTTTTAGCTTTTTATCAAAATTTTCCGTTTTTTGGTTTTTCCTTTATTCTTCTTCGTAAATTTTTTTTTTTATTAAAACAACCCCCTTAAAATAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCct                                                                                     <  1:2317619/150‑1 (MQ=255)
                                                                    gcggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAATGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTg                                                                                  >  1:2872633/1‑150 (MQ=255)
                                                                    gcggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTg                                                                                  >  2:847261/1‑150 (MQ=255)
                                                                      ggcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTATGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTc                                                                                <  2:2454434/150‑1 (MQ=255)
                                                                        cgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGt                                                                              <  2:1996190/150‑1 (MQ=255)
                                                                          cccTCGCGGGTTTTCGCTATTTCTGACAATTTTCCGGTTTACGGTGTTTCCGTTCTTCTTCGTCGTAACTTACTGTATTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGttt                                                                            <  2:3195578/149‑1 (MQ=14)
                                                                          aCCTCGCGTGTTTTCGCTATTTATGAAATTTTTCCGGTGTACGGCGTTACCTTTCTTCTACGCCGTAACTAAATTTATTAATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGttt                                                                            <  2:3335841/150‑1 (MQ=11)
                                                                          aCCTCGCGGGTTTTCGCTATTTATTAAAATTTTACGGTTTAAGGTGTTTCCGTTCTTCTTAGTCGTAACTTAATGTATTTATTTAAAATACCCCCAGAAAAGAAAGGAAACGACAGGTGCTGAAATCGAGCTTTTTGGCCTCTGTCGttt                                                                            >  1:3748446/1‑150 (MQ=14)
                                                                          aCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGttt                                                                            <  2:1957775/150‑1 (MQ=255)
                                                                          aCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTACGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCAGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCATTTTGGCCTCTGTCGttt                                                                            >  1:1336308/1‑150 (MQ=14)
                                                                            cTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCGGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTcc                                                                          >  1:3273065/1‑150 (MQ=255)
                                                                               ggggtttttCGATATTTTTGAAAATTTTCCGGTTTATGGTTTTTCCGTTCTTCTTCGTCTTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCttt                                                                       <  2:45238/145‑1 (MQ=14)
                                                                               gcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCttt                                                                       <  2:2863506/150‑1 (MQ=255)
                                                                                cgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGACGTTTCCTTtc                                                                      >  2:666821/1‑150 (MQ=255)
                                                                                 gggTTTTCGCTATTTATGAAAATTTTCCGGTTTCAGGTGTTTCCGTTCTTCTTCGTCGTAACTTATTGTATTTATTTAAAATACCCCATGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     <  2:406670/150‑1 (MQ=255)
                                                                                   gTTTTCTATTTTTTTTAAAATTTTCCGTTTTATGGTGTTTCCTTTCTTCTTCGTCTTAACTTAATTTATTTATTTAAAATACCCCCTTAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTctct                                                                   <  2:1364339/150‑1 (MQ=14)
                                                                                    ttttCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTg                                                                  >  1:408155/1‑150 (MQ=255)
                                                                                       tcgatttttATTAAATTTTTCGGTTTTATGGTTTTTCCTTTCTTATTCTTCTAAATTTAATGTATTTATTTAAAATCCCCCCTCAAAACAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGttt                                                               <  1:540472/144‑1 (MQ=14)
                                                                                           tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           >  2:3719807/1‑150 (MQ=255)
                                                                                           tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           >  1:1672752/1‑150 (MQ=255)
                                                                                           tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           >  2:2671563/1‑150 (MQ=255)
                                                                                           tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           >  1:1189718/1‑150 (MQ=255)
                                                                                           tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           >  2:2590872/1‑150 (MQ=255)
                                                                                           tATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTACGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                           >  2:3020966/1‑150 (MQ=255)
                                                                                            aTTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTc                                                          >  1:2175434/1‑150 (MQ=255)
                                                                                                  gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACACCCTGAAAAGAAAGGAAACGACAGGTGCTGAAATCGAGCTTTATGGCCTCTGTCGTTTCCTTTCACTGTTTTTGTCCGTGGa                                                    >  1:2599825/1‑150 (MQ=255)
                                                                                                   aaaaTTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTGGTAACTTAATTTATTTATTTAAAATACCCCCTGAAAATAATGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                                   <  2:163460/150‑1 (MQ=255)
                                                                                                    aaaTTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCGGTTTTTGTCCGTGGAAt                                                  >  1:1580102/1‑150 (MQ=255)
                                                                                                      aTTTTCCGGTTTATGGTGTTTCCGTTATTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAATAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgt                                                <  2:1325944/150‑1 (MQ=255)
                                                                                                        tttCCGGTTTCTGGGTTTACCTTTCTCCTTCGTCTAAACTTATTGTATTTCTTTAAACACCCCCCTTAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGc                                              <  1:680244/150‑1 (MQ=14)
                                                                                                        tttCCGGTTCTTGGTGTTTCCGTTCTTCTTCTGCTTAATTTATTTTTTTTATTTAAAATACCCCCTCAAAACAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGc                                              <  1:680184/150‑1 (MQ=14)
                                                                                                          tCCGGTTTTTGTTGTTTCTGTTCTTCTTCTTCTTAAATTATTTTTTTTATTTAAAATACCCCATCAAAACAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                            <  1:3301042/150‑1 (MQ=14)
                                                                                                            gggTTTATGGTGTTTCCGTTCTTCTTCTTCTTAAATTATTTTATTTTTTTAAAATACCCCCTGAAAAGAATGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATg                                          <  1:2539719/149‑1 (MQ=17)
                                                                                                             ggTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATgg                                         >  2:1383495/1‑150 (MQ=255)
                                                                                                               tttATGTTGTTTCCGTTCTTATTCGTCGTAACTTATTTTATTTATTAAAAATACCCCCTCAAAACAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  1:3237752/150‑1 (MQ=14)
                                                                                                               tttAAGGTGTTTCCTTTCTTCTTCGTCGTATCTTATTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGaa                                       <  2:2767041/150‑1 (MQ=255)
                                                                                                                 tAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATAACCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCATTATCTGTTTTTGTCCGTGGAATGTGCAATGGAAGt                                     >  1:3287571/1‑150 (MQ=255)
                                                                                                                 tAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGt                                     >  2:2861745/1‑150 (MQ=255)
                                                                                                                 tAAGGTGTTTCCGTTCTTATTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGt                                     <  2:3281262/150‑1 (MQ=255)
                                                                                                                   aGGTGTTTCCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTca                                   <  2:883920/150‑1 (MQ=255)
                                                                                                                     gtgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAAACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTGTTTTTTTCCGTGGAATGTGAAAATGAAGTcaac                                 >  1:1551584/1‑150 (MQ=14)
                                                                                                                        tttCCGTTCTTCTTGTTCGTAAATTAATTTATTTATTTAAAATACCCCCTTAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaa                              <  1:2491291/150‑1 (MQ=14)
                                                                                                                             gTTCTTCTTCGTCGTAACTTAAAGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                         >  2:1095551/1‑150 (MQ=255)
                                                                                                                                    ttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                  <  1:2623909/150‑1 (MQ=255)
                                                                                                                                    ttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                  >  2:574671/1‑150 (MQ=255)
                                                                                                                                      cgcctaaaCTTAATTTATTTATTTAAAACACCCCCTTAAAATAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGa                <  1:3471517/144‑1 (MQ=14)
                                                                                                                                       gtcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGAc               >  2:2869206/1‑150 (MQ=255)
                                                                                                                                                tAATGTATTTATTTAAAACCCCCCCTGAAAACAATGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGt      <  2:3675997/150‑1 (MQ=14)
                                                                                                                                                  aTGTTTTTATTTAAAATACCCCCTTAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGc    <  1:3623994/150‑1 (MQ=255)
                                                                                                                                                  aTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGc    <  2:441600/150‑1 (MQ=255)
                                                                                                                                                   tGTATTTATTTAAAATACCCCCTGAAAATAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  <  2:357966/150‑1 (MQ=255)
                                                                                                                                                    gTATTTATTTAAAATACCCCCTGAAAATAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  <  1:2617378/150‑1 (MQ=255)
                                                                                                                                                    gTATTTATTTAAAATACCCCCTGAAAAGACAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  <  2:1952020/150‑1 (MQ=255)
                                                                                                                                                    gTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGa  <  1:3256653/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580579‑580876

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: