Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	568,992	A→T	6.8%	E285V (GAA→GTA)	mcsB →	possible phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	568,992	0	A	T	6.8%	209.5 / 4.6	70	E285V (GAA→GTA)	mcsB	possible phosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (37/26);  new base T (3/2);  total (42/28)
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.73e-01

TTCATGAAGAAAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTAAAAGAAAAACGAGCAGATAT  >  CP000730/568854‑569127                                                                                                                                           |                                                                                                                                        ttcatGAAGAAAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAg                                                                                                                                       <  1:599341/141‑1 (MQ=255) ttcatGAAGAAAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAg                                                                                                                                       >  1:466505/1‑141 (MQ=255) ttcatGAAGAAAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAg                                                                                                                                       >  2:515503/1‑141 (MQ=255) ttcatGAAGAAAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAg                                                                                                                                       <  2:101712/141‑1 (MQ=255)  tcatGAAGAAAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGt                                                                                                                                      >  1:212918/1‑141 (MQ=255)   catGAAGAAAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGtt                                                                                                                                     >  1:638488/1‑141 (MQ=255)      gaagaaAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTGTCGCTCGCTAGGTATTTGACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTaaa                                                                                                                                  >  2:534739/1‑141 (MQ=255)               cAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGtata                                                                                                                         >  1:80644/1‑141 (MQ=255)                     cGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTa                                                                                                                   >  2:579657/1‑141 (MQ=255)                     cGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTa                                                                                                                   <  1:171397/141‑1 (MQ=255)                         aaaaGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAAtt                                                                                                               >  1:403742/1‑141 (MQ=255)                           aaGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTAc                                                                                                             <  1:485731/141‑1 (MQ=255)                               tAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTg                                                                                                         <  1:8154/141‑1 (MQ=255)                                aGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTTTTTCTTTTTTTATTTTACAAAACTGTATAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAAATACAttaa                                                                                                        >  2:240598/1‑139 (MQ=255)                                aGACACTGATAATCAACTAGAAATACACGACCGTGTTTTTCGCTCGCTAAGTATTTTACAAAACTGTAGAATGATAACTATAGAAGAGGCTTCTTATCGATTACGGGAAGTTAAACGAGGTATAGATTTAAATTACGTTGa                                                                                                        <  2:518849/141‑1 (MQ=255)                                  acacTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAAt                                                                                                      >  1:190970/1‑141 (MQ=255)                                   cacTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAAtt                                                                                                     >  1:145949/1‑141 (MQ=255)                                   cacTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAAtt                                                                                                     <  1:309389/141‑1 (MQ=255)                                       tataATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACaa                                                                                                 >  1:25780/1‑141 (MQ=255)                                          aaTCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTCGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAAc                                                                                              >  2:402722/1‑141 (MQ=255)                                           aTCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACt                                                                                             <  2:219787/141‑1 (MQ=255)                                            ttcaTTACAAACACAAGTCTGTGTTTGTCTCTCGCTAGGTATTTTACAAAACTGGAGAATGATAATTATGGAAGAGCCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACtt                                                                                            <  2:658222/138‑1 (MQ=255)                                            tCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACtt                                                                                            <  1:454958/141‑1 (MQ=255)                                             ccaTTAGACACACAACACCGTGTTTTTCGCTCGCCAGGTATTTTTCATAACTGTAGAATGATATCTATGGACGCGTCTTCTTATAGAACTAGCGAAGTTAAACATGGTATAGATTTAAATTACATTGAATTACAAAACttt                                                                                           >  2:254094/3‑141 (MQ=255)                                             cAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACttt                                                                                           <  1:671897/141‑1 (MQ=255)                                              aaTTAGAAACACAAGACCGTGTTTTTTTTTCGTTATGTATTTTACAAAACTGTAGAATGATAACTATGGAAGCGGCTTCTTCTAGCTTAAGCGCAGTTAAACTCGGTATAGATTTAAATTACATTGAATTACAAAACTTTa                                                                                          >  2:127152/1‑141 (MQ=255)                                               aTTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGTAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTaa                                                                                         >  2:99829/1‑141 (MQ=255)                                               aTTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGTAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTaa                                                                                         >  2:99789/1‑141 (MQ=255)                                                     aacgccaggccGTGTTTGTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACGTTGAATTACAAAACTTTAAATTTaa                                                                                   <  2:371980/132‑1 (MQ=255)                                                     aaCACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTaa                                                                                   >  1:396547/1‑141 (MQ=255)                                                         caAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGaa                                                                               >  2:312743/1‑141 (MQ=255)                                                         caAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGaa                                                                               >  2:260028/1‑141 (MQ=255)                                                           aGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAAtt                                                                             <  2:208010/141‑1 (MQ=255)                                                               cGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATg                                                                         <  1:260028/141‑1 (MQ=255)                                                                gtgtTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATgg                                                                        >  1:338791/1‑141 (MQ=255)                                                                   tttttCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAg                                                                     <  1:515503/141‑1 (MQ=255)                                                                      ttCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTCACACTTGGTATAGAATTCCAATCCGTTGTATTTCACTATTTTATGTTCAATGAATTGATGGTGGCtc                                                                  >  1:63702/1‑140 (MQ=255)                                                                        ttttttttttttATTTTAAAAAACCGTAGAATGTTAACTATGGTAGATGGTTCTTATAGATTAAGCGTTTTTAAACTTGGTATAGTTTTTTTTTTCATTTAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCtata                                                                >  2:465485/12‑141 (MQ=255)                                                                        cgctcgctAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCtata                                                                >  1:449737/1‑141 (MQ=255)                                                                        cgctcgctAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCtata                                                                >  1:405986/1‑141 (MQ=255)                                                                         gctcgctAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATAc                                                               >  2:513898/1‑141 (MQ=255)                                                                           tcattAGTTATTTTACAAAACTGTAGAAGGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAg                                                             >  2:600457/5‑141 (MQ=255)                                                                            cggtAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGt                                                            >  2:599353/4‑141 (MQ=255)                                                                              ctAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCa                                                          >  1:415523/1‑141 (MQ=255)                                                                               tAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCAc                                                         <  1:523991/141‑1 (MQ=255)                                                                                aGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCAcc                                                        <  1:312743/141‑1 (MQ=255)                                                                                   tATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCAtt                                                     <  1:97119/141‑1 (MQ=255)                                                                                    aacttACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCAttt                                                    <  2:611319/138‑1 (MQ=255)                                                                                      tttACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCAttttt                                                  >  2:587791/1‑141 (MQ=255)                                                                                           aaaagtattgAATGATAATTATGGAGGAGGCTTCTTATAGATCAAACAGAGTTACAATTGAAGTAACTTTCAATAACATTGAATTACAAAACTTTAAATATAATGAATTGATGGTAACTATCCAGTAACCATTTtcataaa                                             >  2:139549/10‑138 (MQ=255)                                                                                            aaaCTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGa                                            >  1:291248/1‑141 (MQ=255)                                                                                            aaaCTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGa                                            >  1:421973/1‑141 (MQ=255)                                                                                             aaCTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGTAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGAt                                           <  1:99829/141‑1 (MQ=255)                                                                                             aaCTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGTAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGAt                                           <  1:99789/141‑1 (MQ=255)                                                                                             aaCTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGAt                                           >  1:515159/1‑141 (MQ=255)                                                                                                    gAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGaagaag                                    >  1:605257/1‑141 (MQ=255)                                                                                                        gATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGa                                >  2:141310/1‑141 (MQ=255)                                                                                                        gATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGa                                >  2:141826/1‑141 (MQ=255)                                                                                                        gATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGa                                <  2:396547/141‑1 (MQ=255)                                                                                                         aTAACAATGGAAGATGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTAAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGAc                               <  2:552763/141‑1 (MQ=255)                                                                                                                   aaGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTaa                     >  2:100427/1‑141 (MQ=255)                                                                                                                    agagGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTaaa                    >  1:428360/1‑141 (MQ=255)                                                                                                                     gagGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTaaaa                   <  2:376141/141‑1 (MQ=255)                                                                                                                            cttATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAAt                                                                                       <  2:579450/66‑1 (MQ=255)                                                                                                                            cttATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAAt                                                                                       >  1:579450/1‑66 (MQ=255)                                                                                                                            cttATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTAAAAGAAAAAc            >  2:283982/1‑141 (MQ=255)                                                                                                                            cttATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTAAAAGAAAAAc            <  2:415523/141‑1 (MQ=255)                                                                                                                              tataGATTAAGCGAAGTTATACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTAAAAGAAAAACGa          <  2:315027/141‑1 (MQ=255)                                                                                                                                    ttAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTAAAAGAAAAACGAGCAGat    <  2:377477/141‑1 (MQ=255)                                                                                                                                      aaGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTAAAAGAAAAACGAGCAGatat  >  1:220064/1‑141 (MQ=255)                                                                                                                                           |                                                                                                                                        TTCATGAAGAAAAACAAATACGACAAAAGTTAGACACTTATAATCAATTAGAAACACAAGACCGTGTTTTTCGCTCGCTAGGTATTTTACAAAACTGTAGAATGATAACTATGGAAGAGGCTTCTTATAGATTAAGCGAAGTTAAACTTGGTATAGATTTAAATTACATTGAATTACAAAACTTTAAATTTAATGAATTGATGGTAGCTATACAGTCACCATTTTTATTAGATGAAGAAGATGACAAATCTGTAAAAGAAAAACGAGCAGATAT  >  CP000730/568854‑569127

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