Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,714,061	A→G	5.6%	C243R (TGT→CGT)	accC ←	biotin carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,714,061	0	A	G	5.6%	210.7 / 3.3	71	C243R (TGT→CGT)	accC	biotin carboxylase
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base G (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

AGCATTCATTTCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGTATAAAAGCCTCAAC  >  CP000730/1713921‑1714196                                                                                                                                             |                                                                                                                                        aGCATTCATTTCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGcaca                                                                                                                                         >  1:625038/1‑141 (MQ=255) aGCATTCATTTCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGcaca                                                                                                                                         >  1:535732/1‑141 (MQ=255)    attcattTCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACGAgg                                                                                                                                      >  2:117092/1‑141 (MQ=255)    attcattTCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACGAgg                                                                                                                                      >  2:553096/1‑141 (MQ=255)    attcattTCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAgg                                                                                                                                      <  1:501827/141‑1 (MQ=255)      tcattTCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGgag                                                                                                                                    <  1:76092/141‑1 (MQ=255)          ttCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGAttc                                                                                                                                >  1:121126/1‑141 (MQ=255)             aaTAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGAttcttc                                                                                                                             <  2:78240/141‑1 (MQ=255)                 aaGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTGATGCTGCACAAGGAGATTCTTCGAtt                                                                                                                         <  2:321221/141‑1 (MQ=255)                 aaGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGAtt                                                                                                                         >  1:231932/1‑141 (MQ=255)                  aGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTa                                                                                                                        >  2:627257/1‑141 (MQ=255)                     aaTGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAAtt                                                                                                                     >  1:664191/1‑141 (MQ=255)                     aaTGTGCAGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTGACATATTTGTTGTCGTCTTTCTTCAGTTAATACTGCACAAGGAGAttctttactatatc                                                                                                                     <  2:543551/141‑10 (MQ=255)                      aTGTGCGGTATCTGCTATTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAAttt                                                                                                                    <  1:572161/141‑1 (MQ=255)                      aTGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAAttt                                                                                                                    <  2:109285/141‑1 (MQ=255)                      aTGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAAttt                                                                                                                    >  2:446059/1‑141 (MQ=255)                        gtgCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAAttttt                                                                                                                  >  1:232201/1‑141 (MQ=255)                         tgCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTg                                                                                                                 <  1:706958/141‑1 (MQ=255)                                  aGTTACTAAAAATTCAATCGTACCTGCACTGCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTAGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCGTCGATTAATTTTTGATTTTTACg                                                                                                        >  2:213054/1‑141 (MQ=255)                                   gTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGt                                                                                                       <  2:75442/141‑1 (MQ=255)                                   gTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGt                                                                                                       <  2:216214/141‑1 (MQ=255)                                     tACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTg                                                                                                     <  2:474961/141‑1 (MQ=255)                                      aCTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACGAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGt                                                                                                    <  1:117092/141‑1 (MQ=255)                                      aCTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACGAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGt                                                                                                    <  1:553096/141‑1 (MQ=255)                                         aaaaaTTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACa                                                                                                 >  1:26533/1‑141 (MQ=255)                                          aaaaTTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAg                                                                                                <  1:612871/141‑1 (MQ=255)                                            aaTTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGaa                                                                                              >  1:265805/1‑141 (MQ=255)                                              ttCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACa                                                                                            <  1:186411/141‑1 (MQ=255)                                                 aaTCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGatcattg                                                                                         >  1:693617/1‑135 (MQ=255)                                                 aaTCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATc                                                                                         <  2:58761/141‑1 (MQ=255)                                                 aaTCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATc                                                                                         >  1:539184/1‑141 (MQ=255)                                                 aaTCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATc                                                                                         >  1:424126/1‑141 (MQ=255)                                                          tGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTcccc                                                                                <  1:627257/141‑1 (MQ=255)                                                          tGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTAGTCTTTCTTCAGTTCATGCTGCACAAGGATATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTcccc                                                                                <  1:242246/141‑1 (MQ=255)                                                           gCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTTTTTTTTTTTTTGTACAGAACAATCACTTTCCCCt                                                                               >  2:434865/1‑141 (MQ=255)                                                                tCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAAt                                                                                          <  2:414784/125‑1 (MQ=255)                                                                tCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAAt                                                                                          >  1:414784/1‑125 (MQ=255)                                                                        ttGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAAt                                                                  >  2:497591/1‑141 (MQ=255)                                                                        ttGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAAt                                                                  >  1:254506/1‑141 (MQ=255)                                                                               aCGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTg                                                           <  2:664191/141‑1 (MQ=255)                                                                                    cAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCAt                                                      <  1:46564/141‑1 (MQ=255)                                                                                     aaCTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCAtc                                                     <  1:446059/141‑1 (MQ=255)                                                                                       cTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCAtctc                                                   >  1:657091/1‑141 (MQ=255)                                                                                              aGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAACTGTACATAATTATCCTGCCCATCTCCCCTAAc                                            <  2:121972/141‑1 (MQ=255)                                                                                                cATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACtt                                          <  1:344567/141‑1 (MQ=255)                                                                                                 aTCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGAtttt                                                                                                            <  1:535926/74‑1 (MQ=255)                                                                                                 aTCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATATTTGAtttt                                                                                                            >  2:535926/1‑74 (MQ=255)                                                                                                      tACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACtt                                    >  2:695861/1‑141 (MQ=255)                                                                                                        cATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCa                                  <  1:90552/141‑1 (MQ=255)                                                                                                         atatTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCaa                                 <  2:653599/141‑1 (MQ=255)                                                                                                             ttgttgTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATg                             <  2:553024/141‑1 (MQ=255)                                                                                                                 tgtCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCGATAACTTGAACTTCAACATGtttt                         >  2:120743/1‑141 (MQ=255)                                                                                                                    cgtcTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCa                      >  1:245837/1‑141 (MQ=255)                                                                                                                    cgtcTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCa                      >  1:292003/1‑141 (MQ=255)                                                                                                                     gtcTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTTTTTTTTACTTTTTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCaa                     >  2:451841/1‑141 (MQ=255)                                                                                                                       ctttcttCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACa                   <  1:497591/141‑1 (MQ=255)                                                                                                                         ttctttAGTTAATGATGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGCACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGGACTTCAACATGTTTTGCAACAgg                 <  2:490105/141‑1 (MQ=255)                                                                                                                         ttcttcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTTTTTTTTTTTTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTTCAACAgg                 >  2:297218/1‑141 (MQ=255)                                                                                                                         ttcttcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAgg                 <  2:535732/141‑1 (MQ=255)                                                                                                                         ttcttcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAgg                 >  1:202402/1‑141 (MQ=255)                                                                                                                           cttcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATc                                                              <  1:145243/94‑1 (MQ=255)                                                                                                                           cttcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATc                                                              >  2:145243/1‑94 (MQ=255)                                                                                                                            ttcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGa                                                                                               >  2:477591/1‑60 (MQ=255)                                                                                                                            ttcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGa                                                                                               <  1:477591/60‑1 (MQ=255)                                                                                                                            ttcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGtat              <  2:3747/141‑1 (MQ=255)                                                                                                                             tcAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGtata             <  2:88893/141‑1 (MQ=255)                                                                                                                              cAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGTATaa            >  2:618784/1‑141 (MQ=255)                                                                                                                              cAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGTATaa            >  2:531481/1‑141 (MQ=255)                                                                                                                                gTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAAt                                                                  <  2:221344/85‑1 (MQ=255)                                                                                                                                gTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAAt                                                                  >  1:221344/1‑85 (MQ=255)                                                                                                                                gTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGTATaaaa          >  1:388872/1‑141 (MQ=255)                                                                                                                                      gctgcACAAGGAGATTCTTCGATTAATTTTTTTTTTTTTCTTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGTATAAAAGCCTCa    >  2:121951/1‑141 (MQ=255)                                                                                                                                        tgcACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGTATAAAAGCCTCAAc  >  2:697721/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                        AGCATTCATTTCAATAAAGTAATGTGCGGTATCAGTTACTAAAAATTCAATCGTACCTGCACTTCTATAATTTGCTGCACGTGCAACTTTAACAGCATCGTTACATATTTGTTGTCGTCTTTCTTCAGTTAATGCTGCACAAGGAGATTCTTCGATTAATTTTTGATTTTTACGTTGTACAGAACAATCACGTTCCCCTAAATGTACATAATTATCCTGCCCATCTCCCATAACTTGAACTTCAACATGTTTTGCAACAGGTATAAAAGCCTCAAC  >  CP000730/1713921‑1714196

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