Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,4980GA57.4% 62.6 / 134.2 101V148I (GTT→ATT) nohBbacteriophage DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (29/29);  total (51/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 109 total aligned reads displayed.

GGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAG  >  NC_000913/580415‑580593
                                                                                   |                                                                                               
ggctCGTAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGa                                                                                     <  1:6176300/93‑1 (MQ=0)
       agttGAAATTGCCACTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGa                                                                              <  1:1176939/93‑1 (MQ=0)
                                   atgccccTGTCGGTGCAGAGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGcc                                                  >  1:8389948/4‑96 (MQ=0)
                                      cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCtgcg                                               >  2:14306642/4‑96 (MQ=0)
                                      cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCtgcg                                               >  1:5898286/4‑96 (MQ=0)
                                                          cttcccGGAGCTGGAAACCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCggg                           <  1:2936405/92‑1 (MQ=0)
                                                          cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCggg                           >  2:9622432/5‑96 (MQ=0)
                                                          cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCggg                           <  2:12302344/92‑1 (MQ=0)
                                                          cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCggg                           >  2:5714961/5‑96 (MQ=0)
                                                          cttcccGGAGCTGGAAAACCGACATATGGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCggg                           >  2:10836076/5‑96 (MQ=0)
                                                           ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCgggg                          <  2:10490969/93‑1 (MQ=0)
                                                           ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCgggg                          <  1:4854400/93‑1 (MQ=0)
                                                           ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCgggg                          >  2:12659087/4‑96 (MQ=0)
                                                           ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCgggc                          <  1:11271892/93‑2 (MQ=0)
                                                                ggcaCTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCt                     <  1:5522607/93‑1 (MQ=1)
                                                                 gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGGTGCTg                    >  2:616927/4‑96 (MQ=0)
                                                                 gaccTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCtt                    <  1:11328336/93‑2 (MQ=1)
                                                                            ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat         >  1:6178562/1‑96 (MQ=31)
                                                                            ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat         <  1:9253283/96‑1 (MQ=31)
                                                                            ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat         <  2:4501099/96‑1 (MQ=31)
                                                                            ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat         <  1:14017982/96‑1 (MQ=31)
                                                                            ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat         >  2:12059625/1‑96 (MQ=11)
                                                                            ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat         <  2:11846428/96‑1 (MQ=11)
                                                                            ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat         <  2:10586886/96‑1 (MQ=11)
                                                                            ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat         >  2:13573492/1‑96 (MQ=11)
                                                                            ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat         <  1:3021734/96‑1 (MQ=11)
                                                                             cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAATCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc        >  2:10300621/1‑96 (MQ=17)
                                                                             cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc        >  1:290336/1‑96 (MQ=32)
                                                                             cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc        <  1:5300820/96‑1 (MQ=32)
                                                                             cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc        <  1:8310594/96‑1 (MQ=32)
                                                                             cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc        >  2:1214579/1‑96 (MQ=32)
                                                                             cGACATATTGATTTCCTCATGAAGGTCATCATTAAAGCCATGATCAAAGCAGCTGCGCTGGCTGCAATAATCCCGGGGTTGCTGAGTGAatatatn        <  1:1205098/96‑2 (MQ=2)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGTAATAATACCGGGGTTGCTGAGTGAATATATc        >  2:4632032/1‑96 (MQ=11)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc        >  1:10856790/1‑96 (MQ=11)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc        >  2:9616420/1‑96 (MQ=11)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc        >  2:4128964/1‑96 (MQ=11)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc        <  2:8955439/96‑1 (MQ=11)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc        >  2:5642696/1‑96 (MQ=11)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc        <  1:5598898/96‑1 (MQ=11)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc        <  2:2625773/96‑1 (MQ=11)
                                                                             cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc        <  2:3490362/96‑1 (MQ=11)
                                                                              ggcATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGCTGCTGAGTGAATATATCg       >  1:11881950/3‑96 (MQ=11)
                                                                              gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg       <  2:8970517/96‑1 (MQ=33)
                                                                              gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg       >  2:10617778/1‑96 (MQ=33)
                                                                              gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg       <  1:8909584/96‑1 (MQ=33)
                                                                              gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg       >  1:8774641/1‑96 (MQ=33)
                                                                              gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg       >  1:4776672/1‑96 (MQ=33)
                                                                              gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGATTATATCg       >  1:4326504/1‑96 (MQ=16)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       <  2:12998471/96‑1 (MQ=11)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       <  1:14334549/96‑1 (MQ=11)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       <  2:5870438/96‑1 (MQ=11)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       >  2:1685750/1‑96 (MQ=11)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       >  2:7208705/1‑96 (MQ=2)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       <  1:6203468/96‑1 (MQ=11)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       <  1:8175934/96‑1 (MQ=11)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       <  2:9774323/96‑1 (MQ=11)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg       <  1:8626059/96‑1 (MQ=11)
                                                                              gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAAGAATACCGGGGGTGCTGAGTGAATATATCg       >  1:10411990/1‑96 (MQ=2)
                                                                               aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa      <  2:9653034/96‑1 (MQ=34)
                                                                               aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATAGATCGa      >  2:12738204/1‑96 (MQ=18)
                                                                               aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa      <  1:11794217/96‑1 (MQ=11)
                                                                               aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa      >  1:8923132/1‑96 (MQ=2)
                                                                               aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa      >  2:5649533/1‑96 (MQ=2)
                                                                               aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa      >  2:336372/1‑96 (MQ=2)
                                                                               aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa      <  2:1741192/96‑1 (MQ=11)
                                                                               aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa      <  2:12801911/96‑1 (MQ=11)
                                                                               aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa      >  2:1225152/1‑96 (MQ=11)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCGGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa     >  2:6773774/1‑96 (MQ=25)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa     >  2:13520927/1‑96 (MQ=35)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa     >  2:706005/1‑96 (MQ=35)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa     <  2:13045783/96‑1 (MQ=35)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa     >  2:6586974/1‑96 (MQ=35)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa     <  2:9404469/96‑1 (MQ=35)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa     >  2:9826294/1‑96 (MQ=35)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa     <  1:1748369/96‑1 (MQ=35)
                                                                                cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCGGAGTGAATATATCGaa     >  1:12807879/1‑96 (MQ=25)
                                                                                cATATTGATTTCGTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     >  1:7373227/1‑96 (MQ=11)
                                                                                cATATTGATTTCCTCAAGAAGGATATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     <  2:13298159/96‑1 (MQ=11)
                                                                                cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     <  2:706063/96‑1 (MQ=11)
                                                                                cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     >  1:9847076/1‑96 (MQ=2)
                                                                                cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     >  1:12213835/1‑96 (MQ=11)
                                                                                cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     >  2:200107/1‑96 (MQ=2)
                                                                                cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     <  1:5543789/96‑1 (MQ=11)
                                                                                cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     >  2:11214691/1‑96 (MQ=2)
                                                                                cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa     <  1:5981300/96‑1 (MQ=11)
                                                                                 aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCTCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc    <  2:12865006/96‑1 (MQ=37)
                                                                                 aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc    >  1:3782442/1‑96 (MQ=37)
                                                                                 aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc    <  2:8758841/96‑1 (MQ=37)
                                                                                 aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAAc    >  2:14283961/1‑96 (MQ=37)
                                                                                 aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAAc    >  2:8727217/1‑96 (MQ=37)
                                                                                  tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa   <  1:5560989/96‑1 (MQ=38)
                                                                                  tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa   <  2:8949719/96‑1 (MQ=38)
                                                                                  tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa   <  1:12314381/96‑1 (MQ=38)
                                                                                  tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa   >  1:7150515/1‑96 (MQ=38)
                                                                                   gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg  <  2:10001663/96‑1 (MQ=39)
                                                                                   gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg  <  1:8666682/96‑1 (MQ=39)
                                                                                   gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg  >  1:11034744/1‑96 (MQ=39)
                                                                                   gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg  >  2:8079739/1‑96 (MQ=39)
                                                                                   gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg  <  1:5243417/96‑1 (MQ=39)
                                                                                   gTTGATTTCATGAAACGGGAAATCATCAACGCCATCAACAAAGCAGCCGAGCTGGATGAACTGATCCCGGGGTTGCTGAGTGAATATATCGAACAg  <  1:4111865/96‑1 (MQ=14)
                                                                                   |                                                                                               
GGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAG  >  NC_000913/580415‑580593

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: