Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	66,465	T→C	5.6%	R14R (AGA→AGG)	ccu ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	66,465	0	T	C	5.6%	235.7 / 3.7	90	R14R (AGA→AGG)	ccu	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/40);  new base C (2/3);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

CATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCCCTCCTTCTT  >  CP000730/66326‑66603                                                                                                                                            |                                                                                                                                           cATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTACTCTGTATTAAGTTGATTGTCTCGTTCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCctc                                                                                                                                           <  2:204660/141‑1 (MQ=255) cATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCctc                                                                                                                                           <  1:97215/141‑1 (MQ=255) cATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCctc                                                                                                                                           <  1:134714/141‑1 (MQ=255)    gtgtAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTtg                                                                                                                                        >  2:739833/1‑141 (MQ=255)      gtAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCCCTtgtg                                                                                                                                      >  2:118824/1‑141 (MQ=255)      gtAATTGACAACATATTATCATCTTCATGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCCCTtgtg                                                                                                                                      >  2:502734/1‑141 (MQ=255)      gtAATTGACAACATATTATCATCTTATTGTAGTTTGAAAATGTGTAGGTATTCTGTATTACGTTGATTGTCTCGATCTTTGACCATGTTCCTATGTAAGCTTTGACGGTCTGGAGCTAGTTGTTCACTCATGCccctactg                                                                                                                                      >  2:300843/1‑137 (MQ=255)           tGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTa                                                                                                                                 >  2:132333/1‑141 (MQ=255)           tGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAAACTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTa                                                                                                                                 >  1:555160/1‑141 (MQ=255)             acaacaTATTATCATCTTCTGGTAGTTTGAAAATGTGTAGGTAATCTGTCTTAAGTTGATTGTCTCGATCTTTTACCATGTTGCAGAGTAAGATTTAAAGGTCTGGAGATAGTTGTTAACTAATGGATCTtgtagtgtttc                                                                                                                               >  1:507087/1‑133 (MQ=37)                acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCgat                                                                                                                            <  2:557612/141‑1 (MQ=255)                acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCgat                                                                                                                            >  1:435105/1‑141 (MQ=255)                acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCgat                                                                                                                            >  2:196754/1‑141 (MQ=255)                acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTAGTGATATATCgat                                                                                                                            >  1:21458/1‑141 (MQ=255)                acaTATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCACGATCTTTTACCATGTTCCACAGTAAGATTTAAAGGTCTAGAGATCTGTCTTCCCTAATCCCTCTTGTGACGTAtggct                                                                                                                            >  2:313871/1‑137 (MQ=25)                       atcatcTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGAtttt                                                                                                                     >  2:718740/1‑141 (MQ=255)                          atcTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCAt                                                                                                                  <  1:95134/141‑1 (MQ=255)                                ttGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCAAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTAtt                                                                                                            >  2:619410/1‑141 (MQ=255)                                  gTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATGGAATGATTTTCATACTAtttt                                                                                                          >  2:125246/1‑141 (MQ=255)                                  gTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTAtttt                                                                                                          >  2:547661/1‑141 (MQ=255)                                  gTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACCAATGCCGCTTGTGATGTATCGATTGATTTTCATACTAtttt                                                                                                          >  2:479612/1‑141 (MQ=255)                                    aGTTTGAAAATGTTTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTATCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTTTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTtcc                                                                                                        <  2:226215/141‑1 (MQ=255)                                    aGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTtcc                                                                                                        >  2:606506/1‑141 (MQ=255)                                          aaaaTGTGTAGGCAATCTGGATTAAGTTGATTGTCTCGATCTTTTTCCATGTTCCAAAGTAATATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCAtt                                                                                                  <  2:455558/141‑1 (MQ=255)                                              tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATTTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  2:689063/141‑1 (MQ=255)                                              tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  2:22067/141‑1 (MQ=255)                                              tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  1:647794/141‑1 (MQ=255)                                              tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  1:749569/141‑1 (MQ=255)                                              tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              >  1:91211/1‑141 (MQ=255)                                              tgtgTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGc                                                                                              <  1:180006/141‑1 (MQ=255)                                               gtgtAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGCTTGATTTTAATACTATTTTCGTCCCTTTTGCt                                                                                             >  1:267010/1‑141 (MQ=255)                                               gtgtAGGTAATATGTATTAAGTTGATTGTCTCGATCTTTTAACATGTTCCAAAGTAAGATTTGAAGGTCTAGCGATAGTTGTTCACTAATGCCTCTTGTGATGTATGGATTGATTTACTTAATGTCTTCCTCCATTTTgat                                                                                             >  2:37737/1‑139 (MQ=37)                                                  tAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCtttt                                                                                          <  1:472388/141‑1 (MQ=255)                                                  tAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCtttt                                                                                          <  1:727097/141‑1 (MQ=255)                                                    ggTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTtttct                                                                                        >  2:291286/1‑141 (MQ=255)                                                         tCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCat                                                                                   <  2:82844/141‑1 (MQ=255)                                                          cTGTCTTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTTAGATTTGAATGTTTATTTAGAGTTTTTCACTAATGCCTCTTGGTATGTATTGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCatg                                                                                  <  2:369174/141‑1 (MQ=255)                                                           tGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCatga                                                                                 >  1:668213/1‑141 (MQ=255)                                                           tGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATGGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTTatga                                                                                 >  2:621178/1‑141 (MQ=255)                                                             tATTAAGTTGATTGTCCCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGGTGTATCGATGGGTTTTCTTACTATTTTCCTCAAGTTTGCTTTTCTTTCatgatg                                                                               >  1:303945/1‑141 (MQ=37)                                                               ttAAGTTGATTGTCTCGATCTTTTCCCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTc                                                                             <  1:492717/141‑1 (MQ=255)                                                               ttAAGTTGATTGTCTCGATCCTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTc                                                                             <  1:323847/141‑1 (MQ=255)                                                                 aaGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATCCTATTTTCCTCCATTTTGTTTCTCTTTCATGATGTCaa                                                                           <  1:452142/141‑1 (MQ=255)                                                                 aaGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCaa                                                                           <  1:629682/141‑1 (MQ=255)                                                                 aaGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCaa                                                                           <  1:251508/141‑1 (MQ=255)                                                                 aaGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCaa                                                                           <  1:312320/141‑1 (MQ=255)                                                                  aGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAAt                                                                          <  1:36629/141‑1 (MQ=255)                                                                        ttGTCTCGATCTTTTACCATGTTCCAAATTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTc                                                                    >  1:512617/1‑141 (MQ=255)                                                                        ttGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTc                                                                    >  1:344528/1‑141 (MQ=255)                                                                           tctcGATCTTCTACCCTGTTCCAAAGTAAGATTTGAAGGCCCAGAGATAGTTGTTCACTAATGCCTTTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGtt                                                                 <  2:3103/141‑1 (MQ=25)                                                                             tcGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTaa                                                               <  1:739833/141‑1 (MQ=255)                                                                              cGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAAt                                                              >  2:418898/1‑141 (MQ=255)                                                                                 tCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATgct                                                           >  2:32402/1‑139 (MQ=255)                                                                                 tCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATACCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCCCTTCGTTAATgcc                                                           >  2:319409/1‑139 (MQ=255)                                                                                     ttACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGCCTGTa                                                       >  2:314580/1‑141 (MQ=255)                                                                                      tatcATGTTCCAAAGTAAGATTTGAGGGTCTATAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTaa                                                      <  1:784975/138‑1 (MQ=255)                                                                                      tACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTCATGACTGTaa                                                      >  1:552731/1‑141 (MQ=255)                                                                                      tACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTaa                                                      <  1:287156/141‑1 (MQ=255)                                                                                      tACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTaa                                                      <  1:210896/141‑1 (MQ=255)                                                                                       tccATGTTCCAAAGTAAGCTTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAAc                                                     <  1:161249/140‑1 (MQ=255)                                                                                        ccATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACa                                                    <  2:707425/141‑1 (MQ=255)                                                                                            gTTCCAAAGTAGGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGatt                                                >  2:95137/1‑140 (MQ=255)                                                                                            gTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGatt                                                >  1:233273/1‑140 (MQ=255)                                                                                              tCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATAtt                                              >  1:466535/1‑141 (MQ=255)                                                                                               ccAAAGTAAGATTTTAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGacagtt                                             <  2:559517/141‑6 (MQ=255)                                                                                               ccAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATAttt                                             >  1:717923/1‑141 (MQ=255)                                                                                                ccaaGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTtg                                            <  2:377693/139‑1 (MQ=255)                                                                                                cAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTtg                                            <  2:377669/141‑1 (MQ=255)                                                                                                   atatagcttttAAGGTCTAGCGATAGTTGTTCACTAATGCCTCTTGTGGTTTTTCGATTGATTTCCATCCTTTTTTCATCCATTTTTCTTTTCTCTTATAATTTCAATAACTTAGTCTATGACTGGATCAGATATTTGTGc                                         <  2:273788/134‑1 (MQ=255)                                                                                                   aGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTATTTCATGATGTCAATCACTTCGTTAATGACTGTAACCGATATTTGTGc                                         >  1:421806/1‑141 (MQ=255)                                                                                                     tAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCa                                       <  2:36895/141‑1 (MQ=255)                                                                                                      aaGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATAACTTCGTTAATGACTGTAACAGATATTTTTGcccc                                      >  2:266008/1‑139 (MQ=255)                                                                                                        gATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTATTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACGTCGTTCATGACTGTAACAGATTTTTTGACCACtt                                    >  2:425040/1‑141 (MQ=25)                                                                                                            tGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCGATCACTTCGTTAATGACTGTAACAGATATTTGTGCAACTTTGAt                                >  1:150832/1‑141 (MQ=255)                                                                                                            tGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGAt                                >  1:762738/1‑141 (MQ=255)                                                                                                              aaGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATcc                              <  1:354420/141‑1 (MQ=255)                                                                                                                ggTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCaa                            <  2:3412/141‑1 (MQ=255)                                                                                                                 gtttAGCGATATGTGTTCACTAACGCCTTTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           <  2:555160/138‑1 (MQ=37)                                                                                                                 gTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           <  2:495259/141‑1 (MQ=255)                                                                                                                 gTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           <  1:132333/141‑1 (MQ=255)                                                                                                                 gTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           >  2:31069/1‑141 (MQ=255)                                                                                                                 gTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAAAGACTGTAACAGATATTTGTGCCACTTTGATCCAat                           >  2:360957/1‑141 (MQ=255)                                                                                                                   cTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGCCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAatta                         <  2:605857/141‑1 (MQ=255)                                                                                                                           aGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTAGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATgat                 >  2:464241/1‑139 (MQ=255)                                                                                                                           aGTTGTTCACAAATGCCGCTTGTGATGTATCGATGGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTCACAGATATTTGTGCCACTTTGATCCAATTATTCATGGt                 >  2:275144/1‑141 (MQ=255)                                                                                                                              tgttCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTCGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTAt              >  1:298654/1‑141 (MQ=255)                                                                                                                               gttCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTGTTTTGCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTAtt             >  2:608161/1‑141 (MQ=255)                                                                                                                                    cTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCCTGGTTATTCcttc        >  2:351518/1‑141 (MQ=255)                                                                                                                                     tAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCcctcc       <  2:346501/141‑1 (MQ=255)                                                                                                                                      aaTGCCCCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCcctcct      <  1:502734/141‑1 (MQ=255)                                                                                                                                      aaTGCCCCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCcctcct      <  1:118824/141‑1 (MQ=255)                                                                                                                                      aaTGCCCCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCcctcct      <  1:300843/141‑1 (MQ=255)                                                                                                                                         gCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTGATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCCCTCcttct   >  1:521960/1‑141 (MQ=255)                                                                                                                                          ccTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCCCTCcttctt  <  2:435105/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                           CATGTGTAATTGACAACATATTATCATCTTCTTGTAGTTTGAAAATGTGTAGGTAATCTGTATTAAGTTGATTGTCTCGATCTTTTACCATGTTCCAAAGTAAGATTTGAAGGTCTAGAGATAGTTGTTCACTAATGCCTCTTGTGATGTATCGATTGATTTTCATACTATTTTCCTCCATTTTGCTTTTCTTTCATGATGTCAATCACTTCGTTAATGACTGTAACAGATATTTGTGCCACTTTGATCCAATTATTCATGGTTATTCCCTCCTTCTT  >  CP000730/66326‑66603

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