Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,497,179	G→A	L59L (TTG→TTA)	nirC →	nitrite transporter NirC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,497,179	0	G	A	97.0%	220.1 / ‑2.0	66	L59L (TTG→TTA)	nirC	nitrite transporter NirC
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (28/36);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGC  >  NC_000913/3497044‑3497310                                                                                                                                        |                                                                                                                                    gcgcATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGa                                                                                                                                 >  1:521696/1‑140 (MQ=255)      ttGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGgc                                                                                                                            <  1:118679/140‑1 (MQ=255)      ttGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGgc                                                                                                                            <  1:375845/140‑1 (MQ=255)       tGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGgcg                                                                                                                           <  2:544462/140‑1 (MQ=255)        gCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGgcgc                                                                                                                          <  2:637594/140‑1 (MQ=255)         cACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGcgcg                                                                                                                         <  1:347353/140‑1 (MQ=255)               tGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCttt                                                                                                                   >  1:365839/1‑140 (MQ=255)               tGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGACGCGACCttt                                                                                                                   >  2:599233/1‑140 (MQ=255)                    gCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTAt                                                                                                              <  1:399755/140‑1 (MQ=255)                        aTAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTAt                                                                                                              >  1:136184/1‑136 (MQ=255)                        aTAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTAt                                                                                                              <  2:136184/136‑1 (MQ=255)                        aTAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGcc                                                                                                          >  1:905900/1‑140 (MQ=255)                          aaCCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCtt                                                                                                        <  1:729465/140‑1 (MQ=255)                              cGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACg                                                                                                    <  2:279040/140‑1 (MQ=255)                              cGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACg                                                                                                    <  1:350781/140‑1 (MQ=255)                              cGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACg                                                                                                    <  2:905900/140‑1 (MQ=255)                                            gTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGc                                                                                      <  2:47934/140‑1 (MQ=255)                                            gTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGc                                                                                      <  1:784375/140‑1 (MQ=255)                                            gTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGc                                                                                      <  1:599233/140‑1 (MQ=255)                                             tCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGcc                                                                                     >  1:527529/1‑140 (MQ=255)                                              cAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCg                                                                                    <  2:521696/140‑1 (MQ=255)                                                         tGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTg                                                                         <  1:703010/140‑1 (MQ=255)                                                           gcgggcgCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGa                                                                             <  1:362657/134‑1 (MQ=255)                                                           gcgggcgCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGa                                                                             >  2:362657/1‑134 (MQ=255)                                                           gcgggcgCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGtt                                                                       >  1:797616/1‑140 (MQ=255)                                                           gcgggcgCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGtt                                                                       >  2:945905/1‑140 (MQ=255)                                                             gggcgCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCa                                                                     <  2:909465/140‑1 (MQ=255)                                                              ggcgCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCAc                                                                    <  2:327099/140‑1 (MQ=255)                                                              ggcgCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCAc                                                                    <  1:292579/140‑1 (MQ=255)                                                                   gTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGac                                                               <  2:527529/140‑1 (MQ=255)                                                                    tATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGaca                                                              <  1:393079/140‑1 (MQ=255)                                                                              ttGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGacaca                                                            <  2:198206/132‑1 (MQ=255)                                                                              ttGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGacaca                                                            >  1:198206/1‑132 (MQ=255)                                                                              ttGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTc                                                    >  1:690040/1‑140 (MQ=255)                                                                              ttGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTc                                                    >  1:618317/1‑140 (MQ=255)                                                                               tGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTcc                                                   >  1:214764/1‑140 (MQ=255)                                                                                gggATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCt                                                  >  1:665081/1‑140 (MQ=255)                                                                                gggATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCt                                                  <  1:887180/140‑1 (MQ=255)                                                                                gggATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCt                                                  <  1:794167/140‑1 (MQ=255)                                                                                gggATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCt                                                  <  2:797616/140‑1 (MQ=255)                                                                                 ggATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTc                                                 >  2:711812/1‑140 (MQ=255)                                                                                    tcatcCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGacaca                                                            <  2:376174/126‑1 (MQ=255)                                                                                    tcatcCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGacaca                                                            >  1:376174/1‑126 (MQ=255)                                                                                     catcCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCt                                             >  2:818702/1‑140 (MQ=255)                                                                                     catcCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCt                                             >  1:241546/1‑140 (MQ=255)                                                                                      atcCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCtt                                            <  1:740442/140‑1 (MQ=255)                                                                                        ccTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTg                                          >  1:486015/1‑140 (MQ=255)                                                                                           gATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTgggg                                       <  2:486015/140‑1 (MQ=255)                                                                                                    gCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCgg                              <  2:241546/140‑1 (MQ=255)                                                                                                      tCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGgc                            <  2:346882/140‑1 (MQ=255)                                                                                                      tCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGgc                            <  2:214764/140‑1 (MQ=255)                                                                                                      tCGGTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGgc                            <  1:945905/140‑1 (MQ=255)                                                                                                        ggTAATTTGCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGgcag                          >  1:928897/1‑140 (MQ=255)                                                                                                               tGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCggg                             <  1:238740/130‑1 (MQ=255)                                                                                                               tGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCggg                             >  2:238740/1‑130 (MQ=255)                                                                                                                gCTCGATCCATCCGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGATCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGcc                  <  2:665081/140‑1 (MQ=255)                                                                                                                     atccatccGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAATGCGGGCAGCATCAGCCACggg             <  2:690040/140‑1 (MQ=255)                                                                                                                           ccGTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGTGCAGCATCAGCCACGGGCAAAtg       >  1:485158/1‑140 (MQ=255)                                                                                                                             gTACGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAAtgtg     <  2:928897/140‑1 (MQ=255)                                                                                                                                cGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGc  >  2:312187/1‑140 (MQ=255)                                                                                                                                cGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGc  >  1:241854/1‑140 (MQ=255)                                                                                                                                cGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGc  >  2:654898/1‑140 (MQ=255)                                                                                                                                cGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGc  >  1:847662/1‑140 (MQ=255)                                                                                                                                cGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGc  >  2:150574/1‑140 (MQ=255)                                                                                                                                cGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGc  >  2:421093/1‑140 (MQ=255)                                                                                                                                cGCCCTTTAGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGc  >  1:63455/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                    GCGCATTGCACGCCTGTCGGCAAATAACCCGCTCGGCTTTTGGGTCAGCTCCGCCATGGCGGGCGCGTATGTGGGTCTTGGGATCATCCTGATTTTCACGCTCGGTAATTTGCTCGATCCATCCGTACGCCCTTTGGTGATGGGCGCGACCTTTGGTATCGCCTTAACGCTGGTGATTATCGCCGGTTCTGAACTGTTCACCGGACACACCATGTTCCTCACCTTTGGGGTAAAAGCGGGCAGCATCAGCCACGGGCAAATGTGGGC  >  NC_000913/3497044‑3497310

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