Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,568,504	0	T	G	33.8%	142.4 / 13.9	81	T6P (ACC→CCC)	mdtM	multidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base T (11/42);  new base G (27/0);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 7.83e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAAC  >  NC_000913/4568367‑4568643                                                                                                                                          |                                                                                                                                                agGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑aa                                                                                                                                              >  1:929148/1‑140 (MQ=255)  ggggCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑aaa                                                                                                                                             <  1:13216/140‑1 (MQ=255)          aCTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGtc                                                                                                                                     <  1:1288603/140‑2 (MQ=255)          aCTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTg                                                                                                                                     <  2:1071900/140‑1 (MQ=255)             gACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCa                                                                                                                                  <  1:718419/140‑1 (MQ=255)               cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACGTGGCaaa                                                                                                                                >  1:562217/1‑138 (MQ=255)                    gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACGTGGCAT‑aacaaa                                                                                                                          >  1:165450/1‑140 (MQ=255)                    gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAATGGGAT‑aaaaaa                                                                                                                          >  2:1118547/1‑140 (MQ=255)                      aTTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGCAT‑aacaaaaa                                                                                                                        >  1:382827/1‑138 (MQ=255)                        tAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAaacaac                                                                                                                      <  1:1017971/140‑1 (MQ=255)                          aaaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTc                                                                                                                    <  2:427532/140‑1 (MQ=255)                          aaaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGCAT‑Aaaaaaaaaaac                                                                                                                    >  2:670673/1‑137 (MQ=255)                           aaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTcc                                                                                                                   <  2:348420/140‑1 (MQ=255)                               aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAATGGGCT‑AAAAaaaaaccactttt                                                                                                               >  2:859513/1‑133 (MQ=255)                                cGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTc                                                                                                              <  2:929148/140‑1 (MQ=255)                                  tACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCg                                                                                                            <  1:502377/140‑1 (MQ=255)                                       agttAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑Gaa                                                                                                      <  1:178542/138‑1 (MQ=255)                                               aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACg                                                                                             >  1:669478/1‑140 (MQ=255)                                               aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAAAGGCAT‑AACAAAAAACACCCTTTCCGCG‑GAACCCCT‑ACg                                                                                             >  2:1263846/1‑140 (MQ=255)                                                 cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAAAGGTGACT‑AACAAAAAAATACGTTTTCGGG‑GAAACTCA‑ACggg                                                                                           >  1:1289465/1‑140 (MQ=255)                                                 cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGGAT‑AAAAAACAAAACCTTGGTCTGA‑GAACCCCA‑ACggg                                                                                           >  2:1202243/1‑140 (MQ=255)                                                  ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAAAAAGGGAT‑AAAAAAAAAATCCTTGTCCCGG‑GAACCCCA‑ACgggg                                                                                          >  2:1061902/1‑140 (MQ=255)                                                  ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACGTGGCAT‑AAAAAAAAACTCCTTGTCCGGA‑GGACCCCA‑ACgggg                                                                                          >  2:723912/1‑140 (MQ=255)                                                    aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAATGGCAT‑AAAAAAAAAAACCTTGTCCCGA‑GGACCCCC‑AAggggcc                                                                                        >  2:178542/1‑138 (MQ=255)                                                      gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGaaaa                                                                                      <  2:669478/140‑1 (MQ=255)                                                      gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGaaaa                                                                                      <  2:347635/140‑1 (MQ=255)                                                      gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGaaaa                                                                                      <  1:1072531/140‑1 (MQ=255)                                                      gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGaaaa                                                                                      <  1:1268633/140‑1 (MQ=255)                                                      gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAATGGCAT‑AAAAAAAAACTCCTTTTCCGGA‑GAACCCCA‑ACGGGGaaaa                                                                                      >  1:80269/1‑140 (MQ=255)                                                        tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACg                                                                                    <  1:1202243/140‑1 (MQ=255)                                                          gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAAGGGGGAA‑AAAAAAAACCCCCTGTTCCGGG‑GAACCCAA‑ACGGGGaaaaaaa                                                                                   >  2:1268633/1‑137 (MQ=255)                                                          gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGCAT‑AAAAAACAAATCCTTGTCCGGA‑GAAACCCA‑AACGGGAAAAAGaa                                                                                  >  2:273575/1‑140 (MQ=255)                                                          gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGAGGCAT‑AAAAAAAAACACCCTGTCCCGG‑GAACCCCC‑ACGGGGGAAACGaa                                                                                  >  2:854760/1‑140 (MQ=255)                                                           aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCGATCGGGAAAAGCAGCGTGGCGGCAGGGGGGGA‑AAAAAAAAGGGGGGTAAAAAAAAAACCCCTTTCCGGGG‑GAGCCGCA‑AGGGGGaaagggaa                                                                                  >  2:926451/1‑135 (MQ=255)                                                           aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGaaa                                                                                 <  1:486452/140‑1 (MQ=255)                                                           aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGaaa                                                                                 >  2:562982/1‑140 (MQ=255)                                                           aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGGCT‑AACAAAAAACTCCCTTTCCGGAGGAACCCCA‑ACGGGG‑AAACGaaa                                                                                 >  2:486452/1‑140 (MQ=255)                                                           aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGCT‑GAAAAAAAGCGGGAT‑AACAAAAAACTCCCTGTCCAGA‑GAAACCCACAAGGGGaaaagaaa                                                                                  >  1:958086/1‑136 (MQ=255)                                                                aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCaa                                                                            >  1:345966/1‑140 (MQ=255)                                                                aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCaa                                                                            >  1:1254640/1‑140 (MQ=255)                                                                  ccGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATa                                                                          <  2:993549/140‑1 (MQ=255)                                                                    gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACg                                                                                    <  2:715495/128‑1 (MQ=255)                                                                    gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATAcc                                                                        <  1:653168/140‑1 (MQ=255)                                                                    gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAAAGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACg                                                                                    >  1:715495/1‑128 (MQ=255)                                                                    gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGC‑AAAAAAACGCGGCAT‑AACAAACAAACACTTGTCCGGA‑GAACCCCC‑ACACGGAAAAAGAAAATCAATAcc                                                                        >  2:1201416/1‑140 (MQ=255)                                                                     tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGTGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCg                                                                       >  1:210268/1‑140 (MQ=255)                                                                     tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCg                                                                       >  1:5108/1‑140 (MQ=255)                                                                     tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCg                                                                       >  2:977213/1‑140 (MQ=255)                                                                     tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACGTGGCAT‑AACAAACAACTCCCTTTCCGGG‑GAACCCCA‑ACGGGGGAAAAGAAAATCAATACCg                                                                       >  1:1190379/1‑140 (MQ=255)                                                                     tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACAAGGCGT‑AAAAAAAAACCCCTTGTCCGGG‑GGACCCCA‑ACAGGGAAAAAGAAATTCAAAAccc                                                                       >  2:615093/1‑139 (MQ=255)                                                                      cGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGt                                                                      <  2:345966/140‑1 (MQ=255)                                                                           gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTa                                                                 <  2:290295/140‑1 (MQ=255)                                                                            aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTaa                                                                <  1:910786/140‑1 (MQ=255)                                                                                 cagcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACg                                                           <  1:615093/140‑1 (MQ=255)                                                                                   gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGaa                                                         <  1:852546/140‑1 (MQ=255)                                                                                    caAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGaaa                                                        <  1:1118547/140‑1 (MQ=255)                                                                                      aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAAAGAGGCAT‑AACAAAAAACTCCTTGTCCGGG‑GAAACCCA‑ACGGGGGAAACGAAAATCAAAAACGTCGGTAAAGACGAAAgg                                                      >  2:48577/1‑140 (MQ=255)                                                                                      aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACAAGGCAT‑AAAAAAAAACTCCTTTTCCGGA‑GAACCCCA‑ACGGGGGAAAAGAAAATCAAAAACCTCGGGAATTAAGaaaag                                                      >  2:1116486/1‑138 (MQ=255)                                                                                      aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGGGGCAT‑AAAAAACAACCCCTTTTCCGGA‑GAACCCCC‑ACGGGGAAAAAGAAAATCAATACCCTCGGGAAAGACGaaaag                                                      >  2:683260/1‑138 (MQ=255)                                                                                      aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGCGGCAT‑AAAAAAAAACTCCTTGTCCGGA‑GAACCCCC‑ACGGGGGAAAAGAAAATCAATAACGTCGGGAATGAAGaaaag                                                      >  2:661585/1‑138 (MQ=255)                                                                                            tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAAAGTGGCAT‑AAAAAAAAACCCCCTTGCCGGA‑GAAACCCA‑ACGGGGAAAAAGAAAATCAATACCGTCGGGAATGACGAAAAGGAGGTg                                                >  2:1103206/1‑140 (MQ=255)                                                                                                      aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATGCCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCaa                                      <  1:1225564/140‑1 (MQ=255)                                                                                                      aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCaa                                      <  1:562982/140‑1 (MQ=255)                                                                                                      aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCaa                                      <  1:926451/140‑1 (MQ=255)                                                                                                            tCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc                                <  2:499228/140‑1 (MQ=255)                                                                                                             cGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGCTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACt                               >  1:1269449/1‑140 (MQ=255)                                                                                                                 aaaaaCAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCaa                           <  1:818819/140‑1 (MQ=255)                                                                                                                  aaaaCAGCGTGGCGGCATGGCGGGG‑AAAAAAACGGGGCAT‑AACAAAAAACTCCTTTTCCGGA‑GAACCCCA‑ACGGGGGAAAAGAAAATCAATACCGTCGGGAATGACGAAAAGCAGGGGAGAAGGGGAACGAGACTGCaaa                          >  2:670209/1‑140 (MQ=255)                                                                                                                   aaaCAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAAt                         <  1:281378/140‑1 (MQ=255)                                                                                                                       ttgctGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTAc                     <  1:160722/136‑1 (MQ=255)                                                                                                                       aGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTAc                     <  1:661585/140‑1 (MQ=255)                                                                                                                       aGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTAc                     <  2:382827/140‑1 (MQ=255)                                                                                                                       aGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTAc                     <  1:48577/140‑1 (MQ=255)                                                                                                                        gCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGCAT‑AAAAAAAAAACCCTTGTCCGGA‑GAACCCCA‑ACGGGGGAAAACAAAATCAAAACCGTCGGGAATGACGAAAAGCAGGGGAGAAAGGGAAAGAGAATGCAAATGTaat                    >  2:818819/1‑138 (MQ=255)                                                                                                                           tGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGaa                 <  1:683260/140‑1 (MQ=255)                                                                                                                                gcATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCa            <  1:1158860/140‑1 (MQ=255)                                                                                                                                gcATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCa            <  1:360394/140‑1 (MQ=255)                                                                                                                                 cATGGCGGGG‑AAAAAAAAGTGGGAT‑AAAAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGGAAACGAAATTCAATACCGTCGGGAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCAt           >  2:17237/1‑140 (MQ=255)                                                                                                                                 cATGGCGGGG‑AAAAAAAAGTGGGAA‑AACAAACAACCCCTTTTCCGGG‑GAACCCCA‑ACGGGGAAAACGAAATTTAAAACCCTCGGGAAAGACGAAAAGGAGGGGAGAAGGGCAACGGGACTGCAAAAGTAATGAATTTTAt           >  2:360394/1‑140 (MQ=255)                                                                                                                                    ggCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa        <  1:670673/140‑1 (MQ=255)                                                                                                                                    ggCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa        <  1:1078843/140‑1 (MQ=255)                                                                                                                                    ggCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa        <  1:859513/140‑1 (MQ=255)                                                                                                                                      cGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAAt      <  2:21535/140‑1 (MQ=255)                                                                                                                                          t‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAAc  <  2:80269/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                AGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAT‑AACAAACAACTCCTTGTCCGGA‑GAACCCCA‑ACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAAC  >  NC_000913/4568367‑4568643

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: ‑