Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC42.6% 231.0 / 26.7 115intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (61/5);  new base C (0/49);  total (61/54)
Fisher's exact test for biased strand distribution p-value = 3.58e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001769‑2002059
                                                                                                                                                       |                                                                                                                                             
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAAc                                                                                                                                              >  1:365040/1‑151 (MQ=255)
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAAc                                                                                                                                              >  2:877990/1‑151 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAAcc                                                                                                                                             <  2:542314/151‑1 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAAcccc                                                                                                                                           <  1:724282/151‑1 (MQ=255)
      ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCccgcc                                                                                                                                        >  2:654933/1‑151 (MQ=255)
      ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCccgcc                                                                                                                                        >  2:483294/1‑151 (MQ=255)
        agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCgg                                                                                                                                      >  2:1076281/1‑151 (MQ=255)
          agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTg                                                                                                                                    >  1:749552/1‑151 (MQ=255)
           gTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTgg                                                                                                                                   >  1:1086074/1‑151 (MQ=255)
            tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGc                                                                                                                                  >  1:1778372/1‑151 (MQ=255)
              aaTTCACAAAAAAAAGCCCGGGTTTTTATGAGTTATCGGCATGATTTTCCGTTTCCGCAGGGTTTTTAATGGGACAGTTAGTGGGTAAAATGAGGGGTTATTTGGGGTTT‑ACAGGAAAATTCC‑GGGCAAAAAAAAACCCCGCCGGTGGCgg                                                                                                                                <  1:1155635/151‑1 (MQ=255)
               aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                               >  2:1145280/1‑151 (MQ=255)
                      aTAAACAGCCCTGCTTTATATGAGTTATCGGATGGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGGAAAAGGGGGGTTTTTTTGGGGGTT‑ACGGGTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                          <  2:386130/149‑1 (MQ=255)
                         aaCAGCCCTGCGTTATATGATTTATCGGCATGATTATCCGTTTCTGCGGGTTTTTTAATGGGACGATTAGTGGGTAAAATGGGGGTTTATTGGGGGGTT‑ACAGGTAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                     <  1:1110768/151‑1 (MQ=255)
                            aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCt                                                                                                                  >  2:1593870/1‑151 (MQ=255)
                              ccttGCTTTATATGAGTTATCGGCATGTTTATCCTTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGTTTATTGGGGGGTT‑ACAGTTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                                    <  1:579247/144‑1 (MQ=255)
                                 tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                              >  2:658235/1‑150 (MQ=255)
                                 tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                              >  2:1180010/1‑150 (MQ=255)
                                   ctttATATGGTTTAGCGGGTTGATGATCCTTTTCGGGGGGTTTTTTATCCGCACGTTTAGGGGGTGAAATGAGGGTTGATTTGGGGTTT‑ACGGGTAATTCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                              <  1:1001592/146‑1 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCTTTTCGGCAGGTTTTTTAAGCGGACGATTAGGGGGTGAAATGGGGGTTTATTGGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                              <  1:1288801/146‑1 (MQ=255)
                                      tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                        >  2:697164/1‑151 (MQ=255)
                                        tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                      >  1:454560/1‑151 (MQ=255)
                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                     >  1:42741/1‑151 (MQ=255)
                                          tGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                    >  2:1345375/1‑151 (MQ=255)
                                           gatttAGCGGCATTATTATCGGTTTCGGCAGGTTTTTAAATCGGACGATTGGGGGGGAAAAGGGGGGGTTTTTTGGGGGTT‑AAAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                      <  2:1548264/145‑1 (MQ=255)
                                            aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   >  1:406107/1‑150 (MQ=255)
                                            aGTTATCGGAATGATTATCTTTTTCTGAGGTTTTTTTAATCGGACGATTAGTGGGAGAAATGGGGGTTTTTTTGGGGTTT‑AAAGGAAAATTCA‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                         <  2:262135/144‑1 (MQ=255)
                                              ttATCGGCATGATTATCCGTTTCTGCGGGTTTTTTAATCGAACAATTGGGGGGTAAAAGGGGGGTTTTTTTGGGGGTT‑AAAGGAAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                    <  2:189913/147‑1 (MQ=255)
                                              ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                 >  2:919271/1‑150 (MQ=255)
                                              ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                                >  1:1519162/1‑151 (MQ=255)
                                              ttATCGGAAGGATTATCCGTTTCTGCGGGTTTTTAAACCGGAGGATTAGGGGTTAAAAGGAGGGTTTTTTTGGGGGTT‑ACAGGTAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                    <  2:1362866/147‑1 (MQ=255)
                                                   ggCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                           >  1:223214/1‑151 (MQ=255)
                                                   ggCATGATTACCCGTTTCTGCGGGTTTTTTAATGGGACGATTGGGGGGTGAAATGGGGGTTTTTTTGGGGGTT‑AAAGGAAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                               <  2:119287/147‑1 (MQ=255)
                                                    gCATGTTTTTCCTTTTCTGCGGGTTTTTTAATGGGACGATTAGGGGGTAAAAGGGGGGTTTTTTGGGGGGTT‑ACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                          <  1:1076344/151‑1 (MQ=255)
                                                    gCATGATTATCCTTTTCTGCAGTTTTTTAAAGGGGACGATTAGTGGGTAAAATGA‑GGGTTTTTTGGGGGTTAACAGTTAAATTCC‑GGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                          <  1:1227308/151‑1 (MQ=255)
                                                     cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         >  1:940214/1‑151 (MQ=255)
                                                      aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         >  1:80442/1‑150 (MQ=255)
                                                      aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         >  1:484406/1‑150 (MQ=255)
                                                        gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                       >  1:990662/1‑150 (MQ=255)
                                                          ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                    >  1:703000/1‑151 (MQ=255)
                                                                gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                               >  2:938144/1‑150 (MQ=255)
                                                                   tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                           >  1:1240634/1‑151 (MQ=255)
                                                                      gCAGGGTTTTTAATCGGAGGATTGGGGGGGGAAATGAGGGTTTATTTGGGGGTT‑ACAGGTAATTCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                            <  1:1584543/147‑1 (MQ=255)
                                                                             ttttAATGGGACGATTGGTGGGTAAAATGGGGGTTTATTTGGGGGTT‑ACGGGTAAATTCC‑AGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                  <  1:277117/150‑1 (MQ=39)
                                                                             ttttAATCGGACGATTGGGGGGTAAAATGGGGGTTTTTTTGGGGGTT‑ACAGGTAATTCCC‑GGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                  <  2:895683/150‑1 (MQ=38)
                                                                             ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                 >  1:543204/1‑151 (MQ=255)
                                                                              tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                >  1:1294789/1‑151 (MQ=255)
                                                                               ttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                               >  1:969587/1‑151 (MQ=255)
                                                                                tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                              >  2:684855/1‑151 (MQ=255)
                                                                                 aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                             <  2:1133541/151‑1 (MQ=255)
                                                                                 aaTCGGACGATTAGGGGGTGAAAGGAGGGTTTTTTTGGGGGTT‑ACAGGTAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                              <  1:1646965/150‑1 (MQ=38)
                                                                                    cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                           >  1:725969/1‑150 (MQ=255)
                                                                                       aggATTGGTGGGTGAAATGAGGGGTTTTTTGGGGTTA‑ACAGGAAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        <  2:195421/148‑1 (MQ=37)
                                                                                       aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        >  2:1753058/1‑150 (MQ=255)
                                                                                       aCGATTAGTGGGTAAAATGGGGGTTTATTGGGGGGTT‑ACAGGTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        <  1:193795/150‑1 (MQ=38)
                                                                                        cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                      >  2:496561/1‑151 (MQ=255)
                                                                                          aTTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                    >  1:1124140/1‑151 (MQ=255)
                                                                                           ttAGTGGGTAAAATGAGGGTTTTTTTGGGGTTT‑ACAGGTAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                     <  2:1643169/149‑1 (MQ=37)
                                                                                           ttAGTGGGTAAAAGGAGGGTTTATTTGGGGTTT‑ACAGGAAAATCCC‑AGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                      <  1:1703640/148‑1 (MQ=38)
                                                                                             aGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  >  1:1691821/1‑150 (MQ=255)
                                                                                              gTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                 >  2:1505192/1‑150 (MQ=255)
                                                                                                gggTAAAAGGGGGGTTTTTTGGGGGGTT‑AAAGGAAAATCCC‑GGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                                <  1:367359/149‑1 (MQ=25)
                                                                                                gggTAAAAGGAGGGTTTTTTTGGGGGTT‑AAAGAAAAATCCC‑AGACAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                 <  2:940175/148‑1 (MQ=25)
                                                                                                     aaaTGGGGGTTTTTTTGGGGGTA‑ACAGATAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                         <  2:1748765/151‑1 (MQ=25)
                                                                                                     aaaGGAGGGGTTTTTTGGGGGTT‑ACATGAAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa                                          <  2:734607/150‑1 (MQ=37)
                                                                                                         ggggggTTTTTGGGGGTTT‑ACAGGTAATTCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                     <  1:1588388/149‑1 (MQ=37)
                                                                                                         gAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                     >  2:1485888/1‑151 (MQ=255)
                                                                                                         gAGGGGTTATTGGGGGTTT‑ACAGGAAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                     <  1:1185653/151‑1 (MQ=37)
                                                                                                           gggtTTTTTTGGGGTTT‑ACGGGAAAATGCA‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                      <  2:1565730/148‑1 (MQ=25)
                                                                                                           gggtTTTTTTGGGGGTA‑ACAGGTAAATTCA‑GGCCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATGGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                     <  2:107407/149‑1 (MQ=25)
                                                                                                            ggtttttttGGGGGTT‑ACAGATAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                  <  2:341650/145‑1 (MQ=37)
                                                                                                            gggTTTTTTGGGGGTT‑AAAGATAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                  <  2:335549/151‑1 (MQ=37)
                                                                                                                tATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                              >  1:1743139/1‑151 (MQ=255)
                                                                                                                 aTTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             >  2:1234851/1‑151 (MQ=255)
                                                                                                                  tttGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             >  1:1059477/1‑150 (MQ=255)
                                                                                                                  tttGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             >  1:581230/1‑150 (MQ=255)
                                                                                                                  tttGGGGGTT‑ACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                              <  1:758473/149‑1 (MQ=255)
                                                                                                                  ttgggggtTT‑ACAGGTAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             <  1:1669478/147‑1 (MQ=37)
                                                                                                                  ttggggggTT‑ACAGGTAAATTCC‑GGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             <  1:799294/147‑1 (MQ=255)
                                                                                                                     ggggtTT‑ACATATAAATGCG‑AGGGAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTACGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                           <  1:1426845/149‑1 (MQ=37)
                                                                                                                     gggggTT‑ACAGTTAAATCCC‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                          <  1:1354704/150‑1 (MQ=255)
                                                                                                                      ggggTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                        >  2:1723796/1‑151 (MQ=255)
                                                                                                                       gggTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACCCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCCACtt                        >  1:142848/1‑150 (MQ=255)
                                                                                                                        ggTT‑ACAGGTAATTCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                      <  2:434546/151‑1 (MQ=255)
                                                                                                                        ggTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                      >  1:1209577/1‑151 (MQ=255)
                                                                                                                          tt‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                      >  2:1708527/1‑149 (MQ=255)
                                                                                                                          tt‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                      >  1:1737969/1‑149 (MQ=255)
                                                                                                                           taaaaGGTAATTTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                     <  1:1307448/146‑1 (MQ=37)
                                                                                                                           t‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                   >  1:1484859/1‑151 (MQ=255)
                                                                                                                             acgggaaaaTCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                  <  1:496582/145‑1 (MQ=37)
                                                                                                                                   aaaTTCC‑AGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa            <  1:790347/151‑1 (MQ=255)
                                                                                                                                   aaaTTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa              <  1:173658/149‑1 (MQ=255)
                                                                                                                                   aaaTCCC‑AGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg               <  2:1007691/148‑1 (MQ=255)
                                                                                                                                     aTTCCAAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa           <  2:456013/151‑1 (MQ=255)
                                                                                                                                      tccc‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag          <  1:730649/148‑1 (MQ=255)
                                                                                                                                           aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc       >  2:641368/1‑149 (MQ=255)
                                                                                                                                            ggaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag     <  1:202123/145‑1 (MQ=255)
                                                                                                                                            ggaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca      <  2:333501/144‑1 (MQ=255)
                                                                                                                                               aaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  <  2:1026601/148‑1 (MQ=255)
                                                                                                                                               aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  1:32336/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTT‑ACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001769‑2002059

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