Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC43.1% 94.8 / 11.5 58intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (28/5);  new base C (0/25);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 8.35e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATC  >  NC_000913/2001769‑2002062
                                                                                                                                                     |                                                                                                                                                
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                 <  2:430420/151‑1 (MQ=255)
     tataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                            >  2:399346/1‑151 (MQ=255)
        agggTGAATTCACGATAAACAGCCCTGCTTTATATGAGTTATCGGCATGATTTTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgccg                                                                                                                                          <  1:982716/147‑1 (MQ=255)
              aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGGGAAAGGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCgg                                                                                                                                   <  2:1005130/151‑1 (MQ=255)
                  cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                >  1:447625/1‑150 (MQ=255)
                   aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGGAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGa                                                                                                                                <  2:672331/149‑1 (MQ=255)
                    cGATAAACACCCCTGCTTTATATGATTTATCGGCATGATTATCCTTTTCTGCGGGTTTTTTAACCGGACGATTGGTGGGTGAAATGGGGGTTTTTTTGGGGTTTACAGTTAATTCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                             <  1:881119/151‑1 (MQ=255)
                      aaaaaCAGCCCTGCGTTATATGAGTTATCGGCATGTTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTGGGGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                           <  1:65058/149‑1 (MQ=255)
                       tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg                                                                                                                          >  1:460378/1‑151 (MQ=255)
                        aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                         >  2:424533/1‑151 (MQ=255)
                         aaCAGCCCTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTAAATCGGACGTTTAGGGGGTAAAATGAGGGTTTTTTTGGGGGTTACGGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                        <  1:261576/151‑1 (MQ=255)
                          aCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                                       >  2:947279/1‑151 (MQ=255)
                           cAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                                      >  1:164237/1‑151 (MQ=255)
                              cccTGCTTTATATGATTTATCGGCATGTTTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                   <  1:97470/151‑1 (MQ=255)
                              cccGGCGTTATATGAGTTAGCGGCATGATTATCCTTTTCTGAAGGGTTTTAAACCGGACGATTGGTGGGTGAAATGAGGGGTTATTGGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                    <  2:902798/150‑1 (MQ=255)
                               ccTGCGTTATATGAGTTAGCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                  <  1:693058/151‑1 (MQ=255)
                                 tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                 >  1:953258/1‑150 (MQ=255)
                                     ttttatGGGTTTTCGGCATGTTTATCCGTTTCTGCGGGTTTTTTAATCGGAGGTTTGGTGGTTAAAATGGGGGGTTTTTTGGGGTTTAAAGTTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                                <  1:361306/144‑1 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGTTTATTGGGGGTTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                               <  2:342981/148‑1 (MQ=255)
                                      tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                             >  1:832/1‑149 (MQ=255)
                                            aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                      >  2:224318/1‑150 (MQ=255)
                                               tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                  >  1:189047/1‑151 (MQ=255)
                                                aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                 >  2:853777/1‑151 (MQ=255)
                                                aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                 >  1:963595/1‑151 (MQ=255)
                                                    gCATGATTTACCGTTTCCGGGAGGTTTTTAATCGGAAGTTTGGGGGGGGAAAAGGGGGTTTTTTTTGGGGTTAAAGGGAAATTCCCAGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                             <  1:588246/151‑1 (MQ=255)
                                                     taTGATTATATGTTTTTGTAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTAAAGGTAAATTATAGGCAGAAAAAAACGACGACGGTGGCGGGGAAGCACGTTGCTGAGAAATTGGGCTTTATGTTgtagg                                                                                            >  1:96058/2‑147 (MQ=255)
                                                          ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                       >  2:831495/1‑151 (MQ=255)
                                                           tgccctttttggcggggTTTTTAACCGAACGATAGGGGGGGAAAAGGAGGGTTTTTTTGGGGGTTAGAGGAAAATACCAGCAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                          <  1:471838/133‑1 (MQ=255)
                                                                   tCTGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                               <  1:86614/150‑1 (MQ=255)
                                                                   tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                              >  1:43349/1‑151 (MQ=255)
                                                                         ggttttttAACCGGACTATTAGGGGGTGAAAGGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                          <  2:545470/146‑1 (MQ=255)
                                                                         gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                        >  1:292628/1‑151 (MQ=255)
                                                                            tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                     >  2:760871/1‑151 (MQ=255)
                                                                              tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                   >  2:409286/1‑151 (MQ=255)
                                                                                tAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTTGGGGTTTACAGTTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                  <  2:878567/150‑1 (MQ=255)
                                                                                 aaTGGGACGTTTGGTGGGTAAAATGAGGGTTTTTTTGGGGGTTACGGTTAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                 <  1:147818/150‑1 (MQ=37)
                                                                                     ggACGATTGGTGGGTGAAATGAGGGTTTATTGGGGGGTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                            <  2:711043/151‑1 (MQ=255)
                                                                                        cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                         >  1:810760/1‑151 (MQ=255)
                                                                                         gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                        >  2:50124/1‑151 (MQ=255)
                                                                                             aGTGGGTGAAATGAGGGGATATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                          >  1:611219/1‑145 (MQ=255)
                                                                                               tGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                   >  2:320107/1‑150 (MQ=255)
                                                                                                gggTGAAATGAGGGTTTATTGGGGGTTTACGGGTAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                    <  2:706576/148‑1 (MQ=38)
                                                                                                gggTGAAATGAGGGTTTATTGGGGGTTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                      <  2:963693/146‑1 (MQ=255)
                                                                                                     aaaTGGGGGGTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                            <  1:399317/151‑1 (MQ=255)
                                                                                                        tGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                          >  1:1017563/1‑150 (MQ=255)
                                                                                                           gggtTTTTTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                          <  1:266027/147‑1 (MQ=255)
                                                                                                           gggtTTATTTGGGGGTTACAGTAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATACGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                      <  1:687961/151‑1 (MQ=37)
                                                                                                           ggggTTATTGGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                      <  2:951561/151‑1 (MQ=255)
                                                                                                            gggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                     <  2:22806/151‑1 (MQ=255)
                                                                                                              gTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                    >  2:966363/1‑150 (MQ=255)
                                                                                                                 aTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                >  1:1020956/1‑151 (MQ=255)
                                                                                                                  tttGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAACTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                <  1:858296/150‑1 (MQ=255)
                                                                                                                            aCAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                     >  2:925299/1‑151 (MQ=255)
                                                                                                                                  aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                >  2:734850/1‑150 (MQ=255)
                                                                                                                                       ccAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg           <  1:229886/150‑1 (MQ=255)
                                                                                                                                       ccAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg           <  2:281511/150‑1 (MQ=255)
                                                                                                                                       ccAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc            <  1:119853/149‑1 (MQ=255)
                                                                                                                                              gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt   >  1:961664/1‑151 (MQ=255)
                                                                                                                                               aaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt   <  2:296776/150‑1 (MQ=255)
                                                                                                                                               aaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc  <  2:864945/151‑1 (MQ=255)
                                                                                                                                                aaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc  >  1:841712/1‑150 (MQ=255)
                                                                                                                                                 aaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc  >  1:663970/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATC  >  NC_000913/2001769‑2002062

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑