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

CAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAG  >  NC_000913/4111597‑4111894
                                                                                                                                                       |                                                                                                                                                    
caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                                                                                                                      >  2:849243/1‑151 (MQ=255)
 aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                                                                                                     >  1:917083/1‑151 (MQ=255)
   aGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgttt                                                                                                                                                   >  2:336768/1‑151 (MQ=255)
   aGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgttt                                                                                                                                                   <  2:672052/151‑1 (MQ=255)
     gctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                                                                                                 >  2:92376/1‑151 (MQ=255)
     gctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                                                                                                 <  1:937090/151‑1 (MQ=255)
       tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTcc                                                                                                                                              >  2:252639/1‑151 (MQ=255)
       tttgctttTTGTCACTTATTGAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCACCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTcc                                                                                                                                              >  1:434920/1‑151 (MQ=255)
        ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTcc                                                                                                                                              <  1:1037521/150‑1 (MQ=255)
          gctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccg                                                                                                                                           <  2:291698/151‑1 (MQ=255)
            tttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgcc                                                                                                                                         >  1:2150/1‑151 (MQ=255)
             ttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgcc                                                                                                                                         <  1:74416/150‑1 (MQ=255)
             ttttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCccgc                                                                                                                                          <  1:1058655/149‑1 (MQ=255)
              tttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAACCGTTTCAACAGCGAAGGCCAGCCCCTTTTCGCCGTTGTCCGCGACAACATACCGCAGTGCACCTTCTGTCCCGGTGTAATAACGTTTGTTTTCCCccgccg                                                                                                                                        <  1:316603/150‑1 (MQ=255)
                     tttATTAAAACCTTCAAGCGATTCAGATAGCCCCAGCTTAATGGGTACAACAGCGAAGGTCACCCCCTTTTCGCTTTTTTCCGCGACAACATAACGCAGTGCCCCTTCTGTCGCGGTGTAATACTTTTTTTTTTTCCCCGCCGTtagtag                                                                                                                                 <  1:471264/150‑1 (MQ=255)
                      ttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCg                                                                                                                               >  2:276892/1‑151 (MQ=255)
                      ttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCg                                                                                                                               >  2:815677/1‑151 (MQ=255)
                        aTTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCTAc                                                                                                                             >  2:685696/1‑151 (MQ=255)
                          tCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                                                                                           >  2:139390/1‑151 (MQ=255)
                           cAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                                                                                          >  1:422816/1‑151 (MQ=255)
                           cAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                                                                                                          >  2:260534/1‑151 (MQ=255)
                            aacccTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTACAACAGCGAAGGCCACCCCCTTTTCGCCGTTGTCCGCGAAAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCa                                                                                                                         <  2:316031/148‑1 (MQ=255)
                               cctaaaaGCGATTCAAATAGCGCCAGCTTAATCCGTTAAACCGCAAAGGTCACCCCTTTTCCGCCTTTGTCCGCAACAACACAACGCAGTGCACCTTCTTTCTCGGTTTAAAAACTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                      <  1:162859/146‑1 (MQ=255)
                               ccTACAAGCTATTCAAATAGCGCCAGCTTAATCGTTACAACAGCGAGGGCCACCCCCTTTTCGCCTTTGCCCGCAACAACATACCGCAGTCCCCTTTCTGTCTCGGTTTAAAAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGt                                                                                                                       <  1:207619/150‑1 (MQ=255)
                                 ttCAAGCGTTTCAGAGACCCCCATCTTAATCGTTACAACAGCGAAGGCCCCCCCTTTTCCGCGGTTGTCCGCGAAAACATAACGCAGTGCACCTTCTGTCTGGGTGAAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGt                                                                                                                       <  1:252581/148‑1 (MQ=255)
                                 ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                                                                                                    >  2:1080161/1‑151 (MQ=255)
                                 ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                                                                                                    >  1:405701/1‑151 (MQ=255)
                                    aaGCAATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGCCACCCCTTTTCCGCCGTTGTCCGCGAAAACATACCGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                                      <  2:964323/146‑1 (MQ=255)
                                      gCGATAAAGATGCCGCCATTTAAATCGTTACAACAGCGAAGGCCGCCCCTTTTTCGCCTTTCTCCGCGAAAACATAACGCGGTCCCCCTTCTGTCTGGGTTTAATACCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTg                                                                                                               <  1:424922/151‑1 (MQ=255)
                                         aTTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                            >  1:1104862/1‑151 (MQ=255)
                                          ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCtt                                                                                                           >  2:368134/1‑151 (MQ=255)
                                          tacAGATAGCCCCAGTTAAACCGTTTAAACAGCAAAGGCCAGCCCCTTTTCGCTTTGACCCGCAAAAACATAACGCATTCCACCTTCTGTCCGGGTGTAATACGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                            <  1:1057796/148‑1 (MQ=255)
                                                tagccccacttaaaTCGGTCCAACAGCGAGGGCCACCCCCTTTCCGCCGTTGCCCGCGACAACATACCGCAGTGCCCCTTCTGTCTGGGTTTAATACCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTgtgt                                                                                                      <  2:518096/138‑1 (MQ=255)
                                                   cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                                  >  2:237964/1‑151 (MQ=255)
                                                   ccccAGTTTAACCGGTTCAACAGCGAAGGTCACCCCCTTTTCGCCGTTGCCCGCGAAAACATAACGCAGTGCCCCTTCTGCCTGGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtt                                                                                                     <  2:1170634/146‑1 (MQ=255)
                                                     ccAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                                >  2:893555/1‑151 (MQ=255)
                                                          ttAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                            >  1:304827/1‑150 (MQ=255)
                                                              tCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                                                                                       >  1:780836/1‑151 (MQ=255)
                                                                gtttaaaCAGCGAGGGCCACCCCTTTTCCGCGGTGGTCCGCGACAACAAACCGCAGTGCCCCTTCTGTCTGGGTGTAATACCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg                                                                                     <  2:1109633/146‑1 (MQ=255)
                                                                gttccAACAGCGAAGGCCCCCCCCTTTTCGTCGTTTCCCGCAACAACATACCGCGTTGCCCCTTCTGTCTGGTTTTAAATGTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                                                                                         <  1:468576/143‑1 (MQ=255)
                                                                 ggtcaaacggcgaGGTCCACCCCTTTTCCGCGTTTGCCCGCAACAAAAAACCGCGTTCCCCCTCCTGTCTGGTTGAAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg                                                                                     <  1:1092807/139‑1 (MQ=255)
                                                                            agggccacccctTTTTCGCCGTTGTCCGCGAAAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATAGGTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGtgtg                                                                             <  1:1029053/139‑1 (MQ=255)
                                                                              ggccAGCCCTTTTTCGCGGTTGCCCGCGACAAGATAGCGCATTTCCCCTTCTTTCTCGGTGTAATATGGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                                       <  1:315643/148‑1 (MQ=255)
                                                                                 cacccccTTTCCGCCTTTGCCCGCGAAAAAATAACGCAGTCCCCCTTCGGTCTGGGTGTAAAACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAcc                                                                    <  2:783102/148‑1 (MQ=255)
                                                                                 cAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAcc                                                                    >  1:540276/1‑151 (MQ=255)
                                                                                 cAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAcc                                                                    >  1:837528/1‑151 (MQ=255)
                                                                                    ccccTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                                   >  1:362882/1‑149 (MQ=255)
                                                                                    ccccTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCt                                                                 >  2:327530/1‑151 (MQ=255)
                                                                                    ccccTTTGCGCTGTTGCCCGGGAAAACATAGCGGAGTGCTCCTTGTCTCCGGGTGTATTAATTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                                   <  2:435959/149‑1 (MQ=255)
                                                                                      ccTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGc                                                               >  2:279937/1‑151 (MQ=255)
                                                                                         tttCGCCTTTGTCCGCGACAACATAACGCATTGCCCTTTCTGTCTGGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCa                                                             <  2:1025950/150‑1 (MQ=255)
                                                                                           tccGCCGTTGTCCGCAACAACATACCGCTGTGCCCTTTCTGTCTGGGTGTAATACCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCAt                                                            <  2:206091/148‑1 (MQ=255)
                                                                                                    gTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTg                                                   >  2:171872/1‑150 (MQ=255)
                                                                                                      ccGCAAAAACATACCGCAGTGCCCCTTCTGTCTCGGTTTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                                  <  2:129438/149‑1 (MQ=255)
                                                                                                           acaacaTAACGCAGTGCACCTTCTTTCTCGGTGTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttc                                            <  2:822253/150‑1 (MQ=255)
                                                                                                           acaacaTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCACTGCCATGTATTCCTGAgctttcg                                           >  1:657084/1‑151 (MQ=255)
                                                                                                             aacaTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgct                                         >  1:380980/1‑151 (MQ=255)
                                                                                                             aaaaaaccgcgGTGCCCCTTTGGTCTGGGTGTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgct                                         <  1:893350/144‑1 (MQ=255)
                                                                                                                aTACCGCGGTGCCCCTTCTGTCTCGGTGTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctttc                                      <  2:358648/151‑1 (MQ=255)
                                                                                                                aTAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctttc                                      <  2:493565/151‑1 (MQ=255)
                                                                                                                 taccgcgGTGCCCTTTTTGTCTCGGTGTAATATTTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                        <  1:926186/145‑1 (MQ=255)
                                                                                                                  accGCAGTGCCCCTTCTGTCTGGGTGAAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctttc                                      <  2:310239/147‑1 (MQ=255)
                                                                                                                      cAGTGCCCCTTCTGTCGCGGTGTAATACCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                                 <  2:1005805/150‑1 (MQ=255)
                                                                                                                         tGCACCTTCTGTCTCGGTTTAATAACTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                              <  2:1096074/150‑1 (MQ=255)
                                                                                                                           ccccTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttttt                            <  2:8211/148‑1 (MQ=255)
                                                                                                                             ccTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                         >  2:826359/1‑151 (MQ=255)
                                                                                                                                   gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc                   >  2:424659/1‑151 (MQ=255)
                                                                                                                                    tctcGGTGAAATAGTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc                   <  1:660277/150‑1 (MQ=255)
                                                                                                                                        gttttaatacttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                       <  1:747047/130‑1 (MQ=255)
                                                                                                                                        ggTGTAATACGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTg                      <  1:939608/143‑1 (MQ=255)
                                                                                                                                           gtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg           >  1:660928/1‑151 (MQ=255)
                                                                                                                                             aataaGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg          <  1:1050997/150‑1 (MQ=255)
                                                                                                                                             aataaCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc                   <  2:480288/141‑1 (MQ=255)
                                                                                                                                              atacttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGt        <  1:574537/145‑1 (MQ=255)
                                                                                                                                                    ttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGt        <  2:479177/145‑1 (MQ=255)
                                                                                                                                                    ttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTa       <  2:389326/146‑1 (MQ=255)
                                                                                                                                                      tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCTCCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGCCAg  <  2:949475/149‑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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: