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

AAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCA  >  NC_000913/4111598‑4111893
                                                                                                                                                         |                                                                                                                                                   
aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                                    <  2:576565/151‑1 (MQ=255)
aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                                    >  2:695189/1‑151 (MQ=255)
aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                                    >  1:246356/1‑151 (MQ=255)
 aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgtt                                                                                                                                                   >  2:662954/1‑151 (MQ=255)
 aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgtt                                                                                                                                                   >  2:890360/1‑151 (MQ=255)
  aGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgttt                                                                                                                                                  <  1:421398/151‑1 (MQ=255)
   ggCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                                 >  2:240756/1‑151 (MQ=255)
   ggCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                                 <  2:357377/151‑1 (MQ=255)
      tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTcc                                                                                                                                             >  2:448371/1‑151 (MQ=255)
          ctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgc                                                                                                                                         >  1:28794/1‑151 (MQ=255)
            ttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgcc                                                                                                                                        <  1:613067/150‑1 (MQ=255)
             tttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgccg                                                                                                                                       <  2:50320/150‑1 (MQ=255)
             tttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGt                                                                                                                                      >  2:451276/1‑151 (MQ=255)
                   aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagt                                                                                                                                  >  1:799634/1‑149 (MQ=255)
                      tATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGATAGTAGCg                                                                                                                              >  1:619052/1‑150 (MQ=255)
                           aaccctgaaagataTTCAAATAGCCCCAGTTTAATCGTTAAAACACAAACGGCCACCCCTTTTCCGCCTTGCCCCCCGACAACATACCGCAGTCCACTTCCTGCCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTg                                                                                                                          <  1:391078/136‑1 (MQ=255)
                              ccTTCAAGCGATTCAGATAGCGCCAGTTTAATCGTTTCAACAGCGAAGGTCACCCCCTTTCCGCCGTTGTCCGCAAAAACATAACGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                     <  2:200254/151‑1 (MQ=255)
                              ccTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                     >  2:157816/1‑151 (MQ=255)
                              ccTTCAAGCGATTAAGATACCGCCAGTTTAATCGTTACAACAGCGAAGCCCCCCCCTTTTTCGCTGTTGCCCGCAACAAAATAACGCAGTCCCCCTTCTGCCTCGGTTTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGt                                                                                                                      <  1:167435/150‑1 (MQ=255)
                              ccTACAAGCGATTCAGATACCCCCCTCTTAACCGGTTAAACAGCGAGGGCCAGCCCCTTTCCGCCGTTGCCCGCAAAAACATAACGCATTGCCCCTTCTGTCTCGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                     <  2:295921/151‑1 (MQ=255)
                               cTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGTGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttt                                                                                                                    >  2:83311/1‑151 (MQ=255)
                               cTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttt                                                                                                                    >  1:567707/1‑151 (MQ=255)
                                ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                                                                                                   >  1:500809/1‑151 (MQ=255)
                                 tCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttttt                                                                                                                  >  1:393352/1‑151 (MQ=255)
                                  cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                                                                                 >  2:6327/1‑151 (MQ=255)
                                   atgCGATTCAGACAGCCCCACTTTAATCGTTTAAACAGCGAAGGTCACCCCCTTTTCCCCGTTGTCCGCAAAAACAAACCGCTGTCCCCCTTTTGTCTCGGTTTAATAAGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGttt                                                                                                                    <  1:224194/145‑1 (MQ=255)
                                   aaGCGATTAAGATAGCCCCAGCTTAATCGGTTCAACAGCGAAGGTCACCCCCTTTCCGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                                                                                                   <  2:1016388/148‑1 (MQ=255)
                                   aaGCAATCCAGAAGGCCCCATTTTAATCGTTTCACCAGGGAGGGAACCCCCTTTTCCCCGGTTGCCCGCGAAACCAAACCGCAGTGCCCCTCTTTTCTGGGTTTAATAGCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                     <  2:347122/146‑1 (MQ=255)
                                       gATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                                                                                             >  1:629972/1‑150 (MQ=255)
                                         ttCAGATAGCCCCATTTAAATGGTTTAAACGGCGAAGCCCCCCCCCTTTCCGCCGTGGCCCGCAAAAACATAACGCATTCCCCCTTCTGTCTCGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                                                                                            <  1:942579/149‑1 (MQ=255)
                                           cAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                                                                                        >  2:815535/1‑151 (MQ=255)
                                            aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                                       >  2:563052/1‑151 (MQ=255)
                                            aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                                       >  2:856786/1‑151 (MQ=255)
                                               tAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtt                                                                                                    >  1:449702/1‑151 (MQ=255)
                                               tAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtt                                                                                                    >  2:949401/1‑151 (MQ=255)
                                                 gcgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                                                                                                  >  2:937590/1‑151 (MQ=255)
                                                  cgcCAGTTTAATCGTTTAAACAGCGAGGGTCACCCCTTTTCCGCGGTTGTCCGCAAAAAAATAACGCAGTGCACTTTCTGTCTCGGTGTAATATGTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                                 <  1:762482/151‑1 (MQ=255)
                                                   gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                                                                                                >  1:378542/1‑151 (MQ=255)
                                                    ccAGTTTAATCGTTCCAACAGCGAAGGTCACCCCCTTTTCGCCTTTGTCCGCAACAACATACCGCAGTGCCCCTTTTGTCTCGGTGTAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                                                                                                <  2:488803/150‑1 (MQ=255)
                                                      agtttAATCGGTCCAACAGCGAAGGCCACCCCTTTTTCGCCGTTGTCCGCAACAACATAACGCGGTCCCCTTTCTGTCTCGTTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTc                                                                                              <  2:493414/147‑1 (MQ=255)
                                                       gtttAATCGGTCCAACAGCGAAGTTCAGCCCTTTTCCGCCTTGTCCCGCGACACAAAAACGAGGTCCCCCTTCTGCCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                             <  2:651088/148‑1 (MQ=255)
                                                        cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                           >  2:695576/1‑151 (MQ=255)
                                                        cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                           >  2:702730/1‑151 (MQ=255)
                                                         ttAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                           >  2:817982/1‑150 (MQ=255)
                                                         taaaTCGGTTCAACAGCGAAGGTCCGCCTTTTTTCGCCGTTTCCCGCAACAACATACCGAAGTGCCCCTTCTGTCTGGGTGTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                                            <  1:756912/147‑1 (MQ=255)
                                                           attgggtcaaaCAGCAAGGGCCGGCCCCTTTCCGCTGTGGCCCGCGACAACATACCGCGGTCCGCCTTCTGTCTGGTTGTATTAGTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                                                                                        <  1:642188/142‑1 (MQ=255)
                                                           attCGTTTCAACAGCGAAGTCCACCCCCTTTCCGCCGTTGTCCGCAACAACAAACCGCAGTGCCCTTTTGGTCTCGGGGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                                                                                          <  2:1009797/147‑1 (MQ=255)
                                                                 ttCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTg                                                                                  >  1:750950/1‑151 (MQ=255)
                                                                   cAACAGCGAAGGCCACCCCCTTTTCGCCGTTGCCCGCGACAACATAACGCAGTCCCCCTTCTTTCTCGTTTTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                                   <  2:405008/148‑1 (MQ=255)
                                                                    aaCAGAGAAGGTCACCCCCTTTCCGCCGTTGCCCGCAAAAACATAACGCAGTGCACCTTTTGTCTCGGTGTAATACGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                                                                                 <  1:650285/149‑1 (MQ=255)
                                                                    aaCAGAGAAGGCCGCCCCTTTTCCCCCGTTTTCCGCAAAAAATTACCGCATTCCCCCTCCTTCCTTGTTAAAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAg                                                                                <  1:726400/150‑1 (MQ=255)
                                                                        gagagggTCCGCCCGTTTCCGCGTTTGTCCGCTACAACATAACGCAGTTCCCCTCGGTTCCGTGTGTAATAGCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtgt                                                                             <  1:510698/144‑1 (MQ=255)
                                                                        gCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc                                                                           >  2:868016/1‑151 (MQ=255)
                                                                         ggaaaggtgtccccTTTTCCGCCGTTTCCCGGGAAAACATACCGCGCTACTCTTTCTGCCTGGGTGTAATAGCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGtgtg                                                                            <  1:40279/140‑1 (MQ=255)
                                                                          gAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCgg                                                                         >  2:764892/1‑151 (MQ=255)
                                                                            agggccaccccTTTTCCGCGGTGGCCCGCGAAAAAATACCGCGTTCCCCTTTCTGTCTGGGTTTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGa                                                                        <  1:422199/144‑1 (MQ=255)
                                                                                 aGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                                  >  1:937560/1‑151 (MQ=255)
                                                                                   ccccTTTTCGCCGTTGTCCGCGAAAACAAAACGCAGTCCCCCTTCTGTCTGGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                                <  2:51378/151‑1 (MQ=255)
                                                                                     ccTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGc                                                              >  2:354564/1‑151 (MQ=255)
                                                                                      cTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                             >  1:451556/1‑151 (MQ=255)
                                                                                            cGCCTTTGTCCGCGACAACATACCGCTGTGCACTTTTTCTCTCGTTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                          <  2:629394/149‑1 (MQ=255)
                                                                                              ccGTTGTCCACGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTAtt                                                      >  2:974547/1‑151 (MQ=255)
                                                                                              ccGTTGCCCGCAACAACATAACGCTGTCCCCCTTCTGTCTCGGTTTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTAt                                                       <  1:240753/150‑1 (MQ=255)
                                                                                                       ccGCGACAACATACCGCAGTGCCCCTTCTGTCTCGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                                 <  2:473445/149‑1 (MQ=255)
                                                                                                       ccGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATATCTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAg                                                <  2:769674/150‑1 (MQ=255)
                                                                                                         gcgACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctt                                             >  1:766010/1‑151 (MQ=255)
                                                                                                          cgAAAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgct                                              <  2:114974/149‑1 (MQ=255)
                                                                                                                 aTAGCGCATTCCCCCTTCTGTCTCGGTGTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgct                                        <  1:718852/148‑1 (MQ=255)
                                                                                                                     cgctgtcccccTTCTGTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                                  <  1:376477/141‑1 (MQ=255)
                                                                                                                     cGCAGTGCCCTTTTTGTCTCGGTGTAATATTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCgg                                 <  1:312215/151‑1 (MQ=255)
                                                                                                                            cACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTg                          >  1:306379/1‑151 (MQ=255)
                                                                                                                             aCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGc                         >  1:634409/1‑151 (MQ=255)
                                                                                                                              ccTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                        >  1:8735/1‑151 (MQ=255)
                                                                                                                                    gccgcGGTATAATACCGTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc                   <  1:971098/146‑1 (MQ=255)
                                                                                                                                     tctcGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttttt                           <  2:446794/141‑1 (MQ=255)
                                                                                                                                       tcGGTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc                  <  1:215018/148‑1 (MQ=255)
                                                                                                                                        tggttgtAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTg                     <  2:713633/141‑1 (MQ=255)
                                                                                                                                        cGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg              >  1:111518/1‑151 (MQ=255)
                                                                                                                                         ggTGTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCt               <  2:435247/149‑1 (MQ=255)
                                                                                                                                         ggTGTAATACCGTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg              <  2:1018211/150‑1 (MQ=255)
                                                                                                                                            tgtAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg                    <  2:498478/142‑1 (MQ=255)
                                                                                                                                            tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAt           >  1:168910/1‑151 (MQ=255)
                                                                                                                                            tataaaatttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcccc                <  1:638875/136‑1 (MQ=255)
                                                                                                                                                ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATggg        >  2:638247/1‑150 (MQ=255)
                                                                                                                                                  aagtttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAg     <  1:958763/147‑1 (MQ=255)
                                                                                                                                                  aaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGt       >  2:350135/1‑149 (MQ=255)
                                                                                                                                                    cGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGcc   >  2:940319/1‑151 (MQ=255)
                                                                                                                                                     gtttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGGGGGTGTGCGGATCACCGCTGCCTTTTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGcc   >  1:873715/1‑150 (MQ=255)
                                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCa  <  2:567720/150‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
AAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCA  >  NC_000913/4111598‑4111893

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