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

CGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  NC_000913/2001771‑2002064
                                                                                                                                                       |                                                                                                                                                   
cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑ccc                                                                                                                                                 >  1:514286/1‑151 (MQ=255)
          tGAATTCACGATAAACAGCCCTGCG‑TAAATGGGGTAACGGCCTTGATTTTCGGTTTCTGAAGGGTTTTTATTGGGCCGTTTGGGGGGGGAATGGGGGGGTTTTTT‑GGGGG‑TT‑CAGGGAAAAT‑TCAGGGAGGAAAAAA‑ACCCGCCGgggggg                                                                                                                                      >  1:277429/1‑149 (MQ=255)
            aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGGGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCgg                                                                                                                                     >  2:243482/1‑151 (MQ=255)
             atccaagaaaaaCGGCCTTGGTTTATATGAGTTATGGG‑CATGTTTTTCCTTTTCTGGGGGTTTTTTAAGGGGACATTTAGGGGGTAAAAGGGGGGTTTATTT‑GGGGG‑TTACAGGTAAATC‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCggg                                                                                                                                    <  2:520980/142‑1 (MQ=255)
               tCACGATAAACAGCCCTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCgggg                                                                                                                                   >  1:452382/1‑150 (MQ=255)
                caggaaaaccAGCCCTGCTTTATATGAGTTATCGG‑CATAATTTTCCTTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTATTG‑GGGGG‑TAAAAGAAAAATA‑CAAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGaa                                                                                                                                 <  1:1102757/142‑1 (MQ=255)
                cacattAAACACCCCTTCTTTAATGGATTTAGCGG‑GATGTTTTTCCTTTTCAGCGGTGTTTTAAATCGAAGGATTGGGGGGGAAAATGAGGGTTTTTTT‑GGGGT‑TTAAAGGAAAATT‑CCAGGCAAAAAAAAA‑CCCCGCCGGTGGCGGGGaa                                                                                                                                 <  2:938099/146‑1 (MQ=255)
                      aaaCAGCCCTGCGTTATATGATTTATCGG‑CATGATTATCCGTTTCTGCAGGTTTTTTAATGGGACGATTAGGGGGTAAAAGGGGGGTTTATTG‑GGGGG‑TAAAAGTTAATTC‑CCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACg                                                                                                                            <  1:109680/150‑1 (MQ=255)
                       aaCAGCCCTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                          >  1:524350/1‑151 (MQ=255)
                         cAGCCCTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCa                                                                                                                              >  1:916228/1‑145 (MQ=255)
                          aGCCCTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCt                                                                                                                       >  2:392351/1‑151 (MQ=255)
                             ccTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                     >  2:1021574/1‑150 (MQ=255)
                                  gTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                               >  2:784888/1‑151 (MQ=255)
                                   ttATATGATTTAGCGG‑CATGATTATCCTTTTCGGCAGGGTTTTAAAGCGGACGATTAGGGGGTAAAATGGGGGTTTTTTG‑GGGGG‑TTACGGGTAAATC‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                <  2:927216/149‑1 (MQ=255)
                                      tatGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                           >  2:990174/1‑151 (MQ=255)
                                      tatGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAATgcgc                                                                                                           >  1:1064992/1‑151 (MQ=255)
                                        ggAGTTATGGG‑CATGGTTAAGCGTTTGTGTAGGGTTTTTAGTGTGGCGATTGGTGGGTGACATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                         >  1:1104682/2‑151 (MQ=255)
                                            ttACCGG‑AATGATTATCCGTTTCTGCGGGTTTTTTAACCGGACGATTAGTGGGTGAAATGGGGGTTTTTTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                             <  2:190966/143‑1 (MQ=255)
                                             tATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                    >  1:163276/1‑151 (MQ=255)
                                             tATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                    >  1:949269/1‑151 (MQ=255)
                                              aTCGG‑CATGTTTCTCCTTTTCTGCGGGTTTTTAAAGGGGACGATTAGGGGGTAAAAGGGGGTGTTTTTG‑GGGGT‑TTACAGGAAAATC‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                    <  1:938439/150‑1 (MQ=255)
                                                  g‑CATGTTTATCGTTTTCTGCGGGTTTTTAAATCGGACGTTTGGGGGGTAAAAGGGGGGTTTTTTG‑GGGGG‑TTAAAGGAAAATT‑CCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                   <  1:590752/147‑1 (MQ=255)
                                                    cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                              >  2:919645/1‑151 (MQ=255)
                                                       taTTATCCGTTTCTGCAGGGTTTTTATTCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGAAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgca                                                                                             >  2:1254811/2‑147 (MQ=255)
                                                         tttcccttttCTGCAGGGTTTTAAAGCGGAAGTTGAGTGGGTGAAAGGAGGGGTT‑TTT‑GGGGG‑TTAAAAGTTAATTACCAGCAAAAAAAACA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                          <  1:1236636/143‑1 (MQ=255)
                                                         ttATCCTTTTCTGCAGGTTTTTTAACCGGACGTTTAGTGGGTAAAAGGGGGGTTTTTTT‑GGGGG‑TTACAGGTAAATT‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                           <  1:803103/149‑1 (MQ=255)
                                                             ccttttCTGCAGGGTTTTTAACCGGACGATTAGTGGGTGAAATGAGGGTTTTTTT‑GGGGG‑TTACAGGTAAATT‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                      <  2:30253/147‑1 (MQ=255)
                                                             ccGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                     >  2:501069/1‑151 (MQ=255)
                                                              cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                     >  2:330075/1‑150 (MQ=255)
                                                               gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                                    >  2:478702/1‑150 (MQ=255)
                                                                tttCGGCGGGGTTTTAAAGCGGACGATTAGGGGGTGAAATGGGGGGTTATTG‑GGGGG‑TTACAGGTAAATC‑CCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                      <  2:441856/147‑1 (MQ=255)
                                                                 ttgtGCGGGGTTTTTAGGCGGAGGATTGGTGGGTGAAGTGGGGGTTTTTTT‑GGGGG‑TTACAGTTAATTC‑CCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                                  <  2:358957/147‑1 (MQ=255)
                                                                        ggttttttAATCGGCCGATTAGGGGGTAAAATGAGGGTTTATTT‑GGGGG‑TTACAGGAAAATC‑CCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          <  2:339609/148‑1 (MQ=255)
                                                                        ggtttttaAATGGGACGATT‑GGGGGTTAAAAGGGGGGTTTTTT‑GGGGGTTTACAGGAAAATC‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               <  1:1135818/143‑1 (MQ=255)
                                                                        gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          >  2:1169603/1‑151 (MQ=255)
                                                                        gggTTTTTAACCGGACGTTGAGGGGATAAAAAGAGGGTTTTTGG‑GGGGG‑TAAGAGGTAAATC‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          <  1:1166988/151‑1 (MQ=255)
                                                                         ggTTTTAAACCGGAGGATGAGTGGGGAAAAGGGGGGTTTTTTT‑GGGGG‑TTACAGGAAAATT‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          <  1:149647/150‑1 (MQ=255)
                                                                            ttttAATCGGCCGATTAGGGGGAAAAATGAGGGGTTATTT‑GGGGT‑TTACAGGTAATTC‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                       <  2:315394/150‑1 (MQ=38)
                                                                               tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                   >  1:772294/1‑151 (MQ=255)
                                                                                aaccGGACGTTTGGGGGGTAAAAGGGGGGTTTATTG‑GGGGG‑TTACAGGAAAATC‑CCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                  <  1:456522/148‑1 (MQ=37)
                                                                                    ggACGTTTAGGGGGTAAAATGGGGGGTTTTTT‑GGGGG‑TTACAGGAAATTC‑CCAGCAAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                              <  1:57637/151‑1 (MQ=37)
                                                                                    ggACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                               >  1:603156/1‑150 (MQ=255)
                                                                                       cGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                            >  2:882191/1‑150 (MQ=255)
                                                                                        gATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                          >  2:323112/1‑151 (MQ=255)
                                                                                           tAGGGGGTAAATGGAGGGGTTTTTT‑GGGGG‑TAAAGGGTAAATT‑CAAGAAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                            <  1:679491/146‑1 (MQ=25)
                                                                                             gTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                      >  1:648719/1‑150 (MQ=255)
                                                                                              tGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                     >  2:668319/1‑150 (MQ=255)
                                                                                               gggTAAAATGAGGGTTTATTG‑GGGGG‑TTACAGGTAAATT‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                   <  2:452431/151‑1 (MQ=38)
                                                                                                     aaaggggggTTTTTTGGGGGG‑TAACAGG‑AAATT‑CCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                <  1:151991/144‑1 (MQ=25)
                                                                                                     aaTGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                              >  1:296214/1‑150 (MQ=255)
                                                                                                        gAGGGGTTTTTG‑GGGTT‑TTAAAGGAAAATC‑CCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                            <  1:1225738/149‑1 (MQ=37)
                                                                                                          gggtTTTTTG‑GGGGG‑TTAAAGGAAAATC‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                         <  1:154788/150‑1 (MQ=25)
                                                                                                          gggtTTATGG‑GGGGG‑TTACAGGTAAATC‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  2:514333/151‑1 (MQ=37)
                                                                                                           ggtttATGT‑GGGGT‑TTA‑AAGTAAATT‑TCTGGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                         <  1:243457/146‑1 (MQ=25)
                                                                                                           gggTTTTTT‑GGGGGTTTAAATGAAAATA‑CCAGGAAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  1:1179603/151‑1 (MQ=25)
                                                                                                           gggTTATTG‑GGGGG‑TTACAGGAAAATC‑CCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  2:1076590/150‑1 (MQ=255)
                                                                                                            gtttATGT‑GGGGT‑TAAGAGTAAAATG‑CCAGCAAAAAAAACA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                      <  1:159606/149‑1 (MQ=21)
                                                                                                                      gggt‑TTACAGGTAAATC‑CCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                 <  1:220797/147‑1 (MQ=37)
                                                                                                                       ggg‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                >  1:148665/1‑147 (MQ=255)
                                                                                                                       ggg‑TTACAGGTAAATC‑CCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                <  2:148677/147‑1 (MQ=255)
                                                                                                                           ttACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                           >  1:387697/1‑149 (MQ=255)
                                                                                                                                   aattt‑CCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                  <  2:277457/147‑1 (MQ=255)
                                                                                                                                   aaaTC‑CCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                  <  2:215599/150‑1 (MQ=255)
                                                                                                                                      tt‑CCAGCCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag               <  2:1023799/150‑1 (MQ=255)
                                                                                                                                         ccAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc            >  1:905187/1‑151 (MQ=255)
                                                                                                                                           aGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc            >  2:284229/1‑149 (MQ=255)
                                                                                                                                            ggCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag          <  2:163291/150‑1 (MQ=255)
                                                                                                                                             gCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag          <  1:219587/149‑1 (MQ=255)
                                                                                                                                                gAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt     >  2:1137806/1‑151 (MQ=255)
                                                                                                                                                   aaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg  <  1:990077/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGG‑CATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  NC_000913/2001771‑2002064

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