Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT41.1% 89.3 / 25.8 56N156K (AAC→AAAyiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (26/7);  new base T (0/23);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 7.85e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAG  >  NC_000913/4111597‑4111894
                                                                                                                                                         |                                                                                                                                                   
caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                                                                                                                     >  2:818785/1‑151 (MQ=255)
  aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgtt                                                                                                                                                   >  2:518464/1‑151 (MQ=255)
     gctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                                                                                                >  2:1055735/1‑151 (MQ=255)
        ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTc                                                                                                                                               <  2:778240/149‑1 (MQ=255)
          gctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccg                                                                                                                                           <  1:264728/151‑1 (MQ=255)
           ctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgc                                                                                                                                          >  2:1035389/1‑151 (MQ=255)
            tttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgcc                                                                                                                                         >  1:896637/1‑151 (MQ=255)
            tttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgcc                                                                                                                                         >  1:26636/1‑151 (MQ=255)
              tttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGt                                                                                                                                       >  1:811543/1‑151 (MQ=255)
                     tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGc                                                                                                                                >  1:718052/1‑151 (MQ=255)
                     tttATAAAATCCTCCAAGCTATTCAGATAGCGCCAGCTTAATCGGTTCAACAACGAGGGTCACCCCTTTTTCGTCTTTGCCCGCAACAACATACCGCAGCGCCCTTTCTGCCTCGTTGAAATAACGTTTTTTTTCCCCCGCCGTtagta                                                                                                                                  <  1:469264/149‑1 (MQ=255)
                      ttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCg                                                                                                                               >  2:395116/1‑151 (MQ=255)
                          tCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                                                                                           >  1:1063342/1‑151 (MQ=255)
                             aacttaaaaGCGATTCAGATAGCGCCAGCTTAATGGCTTCAACAGCGAAGGCCTGCCCCTTTCCGCGTTGGCCCGCAACAACATACCGCTGGGCACTTTCGGCCTGGGTGAAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAg                                                                                                                        <  2:904606/144‑1 (MQ=255)
                              aCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGt                                                                                                                       >  1:832211/1‑151 (MQ=255)
                               ccTTAAAGCGATTCAGATAGCGCCAGTTTAATCGGTTAAACAGCGAAGGTCAGCCCTTTTTCGCCGTTGCCCGCGAAAACATAACGCAGTGCCCCTTCTGTCGCGGTGTAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                      <  2:510854/151‑1 (MQ=255)
                                 ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                                                                                                    >  1:639852/1‑151 (MQ=255)
                                  tCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttttt                                                                                                                   >  2:1084719/1‑151 (MQ=255)
                                  tCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttttt                                                                                                                   >  2:853953/1‑151 (MQ=255)
                                   cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                                                                                  >  2:675336/1‑151 (MQ=255)
                                   cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                                                                                  >  1:127392/1‑151 (MQ=255)
                                    aatcgtttCAGATGGCCCCAGTTTATGGGTTCCAACAGCAAAGGCCGCCCCTTTTTCGCCTGTGTCGGCCACAACATAACGCAGTCCCCCTTCTGTCTCGGTTTAATATCGTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGttttt                                                                                                                   <  1:874670/143‑1 (MQ=255)
                                     aGCAATTCAGAGAGCGCCAGTTTAATCGGTTCAACAGCAAAGGCAACCCCTTTTTCGCCGTTGCCCGCAACAACAAACCGCGTTCCACCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                                                                                                <  2:844338/151‑1 (MQ=255)
                                       cGATTCAGATAGCGCCAGTTTAATCGTTTCAACAGCGAAGGTCACCCCTTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                                                                                              <  2:126919/151‑1 (MQ=255)
                                       cGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                                                                                              >  1:204060/1‑151 (MQ=255)
                                          ttCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGTTCACCCCCTTTTCGCCTTTGTCCGCGACAACATACCGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                            <  1:1070388/150‑1 (MQ=255)
                                             aaaTAGCCCCAGCTAAACCGGTCAAACAGCAAGGGTCACCCCTTTTTCGCCGTTGCCCGCAAAAACATACCGCAGTGCCCCTTCTGTCTCGGTTTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                                                                                         <  1:849306/148‑1 (MQ=255)
                                             aGATAGCTCCATCTTAATCGTTTCAACAGCAACGGCCACCCCCTTTTCGCCGTTGTCCCCAACAACACAACCCGGTCCCCCTTCTGTCTCGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                                                                                         <  2:1105104/150‑1 (MQ=255)
                                               aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTgtgt                                                                                                      >  2:366376/1‑151 (MQ=255)
                                                  gcgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                                                                                                   >  2:552592/1‑151 (MQ=255)
                                                   cgcCAGCTTAATCGGTTCAACAGCAAAGGCCAGCCCCTTTCCGCCTTTGTCCGCAACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                                  <  2:299974/151‑1 (MQ=255)
                                                     ccATCTTATTCGTTACAAAAGCGAAGGTCACCCCTTTGTCCCTGTTCCCCGCGAAAACAAACCGCAGTCCCCCTCTTGCCTTGGTGAATTAATGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                                                                                                    <  1:567392/147‑1 (MQ=255)
                                                     ccAGCTTAATCGGTTCAACAGCGAAGGTCACCCCCTTTCCGCCGTTGCCCGCGAAAACATAACGCAGTGCCCCTTTTGTCTGGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                                <  2:200797/151‑1 (MQ=255)
                                                       aGCTTAACCGTTAAAACAGCGAAGGCCGCCCCTTTTCCGCCTTTGTCCGCAACACCATAACGCATTCCCCTTTTTGTCTCGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                              <  1:243330/151‑1 (MQ=255)
                                                         cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                            >  2:906746/1‑151 (MQ=255)
                                                          ttATTCGGTTCAACGGCGAAGTTCACCCCTTTTTCGCCGTTGTCCGCGACAACATACCGCGGTGCCCCTTCTGTCTCGGTGTAATACCTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                              <  1:914842/148‑1 (MQ=255)
                                                          ttAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                            >  1:586812/1‑150 (MQ=255)
                                                          ttAATCGGTTAAACAGCAAAGGTCACCCCTTTTTCGCCTTTGCCCGCAACAACATACCGCAGTCCCCCTTCTGTCTCGGTTTAAAACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                              <  2:740067/148‑1 (MQ=255)
                                                           tAATCGGTTCACCAGTGATGGCCACCCTCTTTTCTCTTTTGTCCGCGCCAACATAACACAGTCCCCCTTCTGTCTGGGTTTAAAATCTTTTTTTTTCCTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                            <  1:275460/149‑1 (MQ=255)
                                                             accgtttCAACAGCGAGGGTCAGCCCCTTTCCGCCGTTGCCCGCGAAAACATAACGCATTCCCCCTTCTGTCTCGGTTTAATATCTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATg                                                                                        <  2:286002/146‑1 (MQ=255)
                                                                ggTTAAACAGCGAAGGTCACCCCTTTTCCGCCGTTGCCCGCGAAAACATACCGCGGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                                                                                       <  2:918563/149‑1 (MQ=255)
                                                                            gAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCgg                                                                          >  2:187334/1‑151 (MQ=255)
                                                                               ggTCACCCCCTTTCCGCCGTTGTCCGCAACAACATAACGCAGTCCCCCTTCTGTCTCGTTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                                       <  2:30964/151‑1 (MQ=255)
                                                                                   cagcccTTTTCCCCGTTGCGCCGATACACTATAACGCATTGCTCCTTCTGTCGGGGTGAAATGTGTTGTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                                       <  1:176100/145‑1 (MQ=255)
                                                                                     ccctTTTCCGCGTTTCTCCGCAACAACATACCGCGTTCCCCCTTTTGTCTGGGTTTAAAAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGc                                                                  <  2:213200/150‑1 (MQ=255)
                                                                                      ccttttcCGCCGTTGTCCGCGAAAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                                <  2:795137/148‑1 (MQ=255)
                                                                                                  gTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTcc                                                     >  2:469758/1‑151 (MQ=255)
                                                                                                  gTTGCCCCCGACAACAAACCGCAGTCCCCCTTCTTTCTCGGTTTAATATTTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTc                                                      <  1:586713/150‑1 (MQ=255)
                                                                                                        ccGCAACAACATACCGCTGTCCCCTTTCGGTCTGGGTGTAATAACGTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCt                                                    <  2:402746/147‑1 (MQ=255)
                                                                                                           cgACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgcttt                                             >  2:471713/1‑151 (MQ=255)
                                                                                                            aaaaacaTAAGGCGGTCCTCCTTTTGTTTCGGTATAAAATGTTTTTTTTTTCGCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctt                                              <  1:107505/146‑1 (MQ=255)
                                                                                                                 caTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                                       >  2:519561/1‑151 (MQ=255)
                                                                                                                       gCAGTGCCCTTTCTGTCTGGTTTTATTATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                                     <  1:156208/147‑1 (MQ=255)
                                                                                                                        cAGTGCCCCTTTTGTCTCGGTGTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCtttcgctttcgc                                    <  2:87595/147‑1 (MQ=255)
                                                                                                                                    tGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc                    >  1:206965/1‑151 (MQ=255)
                                                                                                                                         cGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg               >  1:569862/1‑151 (MQ=255)
                                                                                                                                            tgtaaaacttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCt                <  1:452731/137‑1 (MQ=255)
                                                                                                                                               aataatGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcccc                 <  1:757636/143‑1 (MQ=255)
                                                                                                                                                atatttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGt        <  1:62801/145‑1 (MQ=255)
                                                                                                                                                ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGt        >  1:391195/1‑151 (MQ=255)
                                                                                                                                                       ttgtttTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAg  <  1:667864/150‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAG  >  NC_000913/4111597‑4111894

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