Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	753,125	0	T	A	43.3%	79.0 / 34.3	60	intergenic (‑330/+60)	ybgD/gltA	putative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (29/5);  new base A (0/26);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 2.48e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCAT  >  NC_000913/752982‑753274                                                                                                                                                |                                                                                                                                                      cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                                >  1:608012/1‑151 (MQ=255) cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                                <  2:420370/151‑1 (MQ=255)   tAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCa                                                                                                                                               >  1:271583/1‑150 (MQ=255)       gTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                                          >  2:88735/1‑151 (MQ=255)              aTGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATg                                                                                                                                   >  1:338611/1‑151 (MQ=255)               tGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATg                                                                                                                                   >  1:148396/1‑150 (MQ=255)               tGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATg                                                                                                                                   >  1:339061/1‑150 (MQ=255)                     tACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcg                                                                                                                            >  2:152448/1‑151 (MQ=255)                       ccATGAATTTTTTAATTAAAAGTATTGGTTTGCTGATATTTTGGGCTTTTTTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTATTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAAcggcgg                                                                                                                           <  1:791109/150‑1 (MQ=255)                       ccATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCggg                                                                                                                          >  2:32863/1‑151 (MQ=255)                        cATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCggg                                                                                                                          >  1:409312/1‑150 (MQ=255)                         agtaattttttAATTGAAAGTTTTGGGTTGCTAATAATTTGGGTTGTTCTTTTTTTTTTAAATATCAATATAGTTCTTTAATTGAATTTTATTTTAACAATTTTTTGTGTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCggg                                                                                                                          <  1:159477/143‑1 (MQ=255)                            aattttttAATTGAAATTTTTGGTTTGCTGATAATTTGAGCTTTTTTTTTTTTTTTAAATTTTTATATAGGTTTTTTAATGGTTTTTTTTTTTACAATTTTTTGTGTTTGGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTa                                                                                                                       <  2:675688/146‑1 (MQ=255)                            aattttttAATTAAAAGTATGGGGTTGCTGATATTTTGAGTTTTTCTTTTCTTTTTAAATATCTATATAGGTTTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTa                                                                                                                       <  2:273839/146‑1 (MQ=255)                            aatttttaAATTAAAATTATTGGTTTGTTAATATTTTGAGCTTTTTTTTTTTTTTTAAATTTTTATATAGGTCTTTTAAGGAATTTTTTTTTTAAAATTTTTTGTTTTAGGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTaa                                                                                                                      <  1:90999/147‑1 (MQ=255)                            aaGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTaa                                                                                                                      >  2:275002/1‑150 (MQ=255)                              tttttttATAAAAATTTTTGTGTTGCTGATATTTTGCTTTTTTTTTTTTTTTTTAAATATTTATTGAGTTTTATAAATGTTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAACCACCCCGCCATATGAAcggcgg                                                                                                                           <  1:594640/142‑1 (MQ=255)                               tttttAATTAAAATTTTTGGTTTGCTAAAAATTTGTGCTTTTTTTTTTTTTTAAAATATTTAAATAGGTCTTTAAATGGTTTTTTTTTTTAAAATTTTGTTTTTTAAGGCTAAAAAAACCCAACCCGCCATATGAACGGCGGGTTAAAat                                                                                                                   <  1:637365/150‑1 (MQ=255)                                  ttaattaaAAGTTTTGGGTTGTTGATAATTTGAGGTTTTTTTTTTTTTTTAAATATATATATGGGTCTTTTAATGGATTTTTTTTTTAAAATTTTTTGTGTTTAGGAATAAAAAAATAACCCCGCCATATGAACGGCGGGTTAAAATAtt                                                                                                                <  2:931265/150‑1 (MQ=255)                                  ttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                               >  1:697858/1‑151 (MQ=255)                                       tGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTATAAAATTTTCa                                                                                                            >  2:246573/1‑149 (MQ=255)                                                tGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                                  >  1:4126/1‑150 (MQ=255)                                                 ggttaggtGATAATGTGAGCTCTTCAATTTTTATTAAATTTTAATATGGGTCTATTAATGATTTTTTTTTTTACATTTTTTTGTTTAAAGCGATAAAAAAAGCATCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAg                                                                                                      <  2:271603/138‑1 (MQ=255)                                                          atatTTTGAGCTGTTTTATTTTTTTAAAATATCTATATAGTTTTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCa                                                                                            <  2:672095/146‑1 (MQ=255)                                                          aTACTTTGAGCTTTTCTATTATTTTAAACTATCTATATAGGTCTGATAATGGATTTTATTTTTACATTTTTTTATGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAAGATTTACCAATAAGCaaaccaaccttta                                                                                        >  1:804660/1‑139 (MQ=255)                                                          aTAATTTGAGCTGTTCTATTCTTTTAAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                                                          <  2:239059/148‑1 (MQ=255)                                                            aaTTTGAGCTTTTTTATTCTTTTAAAATATCTATATAGGTCTGTAAATGGATTTTATTTTTAAAATTTTTTGTTTTTAGGCATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTaa                                                                                        <  2:703565/148‑1 (MQ=255)                                                            aaTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTaaacc                                                                                     >  1:430327/1‑148 (MQ=255)                                                             aTTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAAAATTTACAACTTAGCAATCAAACATTaaag                                                                                      >  2:794948/1‑147 (MQ=255)                                                              tttGAGCTTTTTTTTTTTTTTTAAATATTTATATAGGTCTGTAAATGGTTTTTTTTTTTACATTTTTTTTTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAc                                                                                       <  2:897294/147‑1 (MQ=255)                                                                tGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCtt                                                                                   >  2:303042/1‑149 (MQ=255)                                                                   gttttttttttCTTTTAAAATATTTATATAGGTTTTTTAATGTATTTTTTTTTTACATTTTTTTTGGTTGAGACTAAAAAAACACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGatat                                                                              <  1:742979/142‑1 (MQ=255)                                                                    cTGTTTTATTTTTTTAAAATATTTAAATAGTTTTGTTAAGGGATTTTTTTTTTAAATTTTTTTGTTTTAAGACATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                                                                <  2:8546/148‑1 (MQ=255)                                                                                    aaaTATCTATATAGGTCTGTTAATGGTTTTTTTTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                             <  2:317258/151‑1 (MQ=255)                                                                                     aaTATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGtt                                                            >  1:567575/1‑151 (MQ=255)                                                                                        atCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACCCTTTATATCGCTTTTAAAGTCGCGttttt                                                         >  1:127538/1‑151 (MQ=255)                                                                                         tCTATATGGGTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTTTTTTTAGGAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttt                                                           <  1:165598/148‑1 (MQ=255)                                                                                               taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                                  >  1:210286/1‑151 (MQ=255)                                                                                                 ggtttttaaaTGAATTTTTTTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                                                <  2:281792/143‑1 (MQ=255)                                                                                                 ggTCTGTTAATGGTTTTTATTTTTAAATTTTTTTGTGTTAAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                                     <  2:143976/146‑1 (MQ=255)                                                                                                 ggTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                                  >  2:571971/1‑149 (MQ=255)                                                                                                   tttGTTAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGGTATATAAAATGCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                                                <  2:325777/147‑1 (MQ=255)                                                                                                     tGTTAAGGGTTTTTATTTTTAAATTTTTTGGTTTTAAGGCATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtat                                             <  1:245500/150‑1 (MQ=255)                                                                                                       ttAATGTTTTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                                            <  2:399629/149‑1 (MQ=255)                                                                                                         aaTGTATTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACa                                          <  1:152433/149‑1 (MQ=255)                                                                                                            ggATTTTTTTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACa                                          <  1:149602/146‑1 (MQ=255)                                                                                                               ttttattttTACCATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAAAACTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgc                                  >  2:175493/1‑151 (MQ=255)                                                                                                               ttttattttTACAATTTTTTGTTTTTAGGCATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAc                                    <  2:104375/149‑1 (MQ=255)                                                                                                               ttttattttTAAAATTTTTTGTTTTTAGACATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAc                                    <  1:715002/149‑1 (MQ=255)                                                                                                                 ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCAAATCCCGTATACAGCTGAc                                    >  2:19814/1‑147 (MQ=255)                                                                                                                    tttttACAATTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                              <  2:158750/150‑1 (MQ=255)                                                                                                                     ttttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACgcgggcgg                            >  1:674608/1‑151 (MQ=255)                                                                                                                         acttttttttGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa                         <  1:389161/146‑1 (MQ=255)                                                                                                                         aCATTTTTTTGTGTTTAGCCATAAAAAATCCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                          <  2:263106/149‑1 (MQ=255)                                                                                                                          cAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTTACGCTTGATATCGCTTTTAAAATCGCGTTTTTTATATCCTGTATACAGCTGACGCGGACGGGCaa                        >  1:94880/1‑150 (MQ=255)                                                                                                                           aaTTTTTTGTTTTTAGGCATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATc                      <  1:365678/151‑1 (MQ=255)                                                                                                                                tttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAAAATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGACGCGTTTTTCCTATCCTCTATACAGCTGACACGGACGGGCAATCtt                    >  1:937735/1‑148 (MQ=255)                                                                                                                                 ttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAAAATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAAAGCTGACGCGGACGGGGAATCTTCAt                 >  2:271275/1‑150 (MQ=255)                                                                                                                                   gtgtTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAAACTTCATAcc              >  1:499305/1‑151 (MQ=255)                                                                                                                                       ttAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGTTTTTCATATACTGTATACAGCTGACGCGGACTGGCAATATTCATACCggc           >  2:870997/1‑148 (MQ=255)                                                                                                                                              tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGc    >  2:642672/1‑150 (MQ=255)                                                                                                                                              tataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGc    >  2:38065/1‑150 (MQ=255)                                                                                                                                               ataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAAAGCTTGATATCGGTTTTAAACAGGCGCTTTTCCAATCCATTATACATCTTAAGCGGGCGGGCAATCTTCCTACCGTTACTGTGctt  >  2:421708/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      CATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCAT  >  NC_000913/752982‑753274

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