Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA45.1% 68.4 / 37.9 71intergenic (‑122/+199)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (33/6);  new base A (0/32);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.75e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAT  >  NC_000913/2001766‑2002061
                                                                                                                                                     |                                                                                                                                                       
ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaa                                                                                                                                                     >  2:22245/1‑150 (MQ=255)
 tCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaa                                                                                                                                                    <  1:784261/150‑1 (MQ=255)
  cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaa                                                                                                                                                   >  2:751832/1‑150 (MQ=255)
     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)
                                                                                                                                                     |                                                                                                                                                       
TTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAT  >  NC_000913/2001766‑2002061

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: