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

CGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  NC_000913/2001771‑2002064
                                                                                                                                                       |                                                                                                                                                   
cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAccc                                                                                                                                                 >  1:514286/1‑151 (MQ=255)
          tGAATTCACGATAAACAGCCCTGCGTAAATGGGGTAACGGCCTTGATTTTCGGTTTCTGAAGGGTTTTTATTGGGCCGTTTGGGGGGGGAATGGGGGGGTTTTTTGGGGGTTCAGGGAAAATTCAGGGAGGAAAAAAACCCGCCGgggggg                                                                                                                                      >  1:277429/1‑149 (MQ=255)
            aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGGGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgg                                                                                                                                     >  2:243482/1‑151 (MQ=255)
             atccaagaaaaaCGGCCTTGGTTTATATGAGTTATGGGCATGTTTTTCCTTTTCTGGGGGTTTTTTAAGGGGACATTTAGGGGGTAAAAGGGGGGTTTATTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCggg                                                                                                                                    <  2:520980/142‑1 (MQ=255)
               tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                   >  1:452382/1‑150 (MQ=255)
                caggaaaaccAGCCCTGCTTTATATGAGTTATCGGCATAATTTTCCTTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTATTGGGGGGTAAAAGAAAAATACAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                 <  1:1102757/142‑1 (MQ=255)
                cacattAAACACCCCTTCTTTAATGGATTTAGCGGGATGTTTTTCCTTTTCAGCGGTGTTTTAAATCGAAGGATTGGGGGGGAAAATGAGGGTTTTTTTGGGGTTTAAAGGAAAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGaa                                                                                                                                 <  2:938099/146‑1 (MQ=255)
                      aaaCAGCCCTGCGTTATATGATTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATGGGACGATTAGGGGGTAAAAGGGGGGTTTATTGGGGGGTAAAAGTTAATTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACg                                                                                                                            <  1:109680/150‑1 (MQ=255)
                       aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                          >  1:524350/1‑151 (MQ=255)
                         cAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCa                                                                                                                              >  1:916228/1‑145 (MQ=255)
                          aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCt                                                                                                                       >  2:392351/1‑151 (MQ=255)
                             ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                     >  2:1021574/1‑150 (MQ=255)
                                  gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                               >  2:784888/1‑151 (MQ=255)
                                   ttATATGATTTAGCGGCATGATTATCCTTTTCGGCAGGGTTTTAAAGCGGACGATTAGGGGGTAAAATGGGGGTTTTTTGGGGGGTTACGGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                <  2:927216/149‑1 (MQ=255)
                                      tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                           >  2:990174/1‑151 (MQ=255)
                                      tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAATgcgc                                                                                                           >  1:1064992/1‑151 (MQ=255)
                                        ggAGTTATGGGCATGGTTAAGCGTTTGTGTAGGGTTTTTAGTGTGGCGATTGGTGGGTGACATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                         >  1:1104682/2‑151 (MQ=255)
                                            ttACCGGAATGATTATCCGTTTCTGCGGGTTTTTTAACCGGACGATTAGTGGGTGAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                             <  2:190966/143‑1 (MQ=255)
                                             tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                    >  1:163276/1‑151 (MQ=255)
                                             tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                    >  1:949269/1‑151 (MQ=255)
                                              aTCGGCATGTTTCTCCTTTTCTGCGGGTTTTTAAAGGGGACGATTAGGGGGTAAAAGGGGGTGTTTTTGGGGGTTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                    <  1:938439/150‑1 (MQ=255)
                                                  gCATGTTTATCGTTTTCTGCGGGTTTTTAAATCGGACGTTTGGGGGGTAAAAGGGGGGTTTTTTGGGGGGTTAAAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                   <  1:590752/147‑1 (MQ=255)
                                                    cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                              >  2:919645/1‑151 (MQ=255)
                                                       taTTATCCGTTTCTGCAGGGTTTTTATTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGAAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgca                                                                                             >  2:1254811/2‑147 (MQ=255)
                                                         tttcccttttCTGCAGGGTTTTAAAGCGGAAGTTGAGTGGGTGAAAGGAGGGGTTTTTGGGGGTTAAAAGTTAATTACCAGCAAAAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                          <  1:1236636/143‑1 (MQ=255)
                                                         ttATCCTTTTCTGCAGGTTTTTTAACCGGACGTTTAGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                           <  1:803103/149‑1 (MQ=255)
                                                             ccttttCTGCAGGGTTTTTAACCGGACGATTAGTGGGTGAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                      <  2:30253/147‑1 (MQ=255)
                                                             ccGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                     >  2:501069/1‑151 (MQ=255)
                                                              cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                     >  2:330075/1‑150 (MQ=255)
                                                               gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                                    >  2:478702/1‑150 (MQ=255)
                                                                tttCGGCGGGGTTTTAAAGCGGACGATTAGGGGGTGAAATGGGGGGTTATTGGGGGGTTACAGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                      <  2:441856/147‑1 (MQ=255)
                                                                 ttgtGCGGGGTTTTTAGGCGGAGGATTGGTGGGTGAAGTGGGGGTTTTTTTGGGGGTTACAGTTAATTCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa                                                                                  <  2:358957/147‑1 (MQ=255)
                                                                        ggttttttAATCGGCCGATTAGGGGGTAAAATGAGGGTTTATTTGGGGGTTACAGGAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          <  2:339609/148‑1 (MQ=255)
                                                                        ggtttttaAATGGGACGATTGGGGGTTAAAAGGGGGGTTTTTTGGGGGTTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               <  1:1135818/143‑1 (MQ=255)
                                                                        gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          >  2:1169603/1‑151 (MQ=255)
                                                                        gggTTTTTAACCGGACGTTGAGGGGATAAAAAGAGGGTTTTTGGGGGGGTAAGAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          <  1:1166988/151‑1 (MQ=255)
                                                                         ggTTTTAAACCGGAGGATGAGTGGGGAAAAGGGGGGTTTTTTTGGGGGTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                          <  1:149647/150‑1 (MQ=255)
                                                                            ttttAATCGGCCGATTAGGGGGAAAAATGAGGGGTTATTTGGGGTTTACAGGTAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                       <  2:315394/150‑1 (MQ=38)
                                                                               tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                   >  1:772294/1‑151 (MQ=255)
                                                                                aaccGGACGTTTGGGGGGTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                  <  1:456522/148‑1 (MQ=37)
                                                                                    ggACGTTTAGGGGGTAAAATGGGGGGTTTTTTGGGGGTTACAGGAAATTCCCAGCAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                              <  1:57637/151‑1 (MQ=37)
                                                                                    ggACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                               >  1:603156/1‑150 (MQ=255)
                                                                                       cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                            >  2:882191/1‑150 (MQ=255)
                                                                                        gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                          >  2:323112/1‑151 (MQ=255)
                                                                                           tAGGGGGTAAATGGAGGGGTTTTTTGGGGGTAAAGGGTAAATTCAAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                            <  1:679491/146‑1 (MQ=25)
                                                                                             gTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                      >  1:648719/1‑150 (MQ=255)
                                                                                              tGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                     >  2:668319/1‑150 (MQ=255)
                                                                                               gggTAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                   <  2:452431/151‑1 (MQ=38)
                                                                                                     aaaggggggTTTTTTGGGGGGTAACAGGAAATTCCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                <  1:151991/144‑1 (MQ=25)
                                                                                                     aaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                              >  1:296214/1‑150 (MQ=255)
                                                                                                        gAGGGGTTTTTGGGGTTTTAAAGGAAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                            <  1:1225738/149‑1 (MQ=37)
                                                                                                          gggtTTTTTGGGGGGTTAAAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                         <  1:154788/150‑1 (MQ=25)
                                                                                                          gggtTTATGGGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  2:514333/151‑1 (MQ=37)
                                                                                                           ggtttATGTGGGGTTTAAAGTAAATTTCTGGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                         <  1:243457/146‑1 (MQ=25)
                                                                                                           gggTTTTTTGGGGGTTTAAATGAAAATACCAGGAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  1:1179603/151‑1 (MQ=25)
                                                                                                           gggTTATTGGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  2:1076590/150‑1 (MQ=255)
                                                                                                            gtttATGTGGGGTTAAGAGTAAAATGCCAGCAAAAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                      <  1:159606/149‑1 (MQ=21)
                                                                                                                      gggtTTACAGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                 <  1:220797/147‑1 (MQ=37)
                                                                                                                       gggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                >  1:148665/1‑147 (MQ=255)
                                                                                                                       gggTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                <  2:148677/147‑1 (MQ=255)
                                                                                                                           ttACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                           >  1:387697/1‑149 (MQ=255)
                                                                                                                                   aatttCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                  <  2:277457/147‑1 (MQ=255)
                                                                                                                                   aaaTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                  <  2:215599/150‑1 (MQ=255)
                                                                                                                                      ttCCAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag               <  2:1023799/150‑1 (MQ=255)
                                                                                                                                         ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc            >  1:905187/1‑151 (MQ=255)
                                                                                                                                           aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc            >  2:284229/1‑149 (MQ=255)
                                                                                                                                            ggCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag          <  2:163291/150‑1 (MQ=255)
                                                                                                                                             gCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag          <  1:219587/149‑1 (MQ=255)
                                                                                                                                                gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt     >  2:1137806/1‑151 (MQ=255)
                                                                                                                                                   aaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg  <  1:990077/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  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: