Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,401,459	0	T	C	23.5%	172.6 / 19.8	81	D625G (GAT→GGT)	csrD	regulator of CsrB and CsrC decay
Reads supporting (aligned to +/- strand):  ref base T (48/14);  new base C (0/19);  total (48/33)
Fisher's exact test for biased strand distribution p-value = 5.41e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTT  >  NC_000913/3401324‑3401596                                                                                                                                        |                                                                                                                                          aGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcccc                                                                                                                                       >  1:893074/1‑140 (MQ=255)  gCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCt                                                                                                                                      <  2:1302215/140‑1 (MQ=255)     gACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGc                                                                                                                                   <  1:915310/140‑1 (MQ=255)       cAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTGGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCATTGccc                                                                                                                                 <  1:348573/140‑1 (MQ=255)       cAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGccc                                                                                                                                 <  1:154374/140‑1 (MQ=255)       cAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGccc                                                                                                                                 <  1:1177929/140‑1 (MQ=255)            acacTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTg                                                                                                                            >  1:714634/1‑140 (MQ=255)                tCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAac                                                                                                                        >  1:1317968/1‑140 (MQ=255)                tCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAac                                                                                                                        >  1:1593325/1‑140 (MQ=255)                   cATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACAcc                                                                                                                     <  1:572374/140‑1 (MQ=255)                      gagaCGCAGCGCGCTTTATTCTACGTGAAAACGTATTAAACGGCAGGTTAATCCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACcgcg                                                                                                                  <  2:1345808/140‑1 (MQ=255)                      gagaCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACcgcg                                                                                                                  <  2:28125/140‑1 (MQ=255)                      gagaCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACcgcg                                                                                                                  <  1:1083042/140‑1 (MQ=255)                         aCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTg                                                                                                               >  1:1345808/1‑140 (MQ=255)                         aCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTAAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTg                                                                                                               <  2:858468/140‑1 (MQ=255)                            cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                            >  2:1408022/1‑140 (MQ=255)                              gcgcgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTTAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcgcg                                                                                                                  <  1:965348/132‑1 (MQ=255)                              gcgcgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTTAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcgcg                                                                                                                  >  2:965348/1‑132 (MQ=255)                              gcgcgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCa                                                                                                          >  2:581075/1‑140 (MQ=255)                                gcgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAgg                                                                                                        >  1:616636/1‑140 (MQ=255)                                gcgcATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAgg                                                                                                        >  1:786737/1‑140 (MQ=255)                                           tACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTc                                                                                             >  2:538482/1‑140 (MQ=255)                                              gTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCt                                                                                          >  2:1330716/1‑140 (MQ=255)                                              gTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCt                                                                                          >  1:1071335/1‑140 (MQ=255)                                               tGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCtt                                                                                         <  1:1556775/140‑1 (MQ=255)                                                gAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTc                                                                                        >  2:82574/1‑140 (MQ=255)                                                gAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTc                                                                                        >  1:561377/1‑140 (MQ=255)                                                  aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCa                                                                                                >  1:1494099/1‑130 (MQ=255)                                                  aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCa                                                                                                <  2:1494099/130‑1 (MQ=255)                                                           aaaCGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAacgcacg                                                                             >  1:1033016/1‑140 (MQ=255)                                                               ggCAGTTTAAACCGAGTATCTTTGTAAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCgg                                                                         <  2:461268/140‑1 (MQ=255)                                                                  aGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTgg                                                                      >  1:936690/1‑140 (MQ=255)                                                                   ggTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGc                                                                     <  1:515962/140‑1 (MQ=255)                                                                    gTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCg                                                                    >  1:1122435/1‑140 (MQ=255)                                                                      aaaaCCGGTTATTTTTGTAAATTTTTTTTCACTTTATTATCAGGGGGTTGTGGGGCCGCGAAAAAACCCCTTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTa                                                                  <  2:786737/139‑1 (MQ=255)                                                                          ccGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAAcc                                                              >  2:945293/1‑140 (MQ=255)                                                                              gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                          >  1:839228/1‑140 (MQ=255)                                                                              gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                          >  1:871145/1‑140 (MQ=255)                                                                              gTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTggg                                                          >  2:1344730/1‑140 (MQ=255)                                                                               tATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGt                                                         >  2:947199/1‑140 (MQ=255)                                                                               tATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGt                                                         >  1:458882/1‑140 (MQ=255)                                                                                     tgtgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGt                                                   >  1:649958/1‑140 (MQ=255)                                                                                     tgtgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGt                                                   >  1:551932/1‑140 (MQ=255)                                                                                       tgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                                 >  1:362941/1‑140 (MQ=255)                                                                                        aaaTTTTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTccc                                                <  2:893074/139‑1 (MQ=255)                                                                                            atTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACACGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg                                            >  2:1002458/1‑140 (MQ=255)                                                                                            atTTTTTCACGTTAGTATCAAGGGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAg                                            <  1:945293/140‑1 (MQ=255)                                                                                                tttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCa                                          >  2:1320689/1‑138 (MQ=255)                                                                                                tttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCa                                          <  1:1320689/138‑1 (MQ=255)                                                                                                 ttCACTTTAGTACAAAGTGGTTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGc                                       <  2:267227/140‑1 (MQ=255)                                                                                                 ttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGc                                       >  1:1074018/1‑140 (MQ=255)                                                                                                 ttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGc                                       <  1:1344730/140‑1 (MQ=255)                                                                                                 tccACGTTAGTATCAAGTGGCTGGGAGGACGGGAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGc                                       <  2:1033016/138‑1 (MQ=255)                                                                                                  tCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCt                                      <  1:82574/140‑1 (MQ=255)                                                                                                   cACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAACCCTGGGTGCTGGTCCCGGAGCAGGCtt                                     >  1:1203097/1‑140 (MQ=255)                                                                                                       ttAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccac                                 <  1:1330716/140‑1 (MQ=255)                                                                                                        tAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccacc                                <  2:1122435/140‑1 (MQ=255)                                                                                                                aGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                        >  1:775784/1‑140 (MQ=255)                                                                                                                aGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                        >  2:208952/1‑140 (MQ=255)                                                                                                                 gTGGCTGTGAGGACGCAAAAAAACCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGa                       <  2:1388842/140‑1 (MQ=255)                                                                                                                  tGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGaa                      >  2:865602/1‑140 (MQ=255)                                                                                                                    gCTTTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                    <  1:1431975/140‑1 (MQ=255)                                                                                                                    gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                    >  1:221491/1‑140 (MQ=255)                                                                                                                    gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                    >  2:389002/1‑140 (MQ=255)                                                                                                                    gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                    <  1:947199/140‑1 (MQ=255)                                                                                                                    gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                    >  2:115464/1‑140 (MQ=255)                                                                                                                    gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc                    >  1:1388842/1‑140 (MQ=255)                                                                                                                    gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAc                     >  1:334282/1‑139 (MQ=255)                                                                                                                    gCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAc                     <  2:334282/139‑1 (MQ=255)                                                                                                                     cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca                   >  1:990110/1‑140 (MQ=255)                                                                                                                     cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca                   >  1:256632/1‑140 (MQ=255)                                                                                                                     cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca                   >  2:1205880/1‑140 (MQ=255)                                                                                                                       gtgAGGACGCGAAAAAATCCCCTTGTCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagc                 >  1:1001452/1‑140 (MQ=255)                                                                                                                       gtgAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagc                 >  1:372439/1‑140 (MQ=255)                                                                                                                       gggAGGCCGCAAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagc                 <  2:616636/138‑1 (MQ=255)                                                                                                                               gcgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                        >  2:780421/1‑125 (MQ=255)                                                                                                                               gcgAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTg                        <  1:780421/125‑1 (MQ=255)                                                                                                                                cgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtc        >  2:1147197/1‑140 (MQ=255)                                                                                                                                cgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtc        >  1:1370338/1‑140 (MQ=255)                                                                                                                                     aaaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGt   <  2:551932/140‑1 (MQ=255)                                                                                                                                      aaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGtt  <  1:208952/140‑1 (MQ=255)                                                                                                                                      aaTCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGtt  <  1:1205880/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                          AGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTCTCCGTT  >  NC_000913/3401324‑3401596

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