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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065
                                                                                                                                                       |                                                                                                                                                   
aGCGGTATAGAGTGAATTCACGATAAACAGACGTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTACTTGGGGGTTACAGGTAAATTCAAGGCAGaaaaaaa                                                                                                                                                     >  1:176724/1‑150 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                   <  2:62855/151‑1 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                                 <  1:718945/151‑1 (MQ=255)
    gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                                >  1:1081878/1‑151 (MQ=255)
                 tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGTAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                                     >  1:143281/1‑149 (MQ=255)
                 tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                    >  1:127898/1‑150 (MQ=255)
                   aCGAAAAAAAGCCCTGCGTTATATGAGTTATCGGCATGTTTATCTGTTTCGGGAGGTTTTTAAGGCGGAAAATTAGGGGGTAAAGGGAGGGTTTTTTGGGGGGTTAAAGTAAACCCCCGGGGAGAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                  <  2:325081/150‑1 (MQ=255)
                      aaaaaCAGCCCTGCTTTATATGAGTTCTCGGCTTGTTTATCCTTTTCGGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              <  2:138636/149‑1 (MQ=255)
                        aaaCACCCCTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGAGGATTGGTGGGTAAAAGGAGGGGTTATTGGGGGGTTAAAGGTAATTTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              <  2:472787/149‑1 (MQ=255)
                         aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAAGtt                                                                                                                           >  1:174487/1‑151 (MQ=255)
                              cccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGCAGCACGTTGc                                                                                                                         >  1:719479/1‑148 (MQ=255)
                               ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                       >  2:424728/1‑149 (MQ=255)
                               ccTGCGTTATATGAGTTATCGGAATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGCGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                                         >  1:870877/1‑147 (MQ=255)
                                cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                     >  2:1159717/1‑150 (MQ=255)
                                  gCGTTATATGAGTTATCGGCATGATTATCCTTTTCTGCGGGTTTTTTAATGGGACGATTAGGGGGTGAAATGAGGGTTTATTGGGGGTTTACAGTTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                    <  2:1031216/149‑1 (MQ=255)
                                     ttATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGGTTTTTAAGCGAACGATTAGGGGTTGAAAGGGGGGTTTATTTGGGGGTTACAGGAAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                                <  2:513929/150‑1 (MQ=255)
                                        tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                             <  2:1078212/150‑1 (MQ=255)
                                          tGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                          >  2:939207/1‑151 (MQ=255)
                                              ttATCGGCATGATTATCCTTTTCTGCAGGGTTTTTAATCGGAGGATTAGTGGGTGAAATGAGGGTTTTTTTGGGGTTTACGGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                        <  2:1031405/149‑1 (MQ=255)
                                                aTCGGCATGATTATCCTTTTCGGCGGGTTTTTTAAGCGGACGATTGGTGGGTAAAATGAGGGGTTATTTGGGGGTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                     <  2:1153326/150‑1 (MQ=255)
                                                   ggCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGTTTTTTTGGGGTTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                                  <  2:76762/150‑1 (MQ=255)
                                                          ttATCCGTTTCTGCGGGGTTTTAAAGCGGACGATTAGTGGGTAAAAGGAGGGGTTATTTGGGGTTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                           <  2:942051/150‑1 (MQ=255)
                                                           tATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGTGGGTGAAAGGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                          <  2:1119930/150‑1 (MQ=255)
                                                           tATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTTTTTGGGGGTTAGAGGTAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                          <  1:397482/150‑1 (MQ=255)
                                                             tCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                       >  2:332409/1‑151 (MQ=255)
                                                                 tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                    >  1:361570/1‑150 (MQ=255)
                                                                  ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                  >  1:577365/1‑151 (MQ=255)
                                                                  ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                  >  2:1088529/1‑151 (MQ=255)
                                                                     tGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGATGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                               >  2:1228114/1‑151 (MQ=255)
                                                                      gcgggttttttAATCGGACGATTAGTGGGTGAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                 <  2:755384/142‑1 (MQ=255)
                                                                      gCAGGTTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                               <  2:728921/150‑1 (MQ=255)
                                                                       cAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                             >  1:365986/1‑151 (MQ=255)
                                                                         ggttttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAGTTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                             <  1:77974/146‑1 (MQ=255)
                                                                         ggtttttaAATCGGACGATGAGTGGGTGAAAGGGGGGTTTTTTGGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                               <  2:614988/144‑1 (MQ=255)
                                                                         gggTTTTTAATCGGACGATTGGGGGGTAAAATGGGGGGTTATTGGGGGGTTACAGGTAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                             <  2:826605/149‑1 (MQ=255)
                                                                         gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                  <  1:54100/144‑1 (MQ=255)
                                                                          ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                          >  2:1211499/1‑151 (MQ=255)
                                                                             ttttAATCGGACGCTTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGAGGCGTAAACGCCTTATcc                                                                         >  1:119741/1‑149 (MQ=255)
                                                                             ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                       >  1:296653/1‑151 (MQ=255)
                                                                                tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACCAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGACt                                                                    >  1:281844/1‑151 (MQ=255)
                                                                                tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    >  2:1189993/1‑151 (MQ=255)
                                                                                 atgcGGACGATGAGGGGGTAAAATGGGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                   <  1:1173005/148‑1 (MQ=38)
                                                                                 aaTCGGACGTTTAGGGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                     <  2:722890/149‑1 (MQ=255)
                                                                                 aaTCGGACGATTAGTGGGTGAAATGAGGGGTTTTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                   <  1:223195/151‑1 (MQ=255)
                                                                                 aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    >  1:989636/1‑150 (MQ=255)
                                                                                       atatttAGTGGTTGAAATGAGGGTTTATTTGGGTGTTTGAGTGAAATTCCAGGCAAAAAAAACCTCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTATATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACaaaaa                                                             <  1:66513/147‑1 (MQ=37)
                                                                                        cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                            >  2:274199/1‑151 (MQ=255)
                                                                                           ttggggggTAAAAGGGGGGTTTTTTTGGGGGTTACAGTTAATTCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                           <  1:1139973/144‑1 (MQ=25)
                                                                                            gagggggTAAAATGAGGGTTTTTTTGGGGGTTACAGGTAATTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                         <  1:614197/146‑1 (MQ=37)
                                                                                             aGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                        >  2:142945/1‑150 (MQ=255)
                                                                                              gggggTGAAAGGGGGGGTTATTGGGGGGTTAAAGGTAAATCCCAGGAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                       <  2:274976/148‑1 (MQ=38)
                                                                                                gggTGAAAGGGGGGGTTTTTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                         <  1:15518/146‑1 (MQ=38)
                                                                                                     aaaTGAGGGGTTTTTTGGGGGGTACAGTAAAATCCCGGGAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCGCGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               <  1:631345/151‑1 (MQ=37)
                                                                                                     aaaGGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGGAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               <  1:705655/151‑1 (MQ=255)
                                                                                                     aaaGGAGGGGTTATTTGGGGGTAACAGTAAATTACCAGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               <  1:294053/151‑1 (MQ=25)
                                                                                                      aaTGGGGGTTTTTTTGGGGGTTACAGTAAATTCCGAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                              <  1:884024/151‑1 (MQ=25)
                                                                                                       aTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                              >  1:106040/1‑150 (MQ=255)
                                                                                                         gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                           >  2:193217/1‑151 (MQ=255)
                                                                                                          aGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          >  2:780463/1‑151 (MQ=255)
                                                                                                                  tttGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGAACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                   >  2:353655/1‑150 (MQ=255)
                                                                                                                   ttGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                  >  2:203265/1‑150 (MQ=255)
                                                                                                                   ttGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                  >  1:1134003/1‑150 (MQ=255)
                                                                                                                     gggtgttaaggttAATTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                 <  2:361649/137‑1 (MQ=25)
                                                                                                                          ttatggtgaagttCCAGCCAGAAAACACCCCCGCCGGTGGCTGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                            <  2:996846/138‑1 (MQ=25)
                                                                                                                          ttACAGATAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                           <  1:977061/150‑1 (MQ=255)
                                                                                                                           tACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                         >  2:301097/1‑151 (MQ=255)
                                                                                                                            aCAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                         <  1:803396/150‑1 (MQ=255)
                                                                                                                                ggaaagtcccAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                        <  2:573808/140‑1 (MQ=37)
                                                                                                                                 tAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                   >  1:791695/1‑151 (MQ=255)
                                                                                                                                  aaaTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                        <  1:86377/145‑1 (MQ=255)
                                                                                                                                  aaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                   <  1:938946/150‑1 (MQ=255)
                                                                                                                                  aaaTACCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                       <  2:36937/146‑1 (MQ=255)
                                                                                                                                   aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAACACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACACCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  >  1:1219197/1‑150 (MQ=255)
                                                                                                                                     ttCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag                <  1:320484/150‑1 (MQ=255)
                                                                                                                                               aGAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg         <  1:757421/149‑1 (MQ=255)
                                                                                                                                                 aaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc          <  2:1185638/146‑1 (MQ=255)
                                                                                                                                                   aaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg   <  1:843082/151‑1 (MQ=255)
                                                                                                                                                    aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  >  2:1176161/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065

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