Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,005,037	C→T	6.3%	T389I (ACC→ATC)	fliD →	flagellar filament capping protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,005,037	0	C	T	6.3%	171.6 / 2.4	64	T389I (ACC→ATC)	fliD	flagellar filament capping protein
Reads supporting (aligned to +/- strand):  ref base C (35/24);  new base T (2/2);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

AGATTGGTATCACGACCGATCCCAGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATTGATGCGCAGGTCGCTCGCTAC  >  NC_000913/2004899‑2005179                                                                                                                                            |                                                                                                                                               aGATTGGTATGACGACCGATCCCAGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATc                                                                                                                                            <  1:615614/143‑1 (MQ=255)     tGGTATCACGACCGATCCCAGCGATGGCAAACTGGCACTGGATGCCGACAAACTCCCCGCTGCACTGAAAAAAGATGCCAGCGGCGTCGGTG‑CATTGATTGTTGGCGATGGGAAACAAACCGGCATCACGACCACCATCGgca                                                                                                                                        >  1:806849/1‑143 (MQ=255)       gTATCACGACCGATCCCAGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAACGACGCGTCGGGACCGCCCTCGgcagc                                                                                                                                      >  2:6646/1‑143 (MQ=255)        taccACGACCGAGCCCAGCAATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTAAAAAAAGATGCAAGCGGCGTAGGTG‑CATTGATTGTTGGCGATAGTAAAAAAACCGGCATCACGACCACCATCGgcagca                                                                                                                                     <  2:657553/140‑1 (MQ=255)        tATCACGACCGATCCCAGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGGTTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGgcagca                                                                                                                                     >  1:547624/1‑143 (MQ=255)        tATCACGACCGATCCCAGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCATCATCGgcagca                                                                                                                                     >  1:1241049/1‑143 (MQ=255)        tATCACGACCGATCCCAGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCATCATCGgcagca                                                                                                                                     >  1:897481/1‑143 (MQ=255)            aCGATCAATCCCTGCGATGGTAAACTGGAATTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGTCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCt                                                                                                                                 <  2:1071468/143‑1 (MQ=255)            aCGACCGATCCCAGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCt                                                                                                                                 >  1:1228156/1‑143 (MQ=255)            aCGACCGATCACAGCGATGGAAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATGGCTGGTTGGCGATGGTAAAAGAACAGGCATCACGACCACCACCGGCAGGaaaca                                                                                                                                 >  1:1235839/1‑140 (MQ=255)                    tCCCAGCGATGGCAAACTGGAACTGTATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGtt                                                                                                                         <  1:465096/143‑1 (MQ=255)                        aGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGcg                                                                                                                     >  1:724026/1‑142 (MQ=255)                        aGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCt                                                                                                                     >  1:1092915/1‑143 (MQ=255)                          cGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCttt                                                                                                                   >  1:822526/1‑143 (MQ=255)                           gATGTCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTc                                                                                                                  <  1:1084382/143‑1 (MQ=255)                            aTGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAACAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGCCATCACGACCACCATCGGCAGCAACCTGCCCAGTTGGCTCTCg                                                                                                                 >  1:190797/1‑143 (MQ=255)                             tGGCAAACTGGAACTGGAGGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGa                                                                                                                <  1:689063/143‑1 (MQ=255)                                 aaaCTGGAACTGGATGCCGGCAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAAc                                                                                                            >  1:662700/1‑143 (MQ=255)                                  aactggaactggATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTCTCGACaccc                                                                                                           >  2:708945/1‑140 (MQ=255)                                  aactggaactggATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACAACCACCATCGGCAGCAACCTGACCAGTTGGCACTCGACAACg                                                                                                           >  2:347077/1‑143 (MQ=255)                                    ctggaactggaTGCCGACAAACTCACCGCTGCCCTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGCCTGTTGGCGATGGTAAAAAAACCGGCATCACGCCCCCCATCGGCAGCACCCTGACCAGTCGACCTTCGCCAACggg                                                                                                         >  2:691061/1‑143 (MQ=255)                                     tggaacaggaTGCCGACAAACACACCGCCGCACTTAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGCTGCGAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGc                                                                                                        <  1:969439/143‑1 (MQ=255)                                       gaactggaTGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACAGGCat                                                                                                      >  2:1140583/1‑143 (MQ=255)                                               tGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGCCGAGTATTCACTCGACAACGGGCATTATTAAAg                                                                                              >  1:173740/1‑143 (MQ=255)                                               tGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAg                                                                                              >  2:120642/1‑143 (MQ=255)                                                  cGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCg                                                                                           >  1:281917/1‑143 (MQ=255)                                                  cGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATCGATTGTTGGCGATGGTAAACAAACCGGCATCACGACCACAATCGGCCGCAACCTGGCCAGTTGGCTTCCGCCCCCGGGGATTATCAACGCCg                                                                                           >  2:1048505/1‑143 (MQ=255)                                                   gACAAACTCGCCGCTGCACTGAAAAAAGATGCCAGCGGCGT‑GATGTCATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCATCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGc                                                                                          <  2:1241049/143‑1 (MQ=255)                                                   gACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCATCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGc                                                                                          <  2:897481/143‑1 (MQ=255)                                                    aCAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGGTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCt                                                                                         >  1:1077928/1‑143 (MQ=255)                                                          tCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGa                                                                                                                               <  1:711783/99‑1 (MQ=255)                                                          tCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGa                                                                                                                               >  2:711783/1‑99 (MQ=255)                                                            aCCGCTGCACTGAAAAAAGATGCCAGCGGCGTCGGTC‑CATTCCTTGTTGACGATGATAAAAAAAAAAACATAACGACAAAAAACGGCAGCAACATGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACGGATgg                                                                                 >  1:17550/1‑143 (MQ=255)                                                              ggCTGCCCGGAAAAAAGACGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTCAAAAAAACGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCTACAACGGGCATTATTAAAGCCGCTACCGATGGCg                                                                               <  2:325296/142‑1 (MQ=255)                                                                 tGCACTGAGAAAACGTGCCAGCGACGTGCCTT‑CATTCATTATTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCCCTACCGATGGCGTTa                                                                            <  2:1078593/143‑1 (MQ=255)                                                                     cTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGGTAGTaa                                                                        >  1:97612/1‑143 (MQ=255)                                                                        aaaaaaGATGCCAGCCGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGAc                                                                     <  1:177116/143‑1 (MQ=255)                                                                           aaaGATGCCAGCGGCGTAGGTG‑CATTGCTTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCGGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCt                                                                  <  1:508574/143‑1 (MQ=255)                                                                           aaaGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCt                                                                  <  1:154480/143‑1 (MQ=255)                                                                            aggATGCCGGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCACCCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTg                                                                 <  1:56091/141‑1 (MQ=255)                                                                            aaGATGCCAGCGGCGTAGGTG‑CATTGAGTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTg                                                                 <  1:1012786/143‑1 (MQ=255)                                                                                 gCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCCGTTGGGTTTCGACAACGGGCATTCTTAAAGCCGCTACCGATGGCGTGAGTAACACCCTGAACaa                                                            >  2:447241/1‑143 (MQ=255)                                                                                   cAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCCCGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAACACCCTGTATAAAt                                                          >  1:426075/1‑143 (MQ=255)                                                                                    aGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAAtt                                                         <  1:1189016/143‑1 (MQ=255)                                                                                        gcgTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGaa                                                               >  2:188085/1‑133 (MQ=255)                                                                                        gcgTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGaa                                                               <  1:188085/133‑1 (MQ=255)                                                                                           tAGGTG‑CATTGGTTGGTGGCGAAGGTAAAAATGCCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTaaa                                                  <  2:376958/143‑1 (MQ=255)                                                                                             ggTG‑CATTGATTGTTGGCGATGGTACAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAgaga                                                <  1:991846/143‑1 (MQ=255)                                                                                              gTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGAc                                               >  2:366529/1‑143 (MQ=255)                                                                                                   attgattgTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACa                                           >  1:817267/1‑143 (MQ=255)                                                                                                    ttgattgTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCGTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACaa                                          >  2:644798/1‑143 (MQ=255)                                                                                                     tgattgTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAAACTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAAc                                         >  2:1054377/1‑143 (MQ=255)                                                                                                        ttgttgGCGAGGGTAAAAATTCCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAAcgcc                                      <  2:758203/143‑1 (MQ=255)                                                                                                           ttgGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGAGCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGACGGCGTCAGTAAGACCCTGAATAAATTAACCAAAGACTACCAAgccgcc                                   >  2:1148898/1‑143 (MQ=255)                                                                                                           ttgGCCATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGAACAGTAGGCTTTCGACATCGGGCATGATTAAAGCCGCTACAGATGGCGGTCGTAAGACCCGGAAGAAATTAGCTAAAGACTACAACgcagcc                                   >  2:1064391/1‑143 (MQ=255)                                                                                                            cggCGATCGTAAAAGAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCa                                  <  2:737517/142‑1 (MQ=255)                                                                                                                    gTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGc                          >  2:426274/1‑143 (MQ=255)                                                                                                                       aaaaaaCCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCAtt                       >  2:1265831/1‑143 (MQ=255)                                                                                                                       aaaaaaCCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGAACCTGAATAAATTAACTCAAGACTACACCGCCGCCAGCGATCGCaga                       >  2:60633/1‑141 (MQ=255)                                                                                                                              cGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATTGATgcgc                >  1:746421/1‑143 (MQ=255)                                                                                                                              cGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATTGATgcgc                <  1:334290/143‑1 (MQ=255)                                                                                                                                gCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTACCTAAAGACTAAAACGCCGCCAGCGATCGCATTGATGCGCaa              >  2:78952/1‑142 (MQ=255)                                                                                                                                  aTCACGATCACCATCGTCAGCATCCTTACCACTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATTGATGCGCAGGt            <  1:95033/143‑1 (MQ=255)                                                                                                                                      cGACCACCATCGGCAGCACCCGGACCAGGTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATTGATGCGCAGGtcgct        <  2:1248620/143‑1 (MQ=255)                                                                                                                                          caccaTCGGAAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATTGATGCGCAGGTcgctcgct    <  1:347077/143‑1 (MQ=255)                                                                                                                                            ccaTCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATTGATGCGCAGGTCGCTCGCGAc  >  2:1073271/1‑143 (MQ=255)                                                                                                                                            ccaTCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGCTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATAGATGCGCAGGTCGCCCGCTAc  >  2:725032/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                               AGATTGGTATCACGACCGATCCCAGCGATGGCAAACTGGAACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTG‑CATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACCGATGGCGTTAGTAAGACCCTGAATAAATTAACTAAAGACTACAACGCCGCCAGCGATCGCATTGATGCGCAGGTCGCTCGCTAC  >  NC_000913/2004899‑2005179

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