Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,238,711	A→G	I360V (ATC→GTC)	yjbH →	YjbH family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,238,711	0	A	G	98.6%	273.2 / ‑4.6	73	I360V (ATC→GTC)	yjbH	YjbH family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (43/29);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTGGCGCAAAAACGCGTCG  >  NC_000913/4238574‑4238846                                                                                                                                          |                                                                                                                                        aTGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTc                                                                                                                                       >  2:389886/1‑140 (MQ=255) aTGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTc                                                                                                                                       >  1:1139078/1‑140 (MQ=255)       gACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACg                                                                                                                                 >  2:417514/1‑140 (MQ=255)        aCTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGa                                                                                                                                >  1:949972/1‑140 (MQ=255)          ttGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATc                                                                                                                              <  2:15391/140‑1 (MQ=255)            gCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCg                                                                                                                            >  1:289276/1‑140 (MQ=255)             ccGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGc                                                                                                                           <  2:286859/140‑1 (MQ=255)              cGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCa                                                                                                                          <  1:1073715/140‑1 (MQ=255)                  ccACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTAc                                                                                                                      <  2:730267/140‑1 (MQ=255)                        aTCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGaaaa                                                                                                                <  1:1030468/140‑1 (MQ=255)                             ggCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGaaa                                                                                                                 <  2:10613/134‑1 (MQ=255)                             ggCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGaaa                                                                                                                 >  1:10613/1‑134 (MQ=255)                             ggCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAAt                                                                                                               <  2:400965/136‑1 (MQ=255)                             ggCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAAt                                                                                                               >  1:400965/1‑136 (MQ=255)                             ggCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGcc                                                                                                           >  1:266550/1‑140 (MQ=255)                                 aaaGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTaa                                                                                                       <  2:266550/140‑1 (MQ=255)                                   aGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACa                                                                                                     >  2:1261098/1‑140 (MQ=255)                                   aGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACa                                                                                                     >  2:1049561/1‑140 (MQ=255)                                    ggCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATAGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACAt                                                                                                    >  1:1121965/1‑140 (MQ=255)                                        aTACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCg                                                                                                >  2:651592/1‑140 (MQ=255)                                               gTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGc                                                                                                  >  2:275217/1‑131 (MQ=255)                                               gTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGc                                                                                                  <  1:275217/131‑1 (MQ=255)                                               gTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTga                                                                                         >  2:364478/1‑140 (MQ=255)                                                tATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTgac                                                                                        >  2:245810/1‑140 (MQ=255)                                                tATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCACATTACGGAAAATCGCCTTAACATGCCGCAGGTgac                                                                                        >  1:775098/1‑140 (MQ=255)                                                  tGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTgacga                                                                                      <  2:949972/140‑1 (MQ=255)                                                      aCCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGa                                                                                  >  1:706527/1‑140 (MQ=255)                                                       ccGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGaa                                                                                 <  1:364478/140‑1 (MQ=255)                                                        cGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGaaa                                                                                <  2:1139078/140‑1 (MQ=255)                                                        cGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGaaa                                                                                <  1:389886/140‑1 (MQ=255)                                                        cGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGaaa                                                                                >  2:232336/1‑140 (MQ=255)                                                          gCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCg                                                                                                <  1:52618/122‑1 (MQ=255)                                                          gCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCg                                                                                                >  2:52618/1‑122 (MQ=255)                                                          gCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAAcc                                                                              >  2:195754/1‑140 (MQ=255)                                                          gCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAAcc                                                                              >  1:427838/1‑140 (MQ=255)                                                          gCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAAcc                                                                              >  1:869493/1‑140 (MQ=255)                                                            gAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGa                                                                            >  1:962380/1‑140 (MQ=255)                                                            gAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGa                                                                            >  2:6735/1‑140 (MQ=255)                                                                  gTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGc                                                                      <  1:231761/140‑1 (MQ=255)                                                                   tGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGcc                                                                     <  1:506119/140‑1 (MQ=255)                                                                      aaTATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGc                                                                  <  1:1212834/140‑1 (MQ=255)                                                                            gTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGaaa                                                            >  1:1053525/1‑140 (MQ=255)                                                                                  cgcgcgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAatcat                                                      <  2:289276/140‑1 (MQ=255)                                                                                  cgcgcgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAatcat                                                      <  1:1266628/140‑1 (MQ=255)                                                                                   gcgcgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAatcatc                                                     >  1:1103369/1‑140 (MQ=255)                                                                                    cgcgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCAtct                                                    >  2:48576/1‑140 (MQ=255)                                                                                      cgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCg                                                  <  2:962380/140‑1 (MQ=255)                                                                                       gAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGc                                                 <  1:1125182/140‑1 (MQ=255)                                                                                            gATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGgag                                            <  1:1261098/140‑1 (MQ=255)                                                                                               catcCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGc                                         >  1:446801/1‑140 (MQ=255)                                                                                                 tcCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCg                                       >  1:265606/1‑140 (MQ=255)                                                                                                  ccGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGt                                      <  1:1049561/140‑1 (MQ=255)                                                                                                     tGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTgg                                   >  2:141870/1‑140 (MQ=255)                                                                                                     tGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTgg                                   >  2:1140049/1‑140 (MQ=255)                                                                                                     tGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTgg                                   >  1:1000379/1‑140 (MQ=255)                                                                                                      gCGAATCGGATCGTGATGACCGATCTGCCGGATGGGGTCAAAACGATCCACATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTggg                                  <  2:775098/140‑1 (MQ=255)                                                                                                       cGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGc                                                 <  2:827388/124‑1 (MQ=255)                                                                                                       cGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGc                                                 >  1:827388/1‑124 (MQ=255)                                                                                                       cGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGc                                 >  2:1179902/1‑140 (MQ=255)                                                                                                       cGAATCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGc                                 >  2:1108281/1‑140 (MQ=255)                                                                                                         aaTCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCa                               >  2:344382/1‑140 (MQ=255)                                                                                                         aaTCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCa                               <  2:1053525/140‑1 (MQ=255)                                                                                                          aTCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCAc                              <  1:1062264/140‑1 (MQ=255)                                                                                                           tCGGATCGTGATGAACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACg                             >  1:578711/1‑140 (MQ=255)                                                                                                                    gatgaACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGc                    >  1:229394/1‑140 (MQ=255)                                                                                                                    gatgaACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGc                    <  2:706527/140‑1 (MQ=255)                                                                                                                      tgaACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTg                  >  2:770214/1‑140 (MQ=255)                                                                                                                       gaACGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTgg                 <  2:427838/140‑1 (MQ=255)                                                                                                                         aCGATCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTGgcg               >  1:152643/1‑140 (MQ=255)                                                                                                                            aTCTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTGGCGCaa            >  2:597920/1‑140 (MQ=255)                                                                                                                              cTGCCGGATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTGGCGCaaaa          >  1:658932/1‑140 (MQ=255)                                                                                                                                   ggATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTGGCGCAAAAAcgcg     <  2:1103369/140‑1 (MQ=255)                                                                                                                                    gATGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTGGCGCAAAAACGCGt    >  2:589198/1‑140 (MQ=255)                                                                                                                                      tGGGGTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTGGCGCAAAAACGCGTCg  >  2:707053/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        ATGCCGGACTTGCCGATCCACAGATCCAGGCGAAAGGCGATACGCTGTATGTTACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGACGCTGGCGCAAAAACGCGTCG  >  NC_000913/4238574‑4238846

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