Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,363,0120A.58.4% 4.2 / 265.5 127coding (206/2382 nt)yhcDputative fimbrial usher protein YhcD
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base . (44/29);  total (72/55)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 138 total aligned reads displayed.

ATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCC  >  NC_000913/3362882‑3363142
                                                                                                                                  |                                                                                                                                  
gtcacatctgtcgggttttttcccatcactccagaatctgtcatttattatgtcattgacacgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGaaa                                                                                                                                  <  2:345537‑M1/14‑1 (MQ=255)
 tcacatctgtcgggttttttcccatcactccagaatctgtcatttattatgtcattgacacgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGaaaa                                                                                                                                 <  2:307306‑M1/15‑1 (MQ=255)
  cacatctgtcgggttttttcccatcactccagaatctgtcatttattatgtcattgacacgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGaaaaa                                                                                                                                <  1:343790‑M1/16‑1 (MQ=255)
   gCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATCGTACCACTAAATGCCGaaaaaa                                                                                                                               >  1:232925/1‑133 (MQ=255)
   gCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaa                                                                                                                               >  1:57239/1‑133 (MQ=255)
   gCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaa                                                                                                                               >  1:357494/1‑133 (MQ=255)
   gCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaa                                                                                                                               >  1:286511/1‑133 (MQ=255)
   gCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaa                                                                                                                               >  1:320989/1‑133 (MQ=255)
     atctgtcgggttttttcccatcactccagaatctgtcatttattatgtcattgacacgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGaaaaaaaa                                                                                                                             <  1:68636‑M1/19‑1 (MQ=255)
             ggttttttcccatcactccagaatctgtcatttattatgtcattgacacgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACAt                                                                                                                     >  1:313345‑M1/107‑133 (MQ=255)
                 ttttcccatcactccagaatctgtcatttattatgtcattaacgcgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCt                                                                                                                 <  2:509603‑M1/31‑1 (MQ=255)
                     cGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGTG‑AAAAAAAGGCAACATTGCTAcgcg                                                                                                            <  2:344904/133‑1 (MQ=255)
                                       gtcatttattatgtcattgacacgatcagttatccctgtgtctgcaacatttcctcgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCg                                                                                           <  2:466743‑M1/53‑1 (MQ=255)
                                       gtcatttattatgtcattgacacgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGTAAAAAAAGACAAGATTGCTACGCGCATTGAAGAAAATGGCg                                                                                           <  2:190615‑M1/53‑1 (MQ=255)
                                       cTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGGAACATTGATACGCGCTCTGATGAAAATGGCGa                                                                                          >  2:157672/1‑133 (MQ=255)
                                       cTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAGCATTGCTGCGCGCTATGATGAAAATGGCGa                                                                                          >  1:159369/1‑133 (MQ=255)
                                       cTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGa                                                                                          >  2:459899/1‑133 (MQ=255)
                                          ggTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAt                                                                                       <  2:286511/133‑1 (MQ=255)
                                            ttattatgtcattgacacgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                                                                      <  2:109425‑M1/58‑1 (MQ=255)
                                                  tgtcattgacacgatcagttatccctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCtt                                                                                <  1:363391‑M1/64‑1 (MQ=255)
                                                      aaTATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCa                                                                           >  1:501774/1‑133 (MQ=255)
                                                      aaTATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCa                                                                           >  1:492687/1‑133 (MQ=255)
                                                      aaTATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATCGCTACGCGATTTTATGAAAATGGCGCCATTTGTCTTGATCa                                                                           >  2:228749/1‑133 (MQ=255)
                                                      aaTATATCGCATTATTTTTCCCGTTCAGCCCGTTTTCTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCa                                                                           >  1:258426/1‑133 (MQ=255)
                                                         gacacgatcagttattcctgtgtctgcaacatgtccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAgg                                                                         >  2:209974‑M1/63‑133 (MQ=255)
                                                         atatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGc                                                                        >  1:83385/1‑133 (MQ=255)
                                                           atCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTATGCGCTTTGTTGACAATTGCGACATTTGTCTTGATCAGGCAt                                                                      >  2:180475/1‑133 (MQ=255)
                                                           atCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAt                                                                      <  1:263039/133‑1 (MQ=255)
                                                            tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGc                                                                        <  1:24709/130‑1 (MQ=255)
                                                            tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGc                                                                        >  2:24709/1‑130 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGATAATGGCGCCATATGACTTGATGAGGCAttt                                                                    >  1:166777/1‑133 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                                                                    >  2:114069/1‑133 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                                                                    >  1:471490/1‑133 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                                                                    >  1:20211/1‑133 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                                                                    >  2:325065/1‑133 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                                                                    >  2:365961/1‑133 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                                                                    >  1:100657/1‑133 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGAACAGGCAttt                                                                    >  2:315012/1‑133 (MQ=255)
                                                             cGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGAGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                                                                    >  2:179421/1‑133 (MQ=255)
                                                              gCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTc                                                                   >  1:283885/1‑133 (MQ=255)
                                                               cATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                                                                  >  2:41907/1‑133 (MQ=255)
                                                                        tCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGc                                                         <  2:492687/133‑1 (MQ=255)
                                                                         cctgtgtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGATAATGGCGAAATTGGCCATGATCAGGTATTTGagcaacaatc                                                         <  2:249557‑M1/87‑3 (MQ=255)
                                                                         cccGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGcc                                                        >  2:420453/1‑133 (MQ=255)
                                                                            gTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGt                                                     >  1:352981/1‑133 (MQ=255)
                                                                             ttCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGACGACATTTGTCTAGATGAGGCATTTGTGCAACAAGCCGGtt                                                    >  2:136977/1‑133 (MQ=255)
                                                                              gtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACTCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAg                                                          >  2:308416‑M1/42‑127 (MQ=255)
                                                                              gtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAg                                                          <  1:308416‑M1/86‑1 (MQ=255)
                                                                              gtctgcaacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGtt                                                    <  2:509156‑M1/92‑1 (MQ=255)
                                                                                 gCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa                                                >  2:112666/1‑133 (MQ=255)
                                                                                   caacatttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACTTTTGTCTTTATCTGGCCTTTCTGCATCAAGTCGGTTTaaaa                                               <  2:149675‑M1/97‑1 (MQ=255)
                                                                                       tttttGCCTGGTGAATATCCACGGATAGTAGCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTcc                                          <  1:228749/133‑1 (MQ=255)
                                                                                       atttccacgagaagctatcattgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTCTCTTGATCAGGCATTTCTGCAACAACCCGGATTAAAAATTc                                           <  2:211612‑M1/101‑1 (MQ=255)
                                                                                            gCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTg                                     <  1:365961/133‑1 (MQ=255)
                                                                                                       atTCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAAt                          >  2:266667/1‑133 (MQ=255)
                                                                                                         tCACTGATAGTATCATTAAATGGCG‑AAAAAAATGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATgg                        <  2:232925/133‑1 (MQ=255)
                                                                                                         tCACTGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATgg                        <  1:470245/133‑1 (MQ=255)
                                                                                                           ttgttatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGc                       >  1:190615‑M1/13‑133 (MQ=255)
                                                                                                           aCTGATAGTATCAGTAAAAGGC‑AAAAAAAATGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCt                      <  1:201356/133‑1 (MQ=255)
                                                                                                             tGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGt                    <  2:336144/133‑1 (MQ=255)
                                                                                                             tGATAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGt                    <  1:114069/133‑1 (MQ=255)
                                                                                                               tatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGt                    >  2:343821‑M1/9‑132 (MQ=255)
                                                                                                               tatatctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGt                    <  1:343821‑M1/124‑1 (MQ=255)
                                                                                                               aTAGTATCTGTAAATTGC‑AAAAAAAAGGCATCATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTa                  <  1:180475/133‑1 (MQ=255)
                                                                                                               aTAGTATCAGTAAATTGCG‑AAAAACAGGGAACCTTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTa                  <  2:211661/133‑1 (MQ=255)
                                                                                                               aTAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTa                  <  2:495758/133‑1 (MQ=255)
                                                                                                               aTAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTa                  <  2:471490/133‑1 (MQ=255)
                                                                                                               aTAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTa                  <  1:420453/133‑1 (MQ=255)
                                                                                                               aTAGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTa                  >  2:40243/1‑133 (MQ=255)
                                                                                                                atatctaATTAAATGCCGAAAAGAAAGGCATCATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTa                  >  1:149675‑M1/8‑133 (MQ=255)
                                                                                                                atatctaATTAAATGCCGAAAAGAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTa                  >  1:211612‑M1/8‑133 (MQ=255)
                                                                                                                 aGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACg                >  1:206853/1‑133 (MQ=255)
                                                                                                                 aGTATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACg                <  1:315012/133‑1 (MQ=255)
                                                                                                                  atctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACg                >  2:280645‑M1/6‑133 (MQ=255)
                                                                                                                   tctaAGTAATAGGCGAAAAAAAAGCCTACATTCCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGa               >  2:227937‑M1/5‑133 (MQ=255)
                                                                                                                   tctaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGa               <  1:454747‑M1/129‑1 (MQ=255)
                                                                                                                   tATCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGAc              <  2:515270/133‑1 (MQ=255)
                                                                                                                    ccaAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGAC‑TTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACt             >  1:501496‑M1/4‑133 (MQ=255)
                                                                                                                    aTCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACt             <  2:471260/133‑1 (MQ=255)
                                                                                                                    aTCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACt             <  1:112666/133‑1 (MQ=255)
                                                                                                                    aTCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACt             <  2:20211/133‑1 (MQ=255)
                                                                                                                    aTCAGTAAATGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACt             <  1:40243/133‑1 (MQ=255)
                                                                                                                     taAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACt             >  1:307306‑M1/3‑133 (MQ=255)
                                                                                                                       AGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACtat           <  1:373842‑M1/133‑1 (MQ=255)
                                                                                                                       AGTAAATGGCGAAAAAAA‑GGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTACTGAAGAGAAAGATGGCTGTTACGACtata          <  1:151883‑M1/133‑1 (MQ=255)
                                                                                                                        gTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGATAAATATGGCTGTTTCGACtata          >  1:227937/1‑133 (MQ=255)
                                                                                                                        gTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACtata          >  1:433785/1‑133 (MQ=255)
                                                                                                                        gTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTACTGAAGAAAAAGATGGCTGTTACGACtata          >  1:145824/1‑133 (MQ=255)
                                                                                                                          aaaTGGCGAAGAAAAGGGCAACATGGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACtatata        <  2:203344/133‑1 (MQ=255)
                                                                                                                           aaTGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTGCGACTATATAtt      >  1:307105/1‑133 (MQ=255)
                                                                                                                           aaTGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACt             <  1:41074/126‑1 (MQ=255)
                                                                                                                            aTGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATAtt      >  1:4500/1‑133 (MQ=255)
                                                                                                                            aTGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATAtt      >  2:79043/1‑133 (MQ=255)
                                                                                                                            aTGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGGAAACAATGGCTGATACGACTATATATTg     >  2:156924/1‑133 (MQ=255)
                                                                                                                             tGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAACTGGCTGGTACGACTATATATTg     >  1:203344/1‑133 (MQ=255)
                                                                                                                             tGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTg     >  2:74356/1‑133 (MQ=255)
                                                                                                                             tGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGCCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTg     >  1:437935/1‑133 (MQ=255)
                                                                                                                             tGGCG‑AAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAACGGCAGTAACGACTCTATATTGt    >  1:183159/1‑133 (MQ=255)
                                                                                                                               gCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTc   >  2:11256/1‑133 (MQ=255)
                                                                                                                                cGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTcc  <  1:411231/133‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                  
ATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCC  >  NC_000913/3362882‑3363142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑