Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	419,825	A→C	17.7%	L51R (CTT→CGT)	U712_01910 ←	Uncharacterized protein yczF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	419,825	0	A	C	17.7%	149.8 / 10.1	73	L51R (CTT→CGT)	U712_01910	Uncharacterized protein yczF
Reads supporting (aligned to +/- strand):  ref base A (44/14);  new base C (2/12);  total (46/27)
Fisher's exact test for biased strand distribution p-value = 3.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATCAACAAACACAGCAGTGATACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCGGTCACG  >  CP006881/419689‑419962                                                                                                                                          |                                                                                                                                          aCCAACAACAAACACAGCCGTTAAAAAATAGACGGAAATATGATTGTCTGAAAAACGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCACACAAGTACTAACAAAAAAAATaaaacc                                                                                                                                         <  6:128675/139‑3 (MQ=255)                agcagTGATACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTCTTTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCg                                                                                                                            >  1:266579/1‑138 (MQ=255)                agcagTGATACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGc                                                                                                                           >  3:183325/1‑138 (MQ=255)                agcagTGATACAAAAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCATTAGGGGTTC‑GTTTTTTTCGATTTATAGATCACGGACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGc                                                                                                                           >  3:264953/1‑138 (MQ=255)                      gATACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCAGCAAAtgtt                                                                                                                    >  1:154658/1‑139 (MQ=255)                       aTACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACATGTACTAACACAAAAAATAAACGCACGCCCCATTCCGCAAAtgttg                                                                                                                   >  2:58039/1‑139 (MQ=255)                       aTACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACCCAAAAAATAAAAGCACGACCCATACCGCTAAtgttg                                                                                                                   >  8:296652/1‑139 (MQ=255)                       aTACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAAtgttg                                                                                                                   >  4:45327/1‑139 (MQ=255)                       aTACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAAtgttg                                                                                                                   >  7:117387/1‑139 (MQ=255)                       aTACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAAtgttg                                                                                                                   >  3:181492/1‑139 (MQ=255)                        tACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAAtgttgt                                                                                                                  >  3:36010/1‑139 (MQ=255)                             tAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAAtgttgt                                                                                                                  >  8:57112/1‑134 (MQ=255)                             tAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAAtgttgt                                                                                                                  <  7:57112/134‑1 (MQ=255)                                  ggCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGCCCCATTCCGCAAATGTTGTCATCTTAAAg                                                                                                        >  3:56791/1‑139 (MQ=255)                                  ggCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTaacg                                                                                                        >  1:255680/1‑137 (MQ=255)                                  ggCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCAGCAAATGTTGTCATCTTAAAg                                                                                                        >  5:239368/1‑139 (MQ=255)                                   caaaTATCATTGTCTCAAAAACGCCCCCCCCTTTAGGGGTTC‑TTTTTTTTCTTTTTATAGATCACGTACATTTCTCACACAAGTAAAAACACAAAAAAAAAACCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGc                                                                                                       <  1:102917/137‑1 (MQ=255)                                     aaTATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTaactgt                                                                                                      >  7:84432/1‑134 (MQ=255)                                     aaTATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTg                                                                                                     >  3:268550/1‑139 (MQ=255)                                      atatGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTg                                                                                                     >  5:136955/1‑138 (MQ=255)                                         tGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGCCCCATTCCGCAAATGTTGTCATCTTAAAGCTGGAc                                                                                                  >  6:50390/1‑138 (MQ=255)                                            ttGACTTAAAAACGCCCCCCCTTTTTGGGGTTT‑TTTTTTTTCGATTTATAGATCACTTAATTTTCACAAAAAACAACTAAAAAAAAAAAAAACCCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACata                                                                                               <  5:331967/138‑1 (MQ=255)                                              gTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAACGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTg                                                                                            >  1:208492/1‑139 (MQ=255)                                              gTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTg                                                                                            >  3:30551/1‑139 (MQ=255)                                              gTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTg                                                                                            >  3:83969/1‑139 (MQ=255)                                               tCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGt                                                                                           >  1:305107/1‑139 (MQ=255)                                               tCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGt                                                                                           >  7:283771/1‑139 (MQ=255)                                                 tGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATAttctgt                                                                                         >  3:123884/1‑135 (MQ=255)                                                     aaacGCCCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAAAAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATaaa                                                                                     <  8:255958/139‑1 (MQ=255)                                                         gcgcCCCAATTTTGGGGTTT‑TTTTTTTTTTATTTAAAGAACACGTACATTACTCCCAAAATTACTAACACAAAAAAAAAAACCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGc                                                                                 <  8:328680/139‑1 (MQ=255)                                                           gcCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACACATTCCGCAAATGTTGTCATCTTAATGCTGGACATATTGTTATAAGCGGCtt                                                                               >  1:175711/1‑139 (MQ=255)                                                                    tAGGGGTTC‑TTTTTTTTCGATTTATAGATAACGTTCATTTCTCTCACAATTTCTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATTTTTTCTTCTTAAATCTGGACATATTGTTATAAACGGCTTTTTTTAAAc                                                                      >  1:331952/1‑139 (MQ=255)                                                                      ggtgTTC‑TTTTTTTTCTTTTTAAATACCACTTAAATTTCTCACAAAATTAATAACAAAAAAAAAAAAACCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTg                                                                    <  3:45327/136‑1 (MQ=255)                                                                      ggggTTA‑TTTTTTTTCGATTTATAGACCACGTACATTTCTCCCACAAATACAAACACAAAAAATAACCGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTg                                                                    <  6:177906/139‑1 (MQ=255)                                                                        tgTTT‑TTTTTTTTCTATTTAAAGATCACTTACATAACTAAAAAATCTAAAAACAAAAAAAAAAAACCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGaa                                                                  <  4:183295/138‑1 (MQ=255)                                                                        gtttC‑TTTTTTTTCTTTTTATAGATAAATTACATTTCTCTCACAATAAATAAAAAAAAAAAAAAAACCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGaa                                                                  <  8:214197/137‑1 (MQ=255)                                                                        gtttC‑TTTTTTTTCGTTTTATAGACCACTTACTTTACTCTCACAATTAAAAAAAAAAAAAATAACCCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGaa                                                                  <  1:98543/137‑1 (MQ=255)                                                                        ggTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCACACAACTACAAAAAAAAAAAATAAACGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGaa                                                                  <  4:32382/139‑1 (MQ=255)                                                                         gTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGCCCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAAc                                                                 >  6:91489/1‑139 (MQ=255)                                                                          ttC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAAcc                                                                >  1:196681/1‑139 (MQ=255)                                                                           tC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAAccc                                                               >  3:119911/1‑139 (MQ=255)                                                                            c‑TTTTTTTTCTATTTATAGATCACGTACTTTTCTCTCACAAGAACTAAAAAAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAAccc                                                               <  4:268550/138‑1 (MQ=255)                                                                            c‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCt                                                              >  4:142769/1‑139 (MQ=255)                                                                            c‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCt                                                              >  7:145479/1‑139 (MQ=255)                                                                            c‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCt                                                              >  3:257070/1‑139 (MQ=255)                                                                            c‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCt                                                              >  5:121852/1‑139 (MQ=255)                                                                            c‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCt                                                              >  7:268692/1‑139 (MQ=255)                                                                            c‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCt                                                              >  5:40557/1‑139 (MQ=255)                                                                            c‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGATTTTGTTAAACTGAACACt                                                              >  1:170889/1‑139 (MQ=255)                                                                              ttttttttCGTTTTATCTATAATTTATATTTTTCAAAAAAGTTGAAACAAAAAAATTAATTGCACGACCGATTCCGTAAATGTTGTCATCTTAATGCTTGACATATTGTTATAAACGGCTTTTGTTAAACTGAAcgat                                                              <  2:175711/138‑4 (MQ=255)                                                                                ttttttCTATTTATACATCACGAACATTTCTCACACAACTACAAACACAAAAAATAAACGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCt                                                              <  6:150766/136‑1 (MQ=255)                                                                                            atccaTCACGTAAATTTCCCTCACAACAACAAACACAAAAAAAAACCCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATc                                               <  2:208492/135‑1 (MQ=255)                                                                                             tataTCACTTACATTTCTCACACAATAAAAAAAAAAAAAAAAAAAACCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATc                                               <  4:183325/135‑1 (MQ=255)                                                                                                aTCACTTAAATTTCTCTCACAAGAACAAAAAAAAAAAATAAACCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAaca                                           <  4:36010/139‑1 (MQ=255)                                                                                                   aCGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAaca                                           >  5:151725/1‑136 (MQ=255)                                                                                                   aCGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAaca                                           >  3:210296/1‑136 (MQ=255)                                                                                                       aCATTACACACACAACAACAAACAAAAAAAAAAACCCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAAt                                    <  2:332909/139‑1 (MQ=255)                                                                                                               ctcACAAGTACTAACACAAAAAATAAAAGCACGCCCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTCTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCtaat                             >  7:7268/1‑138 (MQ=255)                                                                                                                 cacaAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATaaa                          <  6:121852/139‑1 (MQ=255)                                                                                                                     aGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACaaa                      >  8:3462/1‑139 (MQ=255)                                                                                                                      ttaatAACACAAAAAATAAAACCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAAt                     <  2:24273/135‑1 (MQ=255)                                                                                                                       tACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTTAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAAt                     >  4:7138/1‑138 (MQ=255)                                                                                                                            acacAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAat                 <  8:36588/137‑1 (MQ=255)                                                                                                                            acacAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAat                 >  7:36588/1‑137 (MQ=255)                                                                                                                            acacAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGACCATCAACACAGAAATAGCTAATAAACAAATTAatat               >  7:298885/1‑139 (MQ=255)                                                                                                                            acacAAAAAATAAAAGCACGACCCAAACCGAAAATGTTGTCATATAAAAGCAGAACATATTGTTATAAACGCCTTTTCAAAAACTGAACCCTTGAAGCATGTACATCAACACAGAAATAACTAATAAACAAAATAatat               >  3:26200/1‑139 (MQ=255)                                                                                                                                aaaaaaaaaaCCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAAt           <  8:145479/139‑1 (MQ=255)                                                                                                                                aaaaaaTAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAAt           <  4:181492/139‑1 (MQ=255)                                                                                                                                aaaaaaTAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAAt           <  8:117387/139‑1 (MQ=255)                                                                                                                                 aaaaaTAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATc          >  2:46191/1‑139 (MQ=255)                                                                                                                                  aaaaTAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATcc         >  4:270746/1‑139 (MQ=255)                                                                                                                                     aTAAACCCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCGGt      <  3:227356/139‑1 (MQ=255)                                                                                                                                     aTAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCggg      <  3:68661/139‑2 (MQ=255)                                                                                                                                     aTAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCGGt      <  3:270746/139‑1 (MQ=255)                                                                                                                                       aaacGCACGACCCATTCAGCAAATGTAGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTTAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCGGTCa    <  7:84392/139‑1 (MQ=255)                                                                                                                                       aaaaGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCGGTCa    <  7:3462/139‑1 (MQ=255)                                                                                                                                         aaGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCGGTCACg  <  3:142769/139‑1 (MQ=255)                                                                                                                                          aGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCGGTCACg  >  1:129129/1‑138 (MQ=255)                                                                                                                                          |                                                                                                                                          ACCATCAACAAACACAGCAGTGATACAATAGACGGCAATATGATTGTCTGAAAAAAGCGCCCCCCTTTAGGGGTTC‑TTTTTTTTCGATTTATAGATCACGTACATTTCTCTCACAAGTACTAACACAAAAAATAAAAGCACGACCCATTCCGCAAATGTTGTCATCTTAAAGCTGGACATATTGTTATAAACGGCTTTTGTTAAACTGAACCCTTGAAGCATGTCCATCAACACAGAAATAGCTAATAAACAAATTAATATGAATCCGGTCACG  >  CP006881/419689‑419962

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑