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

CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001768‑2002065
                                                                                                                                                        |                                                                                                                                                   
cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaaa                                                                                                                                                     <  2:1167315/151‑1 (MQ=255)
 aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                    >  1:225776/1‑151 (MQ=255)
  gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                   <  2:1481258/151‑1 (MQ=255)
   cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAccc                                                                                                                                                  >  1:548539/1‑151 (MQ=255)
     gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                                >  2:1186139/1‑151 (MQ=255)
       ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                               <  2:122720/150‑1 (MQ=255)
       ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                               <  1:796517/150‑1 (MQ=255)
       ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                               <  1:135036/150‑1 (MQ=255)
       ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                               <  2:1285553/150‑1 (MQ=255)
        taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCGGGCAGAAAAAAACCccgcc                                                                                                                                              <  2:600435/150‑1 (MQ=255)
         agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                            <  1:314830/151‑1 (MQ=255)
            gtaaaTCCACGATAAACAGCCCTGTGTTTTAAGGGTTACGGGCAGGATTATCCGTTTCTGCAGGGTTTTTATTCGGACGATGAGGGGGTGAAATGAGGGGTTATTTGGGGGTTAGAGGAAAATTCCAGAAAAAAAAAACCCccgccg                                                                                                                                             <  2:611623/144‑1 (MQ=255)
                aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                                     >  2:341218/1‑151 (MQ=255)
                   cACGAAAAAAAGCCCTGCTTTATAGGAGTTATCGGCATGATTTTCCGTTTCGGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTTTTTGGGGTTTACAGGAAAATCCCAGGCAGAAAAAAACCCCGCCGGTGGCgg                                                                                                                                      <  1:1496521/147‑1 (MQ=255)
                    aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                   >  1:336572/1‑149 (MQ=255)
                      gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                   >  2:1385560/1‑147 (MQ=255)
                       aTAAACAGCCCTGTGTTATATGAGTTATCGGCATGATTTTCCTTTTCTGCAGGTTTTTTAACCGGACGTTTGGTGGGTAAAATGGGGGTTTTTTTGGGGGTTAAAGTTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              <  2:1101829/151‑1 (MQ=255)
                       aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              >  1:978156/1‑151 (MQ=255)
                         aaaCAGCCCTGCTTTATATGATTTACCGGCATGTTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTGAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              <  1:1106381/149‑1 (MQ=255)
                                   gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                                   >  1:1368979/1‑150 (MQ=255)
                                     gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                >  2:428499/1‑151 (MQ=255)
                                      ttATATGATTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                <  1:109474/150‑1 (MQ=255)
                                      ttATATGAGTTACCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                <  1:129091/150‑1 (MQ=255)
                                         tatGAGTTATCGGCATGATTACCCGTTTCTGCAGGTTTTTTAATGGGACGATTAGGGGGTAAAAGGAGGGGTTATTGGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                              <  1:82726/149‑1 (MQ=255)
                                          atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                           >  2:1238759/1‑151 (MQ=255)
                                              gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                       >  2:1014633/1‑151 (MQ=255)
                                               ttATCGGCATTATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTAAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                        <  1:644159/149‑1 (MQ=255)
                                               ttAGCGGCATGATTACCGGTTTCGGCGGGGTTTTTAATCGGACGATTGGTGGGTAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                            <  1:270725/145‑1 (MQ=255)
                                                   cgtcttgtttttCCTTTTCTGCGGGTTTTTTAATGGAACTTTTAGGGGTTAAAAGGGGGGTTTATTGGGGTTTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                    <  1:958946/138‑1 (MQ=255)
                                                    ggCATATTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGGGAAATGGGGGGTTATTTGGGGGTTACAGGAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                  <  1:846194/150‑1 (MQ=255)
                                                     gCATGATTTCCCGTTTCTGCAGGTTTTTTAATCGGAGGATTAGGGGGTAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                <  1:359620/151‑1 (MQ=255)
                                                       aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                               >  2:387281/1‑150 (MQ=255)
                                                           ttATCCGTTTCGGCAGGTTTTTAAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTGGGGGTTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                            <  1:269271/149‑1 (MQ=255)
                                                           ttACCCTTTTCTGCAGGTTTTTAAATCGGACGTTTAGTGGGGAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                           <  1:1094035/150‑1 (MQ=255)
                                                             aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                        >  2:371781/1‑151 (MQ=255)
                                                              tCCTTTTCGGCAGGTTTTTTAAGCGGACGATTGGGGGGTAAAAGGAGGGTTTTTTTGGGGGTAACAGGAAAATCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                        <  2:1249036/150‑1 (MQ=255)
                                                               ccGTTTCTGCGGGTTTTTTAAGCGGACGATTGGGGGGTGAAATGGGGGTTTATTGGGGGGTTACAGATAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                         <  1:340381/148‑1 (MQ=255)
                                                                cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCCTTATGTTGCCGGATGCGGCGt                                                                                      >  2:1268469/1‑150 (MQ=255)
                                                                  tttCTGGAGGGTTTTTAATCGGACGTTTAGTGGGTAAATTGGGGGTTTTTTTGGGGGTTACAGAAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                       <  2:1475927/147‑1 (MQ=255)
                                                                  tttCTGCGGGTTTTTAAATCGAGCGATTAGGGGGAAAAATGAGGTTTTTTTTGGGGTTTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                    <  2:990340/150‑1 (MQ=255)
                                                                  tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                    >  2:796919/1‑150 (MQ=255)
                                                                    tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                  >  2:408032/1‑150 (MQ=255)
                                                                       gcggggTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTTTTTGGGGTTTACAGTAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                 <  2:1257190/145‑1 (MQ=255)
                                                                       gagggttttttAATGGGACGATTAGTGGGTAAAAGGGGGGTTTTTTGGGGGGTTACAGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                                <  1:1185186/143‑1 (MQ=255)
                                                                       gCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                              >  2:1279577/1‑151 (MQ=255)
                                                                          ggttttttAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                            <  1:921534/147‑1 (MQ=255)
                                                                             tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                        >  1:902042/1‑151 (MQ=255)
                                                                             tttttAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                           <  1:1318495/148‑1 (MQ=255)
                                                                             tttttAAGCGGACTTTTAGGGGGTAAAAGGAGGGTTTTTTTGGGGTTCACAGGAAATTTCAAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                         <  1:999810/150‑1 (MQ=255)
                                                                             ttttaAATCGGAGGATTAGTGGGTAAAATAGGGGGTTATTGGGGGTTGACGGATAATTCCCAGCAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGGCGTAAACGCCtt                                                                             <  1:1429689/146‑1 (MQ=255)
                                                                               tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                      >  2:12800/1‑151 (MQ=255)
                                                                               tttAATCGGACGATTAGTGGGTAAAATGAGGGGTTATTTGGGGGTTACAGGAAAATTCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                       <  1:773757/150‑1 (MQ=255)
                                                                                  aaTCGGACGTTTAGTGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGTAAATTCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    <  2:1525645/150‑1 (MQ=37)
                                                                                  aaTCGGACGATTAGGGGGTAAAAGGGGGGGTTTTTTGGGGGTTAAAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                   <  2:879909/151‑1 (MQ=38)
                                                                                   aTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                  >  2:1191909/1‑151 (MQ=255)
                                                                                      ggACGATTAGTGGGTGAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                               <  2:1152971/151‑1 (MQ=38)
                                                                                        aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                              >  1:208267/1‑150 (MQ=255)
                                                                                        aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                              >  1:868228/1‑150 (MQ=255)
                                                                                            ttggTGGGTAAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                 <  2:617009/140‑1 (MQ=37)
                                                                                               gggggTAAAAGGGGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                        <  2:270133/147‑1 (MQ=37)
                                                                                                 gggTAAAATGAGGGTTTATTGGGGGTTAACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                    <  1:147637/151‑1 (MQ=37)
                                                                                                 gggAGAAAGGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                         <  2:946759/146‑1 (MQ=37)
                                                                                                  ggTAAAATGGGGGTTTATTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                   <  2:460501/151‑1 (MQ=37)
                                                                                                     gAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                 >  1:179285/1‑150 (MQ=255)
                                                                                                           aGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          >  1:830791/1‑151 (MQ=255)
                                                                                                           aGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          >  1:1162648/1‑151 (MQ=255)
                                                                                                            ggggTTTTTTGGGGGTTAAAGGTAAATCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          <  2:490769/150‑1 (MQ=37)
                                                                                                            ggggTTATTTGGGGGTTAGAGAAAATTTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          <  1:716452/150‑1 (MQ=37)
                                                                                                             gggTTTTTTGGGGGTTTACAGAAAAATCCCGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          <  1:41832/150‑1 (MQ=25)
                                                                                                              ggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                        >  2:344386/1‑150 (MQ=255)
                                                                                                                 tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                    >  1:1464152/1‑151 (MQ=255)
                                                                                                                   ttgggggtTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                    <  1:436657/146‑1 (MQ=37)
                                                                                                                    ttGGGGTTTACAGGAAACTCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                  <  2:426647/150‑1 (MQ=25)
                                                                                                                    ttGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                  >  1:812632/1‑150 (MQ=255)
                                                                                                                    ttGGGGGCTACAGGGAAATTCAAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGATGAGAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTCATCCGGCCTACAAAAATGAGCAAATTCAATAAATTGCAATTCa                                  >  2:1211672/1‑150 (MQ=37)
                                                                                                                    tggggggTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                  <  1:169189/148‑1 (MQ=255)
                                                                                                                      gggggTTACAGGAAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                  <  1:1501110/148‑1 (MQ=255)
                                                                                                                        gggTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                             <  2:339809/151‑1 (MQ=255)
                                                                                                                               cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                       >  1:1227693/1‑151 (MQ=255)
                                                                                                                                    aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                   >  1:618107/1‑150 (MQ=255)
                                                                                                                                    aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                   >  1:1309265/1‑150 (MQ=255)
                                                                                                                                    aaaTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  <  2:225767/151‑1 (MQ=255)
                                                                                                                                    aaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                    <  1:1190976/149‑1 (MQ=255)
                                                                                                                                     aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  >  2:217706/1‑150 (MQ=255)
                                                                                                                                      aTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  >  1:818634/1‑149 (MQ=255)
                                                                                                                                          ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc             >  1:1130009/1‑151 (MQ=255)
                                                                                                                                           cAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc             >  1:23696/1‑150 (MQ=255)
                                                                                                                                            aGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca            <  1:960007/150‑1 (MQ=255)
                                                                                                                                             ggCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca            <  2:1238258/149‑1 (MQ=255)
                                                                                                                                             ggCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag           <  2:60393/150‑1 (MQ=255)
                                                                                                                                                 gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt      >  1:781633/1‑151 (MQ=255)
                                                                                                                                                  aaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt      <  1:383131/150‑1 (MQ=255)
                                                                                                                                                   aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa    <  2:387676/151‑1 (MQ=255)
                                                                                                                                                     aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  >  2:127077/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001768‑2002065

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